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BEHAVIORAL  AND  BIOLOGICAL  FFFECTS  OF  RK 503 ANT 
ELECTROMAGNETIC  1’UWF.R  ABSORPTION  IN  RATS 


Abstract 

Past  experiments  of  this  project  utilizing  three  separate  and 
distinct  radiation  facilities  have  determined  the  whole  body  and 
distribution  within  the  body  of  resonant  electromagnetic  power  absorp¬ 
tion  for  both  man  models  and  laboratory  rodents.  The  absorption  of 
such  energy  is  reliably  determined  by  both  the  frequency  of  radiation 
and  the  orientation  of  the  model  or  rodent  in  electromagnetic  energy 
fields.  For  an  ungrounded  1.75  m  tall  man,  power  absorption  has  been 
determined  for  e|]  L,  using  scaled-down  models,  for  the  23  to  570  MHz 
band  with  resonant  power  absorption  at  62-68  MHz.  Power  absorption  for 
the  laboratory  rat  has  been  determined  for  the  300-800  MHz  band  with 
resonant  absorption  being  at  600  MHz  for  the  E ] |  L  orientation.  At 
respective  resonance  frequencies  in  the  E | [  L  orientation  for  both  man 
and  laboratory  rodent,  maximum  power  absorption  is  found  in  the  neck 
region  of  the  body.  For  a  grounded  1.75  m  tall  man,  the  resonance 
region  shifts  to  approximately  one-half  the  free  space  resonance  or 
30-34  MHz.  Man  models  or  laboratory  rats  placed  in  proximity  of  metal¬ 
lic  reflecting  surfaces  experience  considerably  enhanced  electromag¬ 
netic  energy  depositions.  Behavioral  experiments  with  laboratory  rats 
have  shown  that  frequency  of  radiation  and  animal  orientation  in  the 
fields  are  highly  significant  parameters  in  determining  electromagnetic 
energy  absorption  and  consequential  disruption  of  behavioral  performance. 
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Objective 

Ti  e  object!'.  -s  of  the  project  are: 

a.  To  quantity  the  elec tromsgnetic  power  absorption  and  its  dis¬ 
tribution  for  nan  when  subjected  to  radiation  at  different 
frequencies  and  under  various  exposure  conditions. 

b.  To  verify  the  important  observations  of  (a)  by  exposure  of 
living  organisms. 

The  experiments  for  man  are  performed  with  reduced  scale  models 
filled  with  biological-phantom  materials  which  simulate  the  electromag¬ 
netic  properties  of  human  tissue.  The  animal  experiments  consist  of 
behavioral  and  biological  effects  of  electromagnetic  power  absorption  in 
rats. 

The  purpose  of  the  project  is  to  develop  an  understanding  that 
would  lead  to  projections  for  humans. 

Highlights  of  the  Work  Done 
During  the  Last  Year 

1.  The  parallel  plate  radiation  chamber  was  used  to  determine 
the  distribution  of  energy  deposition  for  various  parts  of 
the  body  for  different  sizes  of  biological  phantom  figurines. 
This  work  is  detailed  in  Quarterly  Report  No.  4.  The  high¬ 
lights  of  the  results  are: 

a.  Maximum  rate  of  energy  deposition  is  observed  for  the 
neck  for  E | [  L  resonance  and  near-resonance  conditions. 

For  this  region,  an  energy  deposition  of  nearly  12  times 
the  whole  body  average  is  observed.  The  pattern  of  energy 

-  2  - 

Li 


a  r 


deposition  is  very  sin'lrr  to  mat  obtained  from  fr^;- 

space  Irradiation  experiments .  The  higher  rate  of  energy 

depcs '  v ■•■..*!!  of  2°)  cir.es  the  who]  e  body  average  oi>  -ervnd 

for  the  neck  region  for  frc ;  space  exper ir:en : is  ascribed 

2 

to  the  larger  power  densities  of  100  nV.’/cm  used  there  as 

2 

compared  to  a  maximum  field  intensity  of  26.5  mW/cm  avail¬ 
able  in  the  parallel  plate  radiation  chamber.  The  lower 
field  intensity  needs  a  longer  exposure  time,  which  re¬ 
sults  in  the  heat  being  carried  away  from  the  "hot"  neck 
region  to  the  cooler  torso  and  head  regions  of  the  body. 

b.  For  E||  L  orientation,  a  detailed  examination  of  the  leg 
shows  a  higher  rate  of  energy  deposition  in  the  knee  and 
ankle  regions.  A  power  deposition  varying  as  (cross- 
sectional  area)  ^  is  observed.  A  similar  result  is 
observed  for  the  arm  where  the  maximum  rate  of  deposi¬ 
tion  is  for  the  elbow  region. 

c.  For  lc||  L  orientation  (power  propagating  from  head  to  toe; 

E  from  arm  to  arm) ,  the  maximum  rate  of  energy  deposition 
is  observed  for  the  elbow  region.  The  second  highest 
absorption  rate  is  measured  for  the  neck  region. 

2.  The  monopole-abcve-ground  chamber  was  designed  and  installed 
(Fig.  1).  A  new  computer  program  was  written  to  extend  the 
field  intensity  calculations  to  the  near  field  region.  This 
is  necessitated  by  the  fact  that  high  field  intensities  are 
generally  needed  for  modeling,  animal  heating  and  lethality 
experiments.  This  may  therefore  require  that  the  working 
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region  l>a  located  no  i.-ore  Chan  a  few  wavelengths  from  the 
antenna.  T'.  a  computer  program  calculates  the  vector  sum  of 
near  fie.'  due  the  J  . •  ■  1  •  - ;  and  it.,  image:-;  (in  assumed 
infinite  red J . c tors) .  1:.  order  to  apply  the  results  to  the 

case  of  the  n:..>nopoJ  o-.'bove-ground ,  the  field  components  cal¬ 
culated  for  the  half-space  above  ground  are  used  (because  of 
their  physical  realizability). 

The  field  intensities  calculated  at  various  distances  Z 
from  the  corner  of  the  reflectors  for  some  values  of  the  spac¬ 
ing  d  of  the  monopole  from  the  corner  are  given  in  Table  1. 

Table  1.  Forty-five  degree  corner  reflector 
monopole-above-ground  antenna. 

Power  input  =  100  W 

I  Distance  d  =  33.75  cm 

2 

Ground  plane  field  intensity  in  mW/cm 


Distance  Z 
in  meters 

400  MHz 

d/\  =  0.45 

500  MHz 

0.56 

600  MHz 

0.675 

700  MHz 

0.788 

1.5 

26.0 

26.1 

28.3 

24.7 

2.0 

14.7 

14.7 

15.8 

13.8 

2.5 

9.5 

9.4 

10.1 

8.8 

3.0 

6. 6 

6.5 

7.0 

6.1 

|  4.0 

1 

3.7 

3.7 

3.9 

3.4 

In  each  case  a  transmitted  power  of  100  watts  is  assumed.  A 
point  to  note  is  that  for  distances  beyond  2  meters,  the 
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fields  already  satisfy  the  1 / r ^  dependence  that  is  character¬ 
istic  of  far-1-’ ' .  Ids.  This  is  ascribed  to  the  fact  thac  for 

9 

l  is  ..•'on.ii,  the  2d  >  \  c  ■  lit ion  for  i:  nr-f  i  e '  ds  is  satis- 

2 

fied  for  C .aces  on  the  order  of  2  (id)  /X  or  1. 5-2.0  in  for 
r.ost  of  the  operating  band. 

The  calculated  radiation  patterns  of  the  antenna  in 
the  yz(E-)  and  xz(H-)  planes  are  shown  in  Fig.  2  for  d/X  = 
0.675.  This  value  of  d/X  =  0.675  is  selected  to  obtain 
optimum  input  impedance  and  gain  characteristics2  of  the  radi¬ 
ator.  The  1  dB  beam  widths  in  the  yz  (the  E-plane)  and  xz 
(the  H-plane  or  the  plane  of  the  ground)  planes  are  on  the 
order  of  13.08°  and  ±6.54°,  respectively.  The  targets  oc¬ 
cupying  these  regions  are  consequently  exposed  to  power  den¬ 
sities  varying  by  no  more  than  ±10  percent. 

The  radiation  pattern  of  the  corner  reflector  antenna 
is  relatively  insensitive  to  small  changes  in  d/X  and  monopole 
lengths  h/X.  Table  2  gives  the  beam  widths  for  various  values 
of  d/X  and  h/X.  On  account  of  a  nearly  50  ft  input  impedance 
of  d/X  =  0.675  quarter  wave  monopole  antenna,  this  value  of 
corner  to  antenna  spacing  is  used,  as  far  as  p  -  -  •  i  ' 
convenience  of  matching  over  the  fr  1  ‘!t<  • . 

Two  input  p..rrs  to  w»> :  •  >  may  be  screwed 

.  tei  plates  are  movable  back  and 
J/ •  0.675  at  each  of  the  measurement 


are  i  .•  i 


places  at  a  couple  oC  representative  f  reqne  a  re  o-.:n  ic 

Figs.  3  and  4.  Measurements  confirm  the  similarity  of  pat¬ 
terns  at  other  frequencies. 

3.  Proportionately  scaled  man-shaped  cavities  of  height  3,  4,  5, 
6,  8,  10,  13,  and  16  inches  were  made  in  polyurethane.  Start¬ 
ing  with  wire  frames,  sculptured  figurines  were  made  of  clay 
in  roughly  the  same  proportions  as  those  given  in  reference  3. 
Each  clay  model  was  then  molded  of  Silastic  RTV  silicone  rub¬ 
ber  (Dow  Corning  type  3110).  The  reusable  molds  so  formed 
were  used  to  cast  epoxy  figures.  The  epoxy  figures  were  sub¬ 
sequently  cast  in  polyurethane  blocks.  After  polyurethane 
had  been  formed,  the  epoxy  castings  were  cut  in  half  around 
the  central  frontal  plane  for  their  easy  removal.  Two 
cavities,  one  resealed  for  filling  with  saline  and  the  other 
for  stuffing  with  biological  phantom  materials,  were  formed 
for  each  of  the  eight  sizes. 

The  weights  of  the  epoxy  figures  and  the  volume  of  the 

final  polyurethane  cavities  are  given  in  Table  3.  The  volume 

3 

varies  approximately  as  (height)  ,  which  demonstrates  the  near 
proportionality  of  these  figurines. 

4.  Biological  phantom  materials  were  developed  for  simulating  man 
over  the  frequency  range  13-230  MHz  (including  the  important 

-  9  - 


Measured  field  variation  in  the  E-plane  at  350  MHz. 
d/A  ■  0.675,  Z  =  2.15  m  (fields  normalized  relative 


Measured  field  variation  in 
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Fig.  4  (a).  Measured  field  variation  in  the  E-plane  at  600  MHz.  P  =  50  W,  h/A 
d/A  =*  0.675,  Z  =  2.15  m  (fields  normalized  relative  to  10.6  mW/cm2 


Distance  x  from  the  plane  of  symmetry  in  cms 


Table  3.  bo-.-  measured  parameters  of  the 

turn  models  used  for  experiments. 

Vo’ur.e  of  the 

the  “i-i.:;  t  r  pel /uretha  * 

epux  '  old  cavi  t  •• 

e-'is  err.  •• 


4” 

1 

21.7 

15.0 

5" 

i 

45.7 

33.2 

6" 

59.3 

50.0 

8" 

139.5 

111.4 

10" 

267.5 

208.5 

13" 

548.5 

410.5 

16" 

1063.0 

842.0 

resonance  regions  of  33  and  65  MHz).  These  materials  are  to 

be  used  to  fill  the  eight  proportionately  scaled  figurines  for 

dosimetric  measurements  at  300,  400,  600,  915,  and  985  MHz. 

The  compositions  and  their  measured  properties  were  described 

in  a  paper  presented  at  the  1976  USNC/URSI  meeting,  and  a 

copy  of  the  same  is  attached  here  as  Appendix  A. 

5.  As  a  result  of  our  extensive  experiments  spanning  a  period  of 

almost  two  years,  the  curves  giving  the  whole  body  absorption 

for  man  exposed  to  free  space  electromagnetic  radiation  have 

been  generated  and  are  sketched  in  Figs.  5  and  6.  For  each  of 

u* 

the  indicated  polarizations,  the  orientation4  of  maximum 


fi-;urir.& 
i nches 


A  copy  of  this  paper  (Reference  4)  is  attached  here  as  Appendix  B  for 
ready  reference. 
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E|[  L  orientation 

Direction  of  propagation  from  arm  to  arm 
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Fig.  5.  Whole  body  absorbed  power  density  and  relative  absorption  coefficient  S  fur  human 
electric  field  polarized  along  the  major  length  L  of  the  body  (L’|j  L  orieutnt  •  ->n) . 
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To  date  the  e.\;-  •  rinantal  data  have  been  obtained  to 
8.7  times  the  e||  L  resonance  frequency  of  62-68  MHz.  The 
quantities  plotted  in  Figs.  5  and  6  are: 

a.  The  relative  absorption  coefficient  S  which  is  a  measure 
of  the  efficiency  of  power  absorption 


Microwave  absorption 

_ cross  section _ 

Shadow  area  of  the  body 


(Total  power  absorbed  /(Field  intensity 

_ in  watts)  /  in  watts/cra2) 

2 

Shadow  area  of  the  body  in  cm 


b.  The  specific  absorption  rate  (SAR)  in  watts/kgm. 

c.  The  SAR  as  units  of  basal  metabolic  rate  (BMR) . 

For  the  most  absorbing  E  j ]  L  orientation,  the  whole  body 
absorption  curve  may  be  discussed  in  terms  of  five  frequency 
regions : 

Region  I  —  Frequencies  well  below  resonance  (L/A  <  0.1-0. 2). 

2 

An  f  tyne  depenuvr.ee  derived  theoretically  and  checked 

experimentally  by  barney  and  coworkers/ 

Region  II  —  Subresonance  region  (0.2  <  L/A  <  0.36). 

2.75  3 

An  f  to  f  dependence  of  total  power  deposition  has 

been  experimentally  observed  for  this  region. 

Region  ITT  —  Resonance  region  (L/\  -  0.36-0.4). 
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izations.  For  humans,  this  corresponds  to  (frequency  in 
MHz)  x  (height  in  meters/1.75)  of  about  250. 
b.  The  resonances  for  k||  L  and  h||  L  orientations  are  not 
very  sharp.  In  fact,  for  h]{  L  orientation,  the  value 


gradually  reaches  a  peak  value  and  stays  at  that  value 
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7. 

Also  si  -w a  in  the  same 

figure  are  the  values  calculated  from  the  lethality  data  of 
Schrot  and  Hawkins.7  In  calculating  SAR's  from  their  data, 
it  is  assumed  that: 

a.  A  temperature  elevation  of  7.4°  C  results  in  convulsion. 

b.  The  temperature  increase  occurs  in  an  adiabatic  manner; 
i.e.,  the  rate  of  deposition  is  high  enough  so  that  there 
is  no  other  mechanism  for  gain  or  loss  of  heat  from  the 
animal  body. 

2 

On  account  of  the  high  field  intensities  of  150  mW/cm 
used  in  the  experiments  by  Schrot  and  Hawkins,7  the  second 
assumption  may  be  quite  justified.  From  Fig.  7,  the  correla¬ 
tion  of  the  empirical  graphs  with  experimental  values  is  good. 
7.  The  total  power  absorption  for  a  human  with  feet  touching  the 
ground  is  shown  in  Fig.  8.  These  results  are  projected  from 
measurements  with  4",  5",  6",  and  8"  tall  saline-filled 


In  order  to  determine  the  aspect  ratio  of  the  equivalent  prolate  spher¬ 
oid,  the  bulk  of  the  body  is  considered  representable  by  a  prolate 
spheroid  of  length  L  and  diameter  2b.  The  diameter  2b  is  calculated 
from  the  average  circumference  of  the  animal  body. 
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Fig.  7.  Projected  values  of  whole  body  absorbed  power  density  for  different  size 
without  ground  effects  (L/2b  =  3.25  where  2rb  is  the  "average"  circumfer* 


Fig.  8.  Projected  SAR  in  a  human  of  height  L  (in  meters)  with  feet  touching  ground.  (For 
each  of  the  points  shown  on  the  graph,  three  measurements  were  taken.  The  nver.-r  o 
of  nine  measurements  gives  the  peak  value  for  the  8"  size  at  L/A  =  0.2.) 
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-V  ,  .  A 

figurines  exposed  to  l.,,  I.  radiation  at  different  frequencies 
in  the  nonopole-above-grjund  radiation  chamber.  Comparing  the 
results,  to  those  of  tree  rad  .cion  (Fig.  V),  peak  p.n-  .-r 

ah.-  jrpti.  '  i.i  fills  c-~.  -e  is  o'se".  l\>r  a  f qa  aiv-  if  cn«~ 
half  tV.c-  value  for  ungreu tied  conditions.  The  new  resonance 
frequency  is  projected  to  be  (30-34)  x  (1.75/L  )  MHz.  At 

v  in; 

resonance,  a  value  about  two  times  the  resonance  value  for 
ungrounded  conditions  is  observed. 

A  point  of  concern  is  that  the  graph  of  Fig.  8  does 
not  completely  agree  with  the  data  obtained  with  grounded 
figurines  in  the  parallel  plate  chamber  (Fig.  16,  Appendix  B) . 
While  the  peak  absorption  frequency  is  very  similar,  the  maxi¬ 
mum  absorption  density  observed  in  the  monopole-above-ground 
chamber  is  about  half  as  large  as  that  observed  in  the  paral¬ 
lel  plate  radiation  chamber.  This  is  unlike  the  observations 
with  ungrounded  figurines  where  the  parallel  plate  radiation 
chamber  measurements  were  in  excellent  agreement  with  free 
space  irradiation  experiments.  Further  experiments  are  planned 
to  resolve  the  issue  of  2  versus  4  times  enhancement  factors 
observed  with  grounded  figurines  under  respective  resonant 
conditions. 

8.  Judging  from  the  success  of  the  antenna  theory  in  explaining 
the  resonance  frequencies  and  absorption  cross  sections  under 
grounded  and  ungrounded  conditions,  highly  enhanced  values  of 
SAR  are  expected  for  resonant  bodies  in  the  presence  of  180°, 
90°,  and  otner  corner  reflectors."  Experiments  detailed  in 
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Appendix  C  have  indeed  g i*. en  S'R  enhancement  factors  a.-;  high 
as  27  in  proximity  cf  reflecting,  surfaces  under  resonant 
corn!  i  t  Lo  ■ . 

:  •  ;ul  t  of  r '  ,-e  experi..  )  cv.-  tote.’,  p-.,  vor 

c:  a  .Lie!  •  •  man  '  !  nni.n  •. !  j  to  10  mV!/ incident 

fields  ha  -,  lean  cal.  ulnled  for  various  exposure  conditions 
and  is  t\i>  slated  In  Table  4.  The  rates  of  energy  deposition 
for  resonance  conditions  in  the  presence  of  ground  effects 
and/or  reflectors  are  truly  staggering.  Resonance  SAR  values 
as  high  as  56.5-103  tines  the  BMR  are  predicted  for  adult 
humans.  The  times  to  convulsion  of  100  gm  rats  (Fig.  9, 

2 

Appendix  C)  tor  incident  field  intensities  of  3-20  mW/cm 
confirm  some  of  the  highest  projections  in  the  presence  of 
reflecting  surfaces. 

In  order  to  fully  appreciate  the  reflector-caused 
enhancements  in  SAR's,  it  should  be  mentioned  that,  for  target 
length  L,  reflectors  with  dimensions  2d  x  1.2  L  (in  the  E 
direction)  are  considered  adequate2  to  simulate  infinite 
sheets  for  a  target  to  corner  distance  d.  Furthermore  the 
reflecting  surfaces  need  neither  be  of  good  conducting  materials 
nor  solid  In  construction  to  cause  considerable  enhancements  in 
SAR’s.  In  fact,  surfaces  of  insulating  material  with  conduct¬ 
ing  rods  (oriented  along  incident  fields)  that  are  spaced  <0.1  A 
act  effectively  as  solid  conducting  surfaces.  Experiments, 
nevertheless,  need  to  be  done  with  reduced  dimensions  of 
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1 2 b l e  Total  power  deposited  in  mar.  and  animals  at  10  nW/cm 


Man  (1  .75  rn  height) 


Prior  know! edge 


100  gm  rat 


400  am  rat 


38  v  :  tt$  ass’v::  ro  reflec.  .ee 

19  watts  assuring  50  percent  rerlectan 


I.  At  resonance  for  free  space. 


151  watts 
f  a  62-68  MHz 


0.8  watts 
f  *  987  MHz 


2  watts 
f  «  600  MHz 


II.  At  resonance  for  conditions  of  elec¬ 
trical  contact  with  the  gound  plane. 


2  x  151  =302  watts 
f  *  31-34  MHz 


At  resonance  for  placement  in 
front  of  a  flat  reflector. 


4.7  x  151  =710  watts 
f  «  62-68  MHz 


d  «  0.125  X 


3.8  watts 
f  *  987  MHz 


9.4  watts 
f  =  600  MHz 


Table  4  (continued) 


Man  (1.75  n  height)  100  gin  rat  400  cm  rat 

IV.  At  resonance  roc  S0°  corner  reflector 


27  x  151  =  4077  watts 
f  *  62-68  MHz 


21 .6  watts 
f  =  987  MHz 


54  watts 
f  =  600  MHz 


V.  At  resonance  in  electrical  contact  with 
ground  plane,  in  front  of  a  flat  reflector. 


d  — 


7777 7  /~T?  - 

d  =  0.125  X 


2  x  710  =  1420  watts 
f  =  31-34  MHz 


VI.  At  resonance  in  electrical  contact  with 
ground  plane,  in  a  90°  corner  reflector. 


2  x  4077  =  8154  watts 
f  =  31-34  MHz 
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reflector.-;  and  for  imperfect  reflecting  surfaces  to  determine 
the  enhancement  factors  that  actually  result  in  such  sirua- 


t ions. 


The  re 


e  one ; 


cist-  fr 


run  in  the  mo  no  pole-above-'-.  round  chamLe.. 


ruts  1:.. le-n 

re.-.onance  fre¬ 


quency  for  these  animals  was  pinpointed  in  the  first  instance 
to  be  600  MHz  in  agreement  with  the  graphs  of  Fig.  7.  The 
mean  tiraes-to-convulsion  together  with  the  standard  deviation 
for  four  animals  is  plotted  in  Fig.  9  for  irradiation  fre¬ 
quencies  of  500,  550,  600,  650,  and  700  MHz.  Minimum  time- 
to-convulsion  is  observed  for  irradiation  frequency  of  600 
MHz.  Below  resonance,  the  time-to-convulsion  increases  more 
rapidly  than  that  for  frequencies  above  resonance.  This  ob¬ 
servation  is  consistent  with  the  graphs  in  Figs.  5  and  7,  and 
is  in  agreement  with  the  lethality  data  of  Schrot  and  Hawkins.7 

The  tiraes-to-convulsion  of  400  gm  rats  for  different 
locations  in  front  of  a  90°  corner  reflector  are  shown  in 
Fig.  10.  Of  the  three  locations  that  were  evaluated,  d  = 

0.4  A  is  the  most  absorbing.  For  this  location,  the  times-to- 
convulsion  for  different  incident  field  intensities  are  plotted 


In  actual  experiments,  animals  within  the  weight  range  of  370-430 
grams  were  used. 

On  account  of  the  fairly  limited  length  (=2.5  m)  of  the  monopole- 
above-ground  chamber,  the  more  lethal  placement  d  =  1.5  A  (see 
Fig.  5,  Appendix  C)  could  not  be  tried  under  these  exposure  con¬ 
ditions. 
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convulsion  of  400  gm  rats  for  10  and  20  mW/cm  incident  f  to] 
(Animals  placed  d  =  0.4  X  from  a  90°  corner  of  reflecting 


\  f  A  .1 


in  Fig.  1  !  . 

Cnnipar  !;•;  in  to  is  si".::;ti^n  mean  time— to-conv.Isi^n 
of  730  aecc.  ;  o  -  20  m e...”  iuc  ;er.:  field:,  to  the  co: ve- 

f/pond  L.\  .  t  '  :.e  of  4j'-  seconds  ‘  t  710  v.. .)  ...  rved  by  S  .vt 

.  .  .  .  2  .  . 

ana  nr. '..'nine.  t  *  r-.-a  Space  i:\a_  : :  ,>n  a-:  1..  .•  m.v/cm  ,  it  is 

apparent  that  the  depositin',  enhancement  factor  in  the  mono¬ 
pole-above-ground  chamber  on  the  order  of  7  is  not  as  large  as 
11  observed  for  100  gm  rats  for  d  =  0.4  X  placement.  It 
should  also  be  noted  that  Schrot  and  Hawkins'  experiments7 
were  done  at  710  MHz  and  they  would  have  observed  a  slightly 
lower  convulsion  time  than  455  seconds  had  they  worked  at  the 
resonance  frequency  of  600  MHz  for  these  animals.  As  dis¬ 
cussed  in  the  next  paragraph,  a  lower  enhancement  factor  has 
also  been  observed  in  the  monopole-above-ground  chamber  for 
man  models. 

9.  The  average  values  of  SAR  for  ungrounded  resonant  size  saline- 
filled  figurines  with  and  without  reflectors  have  been  mea¬ 
sured  in  the  monopole-above-ground  chamber  and  are  tabulated 
in  Table  5.  Consistent  with  the  400  gm  rat  lethality  data, 
somewhat  smaller  enhancement  factors  of  2.5  and  6.9,  as  compared 
to  4.7  ar.d  11  observed  under  respective  exposure  conditions  at 
Walter  Reed  Laboratories  (see  Appendix  C,  Figs.  6  and  8),  are 
noted.  Further  theoretical  and  experimental  work  is  needed  to 
understand  this  observation.  The  reduced  effectiveness  of 
in  the  monopcle-above-ground  chambers  is  most 
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Fig.  11.  Time  to  convulsion  of  400  gm  rats  at  their  resonance  frequency  of  600  MHz. 


Table  3.  O. -.on  of  the  whole-h.  Jv-average  "ungrounded" 
SAT.  at  rvsor  mi  e  with  and  without  reflectors. 

Data  scaled  for  1.75  n  h.  an  frcm  the  measurements  on  un¬ 
round -d  8"  «.•  line-filled  figurines  at  600  MHz  (L/X  -  0.406) 

.  ? 

SAR  i  .  ■  .  to..'  f  no  idea',  fields  of  10  iaW/'t 


SARI  _ 

| no  reflectoc 

SAR  C1 
f  lat 

reflector 

SA.f 

v  i  o  j 

| 90  reflector 

1 

CL 

1! 

0.125 

i 

d/A  =  0.4 

2.04 

_ : _ 

5. 

13 

14.02 

probably  due  to  the  spherical  nature  of  the  wavefront  as  it 
impinges  on  the  planar  reflectors.  This  problem  was  not  as 
severe  for  Walter  Reed  experiments  on  account  of  a  larger 
physical  separation  of  14'  of  the  transmitting  antenna  from 
the  experimental  area  which  would  imply  that  the  incident 
waves  were  closer  to  plane  waves  for  those  experiments. 

An  implication  of  the  above  measurements  is  that  while 
higher  enhancement  factors  are  observed  for  plane  waves,  sig¬ 
nificantly  increased  SAR's  are  observed  also  for  nonideal  wave 
incidence  conditions. 

10.  The  average  values  of  SAR  for  grounded  resonant-size  saline- 
filled  figurines  with  and  without  reflectors  have  been  mea¬ 
sured  in  the  monopole-above-ground  chamber  and  are  tabulated 
in  Table  6.  The  results  of  these  experiments  were  used  in 
arriving  at  seme  of  the  numbers  given  in  Table  4. 

strong  absorption  in  the  resonance  region,  our 
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work  to  date  has  concentrated  on  this  region.  Only  recently 
have  we  extended  the  work  to  about  8.7  times  the  resonance 
frequency  which  for  man  corresponds  to  a  frequency  on  the 
order  of  570  MHz.  While  this  has  allowed  us  to  generate 
whole  body  absorption  curves  for  man  (Figs.  5  and  6)  through 
the  important  resonance  region,  we  have  not  addressed  the 
question  of  power  deposition  in  man  at  microwave  frequencies. 
To  extend  the  measurements  on  whole  body  absorption  and  to 
obtain  distribution  for  man  to  3300  MHz,  covering  thereby  some 
fairly  important  microwave  bands  such  as  the  L-  and  S-bands 
and  the  industrial  heating  frequency  of  2450  MHz,  it  will  be 
necessary  to  use  X-band  irradiation  if  reduced  scale  models 
are  to  be  used  for  convenience  of  experimentation. 

High  power  sources  at  X-band  are  not  only  expensive 
but  also  are  not  readily  available  for  power  outputs  larger 
than  250  watts.  A  high  gain  waveguide,  slot  array  has 
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cense  .  irntly  been  designed,  fabclcav  d,  and  tested  for  u.-a 
at  these  frequencies.  The  r  ‘diet  or  descri'  .  ‘  in  App*-.:.  !  .  D 
gi>.  s  an  intensity  that  :  ;  -  ,st.a;U  ,o  within  -0.5  CZ  owv  a 

tori’.  *.i  of  •'  .6°  '  .  the  tad  !  3'T  1-.  the  E-f.me. 

Or.  account  of  a  sharp  reduction  "i  fields  ir  the  H-plar..- 
beyond  the  work  area  (±7.8°),  an  antenna  gain  which  is  2.45 
times  higher  than  that  of  a  horn  with  a  similar  work  area  is 
observed.  This  is  on  account  of  the  saving  in  the  power  that 
would  otherwise  have  been  radiated  in  undesired  directions 

because  of  the  smooth  roll  off  of  the  horn  antenna. 

2 

For  a  working  area  of  55  x  55  cm  ,  a  distance  to  the 

radiator  of  2  meters  is  needed.  From  the  measured  gain,  a 

2 

field  intensity  of  46.4  mW/cm  constant  to  within  ±0.5  dB  is 
calculated  for  a  commercially  available  X-band  power  source 
of  250  W.  This  field  intensity  is  quite  adequate  for  model¬ 
ling  experiments. 

12.  Numerical  studies  have  been  initiated  to  understand  some  of 
the  salient  features  of  our  observations.  Among  these  are 
the  enhancements  in  SAR  caused  by  ground  and  reflecting 
surfaces,  the  1/f  reduction  of  the  absorbed  power  in  the 
post-resonance  region,  etc. 

Green's  function  methods  such  as  the  tensor  integral 
equation  used  by  Chen8  appear  to  hold  the  greatest  promise  for 
calculation  of  SAR  for  models  of  man  at  frequencies  near  whole- 
body  resonance.  Chen  claims  to  have  solved  the  problem  and 
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has  publ resit1. •:  s  or.  the  '.'‘•'■rnal  electric  field  and 
absorbed  p'w.n  demH ty  Inside  an  adult  torso  a  id  a  obi  id’s 
torso  i:  diK.vd  cl.vtiv  \e  •.  s  of  f  i  -:(;ue  .  Les  rang- 

1-  ,  FiTO.'.  :  o>.-  j 

In  order  to  -7 -l  ve  to.-  5  :;t c\.  ral  equation  f  *  ...I  in 
Green1;,  function  methods,  it  is  necessary  to  discretize  by 
representing  the  fields  in  terms  of  some  finite  basis.  For 
a  pulse  function  basis,  the  scatterer  is  partitioned  into  a 
number  of  cells  N,  where  N  is  large  enough  that  complex 
permittivity  and  the  complex  time-independent  electric  field 
may  be  assumed  constant  within  each  subvolume.  Pulse  func¬ 
tions  are  commonly  chosen  as  a  basis  and  have  been  used  in 
Chen's  calculations.  If  linear  or  higher  order  variation  of 
the  electric  field  were  allowed  within  each  cell  by  a  more 
accurate  basis,  the  storage  limitations  of  computers  would 
require  that  fewer  cells  be  used.  The  number  of  cells 
possible  with  pulse  functions  is  adequate  to  represent  the 
torso  with  arms,  legs,  etc.,  but  such  detail  would  not  be 
possible  with  a  more  elaborate  basis. 

Strictly  speaking,  pulse  functions  are  not  in  the 
domain  of  the  operator  and  should  not  be  used  with  the  tensor 
integral  equation  or  related  Green's  function  methods.  In 
practice  usable  results  are  obtained  if  enough  cells  are  used.^ 
Since  pulse  functions  fail  to  allow  for  any  of  the  variation 
of  the  electric  field  within  each  cell,  an  upper  bound  on 
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usable  cell  size  may  be  found  by  determining  the  smallest 
size  guaranteed  to  have  internal  fields  with  significant 
voriat : 

pendix  K  ..  a  -  ;>y  o*  .  pa?  r abmi tt^d  for  ;  cb:  L-:  •- 
tion  i: the  IF  . I'f..-.-  can  t ' -d*  Tk'ory  a-.. I  Tech¬ 

niques.  Theory  is  used  to  establish  an  upper  hound  on  usable 
cell  size  for  pulse  function  solutions.  The  bound  has  been 
demonstrated  to  be  useful  in  our  computations.  Calculations 
made  using  the  bound  on  cell  size  suggest  that  significant 
error  is  present  in  Chen's  calculations  for  the  adult  torso 
beyond  100  MHz  and  the  child  torso  beyond  200  MHz. 

Two  new  discretization  procedures  have  been  developed 
which  give  improved  convergence  by  approximating  the  varia¬ 
tion  of  the  fields  within  each  cell  but  do  not  require 
increased  storage  or  significant  increase  in  computational 
effort.  The  procedures  are  described  in  Appendix  F  which  is 
a  copy  of  a  paper  submitted  for  publication  to  the  IEEE 
Transactions  on  Antennas  and  Propagation.  The  new  methods 
have  been  used  successfully  in  the  two-dimensional  problem  of 
TM  excitation  of  an  infinite  cylinder  of  arbitrary  cross  sec¬ 
tion  shape.  Tha  extension  to  the  three-dimensional  problem 
with  and  without  ground  and  reflector  effects  is  planned. 

To  verify  important  observations  by  animal  experiments, 
attention  was  devoted  to  proper  design  of  the  rat  chambers. 

To  improve  microwave  trar  roarency,  the  animal  holders  were 
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iabrif -1  with  ottom  ar.d  lap  ma  !  •  of  rads  perpendicular 
to  the  tk-ctri:  field  to  ■••'cuce  possible  reflections.  The 
;  rout  nr’  F  .  pl.tt^s  in  ere  p". 'no  of  tk.  einc.ric  fit-id 
t-._re  :■••;’  ■  of  m  .ft. ...  .  '‘rnizc-  .  of  d  -;ign 

:>£  the  c:  nabs.-  was  r._por Co. ...  tori;  hop  .  f  L’o.  5  cl 
this  contract,  and  photographs  of  the  same  are  reproduced 
in  Fig.  12  for  ready  reference.  We  have  found  the  holder 
with  styrofoam  sides  (holder  B)  to  be  adequate  in  contain¬ 
ing  400  gram  Long  Evans  rats  for  short  time  periods.  We 
have  attempted  to  use  this  holder  design  for  the  behavioral 
experiments.  We  find,  however,  that  the  styrofoam  sides  do 
not  provide  adequate  structural  support  for  mounting  a  plastic 
response  lever  and  glass  food  pellet  tube.  To  alleviate  this 
difficulty  we  have  replaced  the  styrofoam  sides  of  this  holder 
design  with  1/16- inch  Plexiglas  panels.  We  have  contained 

400  gram  Long  Evans  rats  in  this  holder  (Fig.  13)  for  up  to 

2 

eight-hour  periods  in  5  mW /cm  ,  600  MHz,  CW  radiation  fields. 
The  addition  of  the  plastic  response  lever  and  glass  food 
pellet  tube  (Fig.  14)  was  easily  accomplished  on  this  holder. 
We  have  used  this  chamber  daily  in  training  stable  VI-30" 
baseline  responding  in  Long  Evans  rats  and  find  this  a 
superior  holder. 

14.  The  first  series  of  dosimetric  measurements  on  the  Long  Evans 
rat  have  been  completed  in  the  monopole-above-ground  radiation 
chamber.  Sixteen  Long  Evans  male  rats  (380-400  grams)  were 
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Holder  A 


Holder  B 

Photographs  of  two  rat  holders  that  should  be  "relatively 
transparent"  to  microwaves  of  Ijj  L  polarization. 


response  leve 


anesthetized  vita  Nembuta]  ->5  rag/kg).  The  LOOK  temperature 


probe1 


placed  7  c:.:  in  t'r 


each  rat.  The  rats 


were  l<  ii»  idually  irrad  -i  once  to  400,  500,  600,  or  7o0 


:;hz,  io 


lei...  .irature  r  ise  of  t* 


:z\  t-  r  ‘  on  ./iaut  :s. 


•>i.-  area  t;as  recorded  using  t i .  a 


LC0F  temperature  probe  end  a  Hewlett-Packard  stripcnart 
recorder.  The  results  of  this  experiment  are  presented  in 
Table  7.  The  average  a  values  indicate  an  enhanced  absorp¬ 
tion  of  energy  at  600  MHz  [o.67  (mW/gm) /( mW/cm2)]  compared 
with  400  [o.  39  (mW/gm) /( mW/cm2  )j  ,  500  fo.46  (mW/gra)  /(.mW/cm2 )j , 
and  700  MHz  [o.50  (mW/gm) / (mW/cm2 )]  .  These  values  of  absorbed 
energy  compare  quite  favorably  with  calculated  power  absorp¬ 


tion  for  a  400  gm  rat  at 


650  MHz  [o.60  (mW/gm) /(mW/cm2)]  . 11 


Colonic  temperature  rise  for  individual  rats  is  shown  in 
Fig.  15.  A  stripchart  recording  of  colonic  temperature  rise 
is  presented  for  each  of  the  frequencies  studied  to  date 
(400,  500,  600,  700  MHz).  Visual  inspection  of  each  record¬ 
ing  shows  that  the  temperature  rise  during  irradiation  seems 

2 

a  linear  function.  Use  of  10  mW/cm  field  intensities  is 
adequate  to  demonstrate  the  resonance  phenomenon  in  Long 
Evans  rats.  We  will,  however,  continue  the  dosimetric  measure¬ 
ments  using  higher  field  intensities  and  whole  body  calorimetry. 
15.  Three  experiments  have  been  completed  to  investigate  thyroid 


(T^)  response  to  resonant  electromagnetic  power  absorption  in 
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Fig.  15.  Rat  colonic  temperature  rise  during  irradiation  at  400, 
500,  600,  and  700  MHz  (l0  mW/cm2). 
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Fig.  15  (continued) 


liu  first  experiment  vjs  designed  to  investigate  the 
effect  of  short  deration  chronic  exposure  of  rats  to  low 

le.  -I  >-.  ion  at  f’.-  •  i  fre"  t  .■  >zy.  For  Long  Evens 

r.i^s  -.It'  gras  !  walf-.t)  •-  ere  c-xp  ■  ~ed  to  5t  MHz  ■ :: 

-v  , 

23.85  mW/ct  r*diac  ion  at  E;,  L  orientation  ir.  the  parallel 
plate  racliat  jn  chamber.  Exposures  'were  15  minutes  duration, 
daily,  five  days  per  week  for  three  weeks.  Colonic  temperature 
measurements  were  made  with  a  BAT-8  digital  electronic 
temperature  probe  before  and  after  each  radiation  exposure. 
Average  colonic  temperature  rise  data  during  irradiation  and 
sham  irradiation  are  presented  in  Table  8. 

Table  8.  Average  weekly  temperature  rise  during  irradiation  and 
sham  irradiation. 
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Animals 
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0.99 

1.26 

0.99 

Standard 
Deviation  °C 

0.51 

0.46 

0.65 

0.46 

i 

0.47 

1 

0.46 

Four  rats  in  the  same  weight  range  underwent  identical  handling 
procedures  and  15-minute  confinement  in  the  radiation  chamber. 
All  animals  were  housed  within  the  radiation  chamber  in  the 
styrofoam  holder  B  shown  in  Fig.  12.  Blood  samples  were  drawn 


from  all  animals  under  ether  anesthesia  by  tail  snip  at  the 
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c  •  12 
sis  of  variance. 

The  results  of 

this  test 

indicate  no  sig- 

nificant  difference  between  the  control  and  irradiated  animals 
(P  >  .10).  There  was  a  difference  across  weeks  (P  <  .01) 
which  was  significant.  This  difference  was  most  likely  due  to 
the  low  measured  values  of  on  week  2.  Since  the  values  for 
control  animals  were  also  low,  it  seems  evident  that  this 
change  is  due  not  to  microwave  exposure  but  perhaps  to  some 
other  factor  such  as  stress  during  the  time  the  exposures  were 
carried  out. 

In  the  above  study,  blood  samples  were  drawn  2  to  4 
hours  after  radiation  exposure  and  under  ether  anesthesia. 

Two  additional  studies  were  conducted  to: 

a.  Withdraw  blood  both  as  soon  as  possible  and  at 
extended  times  after  a  single  radiation  exposure. 

b.  Collect  blood  samples  by  a  method  least  stressful 
to  the  animal. 

Toward  this  goal  twenty-six  long  Evans  rats  (380-430  gins  body 
weight)  had,  under  Nembutal  anesthesia,  in-dwelling  catheters 


Fig.  16.  T,  levels  of  control  and  UHF  radiated  rats 
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I'l.th  8-  ■  veil.  r-  •  ch  :  ••  it  began  3  to  5  d‘;s 

surgery. 

To  determine  the  effect  of  resonance  frequency  radiation 
on  thyroid  function  (T^  release),  seventeen  animals  were 
divided  into  five  groups.  A  blood  sample  (.4  ml)  was  withdrawn 
from  each  of  the  seventeen  animals  and  was  considered  a  base¬ 
line  value.  Each  group  of  animals  was  then  exposed  once  to 

600  MHz  radiation  in  E||  L  orientation  at  100  (N  =  3),  75  (N  = 

2 

4),  50  (N  =  4),  and  25  (N  =  2)  mW/cm  ,  respectively,  for  five 
minutes.  The  fifth  group  underwent  sham  irradiation  procedures. 
All  radiation  exposures  were  done  in  the  monopole-above-ground 
radiation  chamber  immediately  after  the  withdrawal  of  blood 
for  a  baseline  value.  During  irradiation  and  sham  irradiation 
within  the  monopole-above-ground  chamber,  rats  were  housed  in 
the  Plexiglas  holder  (Fig.  13)  described  in  this  report.  Ad¬ 
ditional  samples  were  ti..  ■’vn  from  the  animals  at  15  and  90 
minutes  after  the  initial  baseline  sample  was  drawn.  Average 
body  temperature  rise  during  radiation  and  sham  irradiation  is 
presented  in  Table  9.  All  blood  samples  were  analyzed  for  T^ 
by  Pantex  radioimmuno  assay  of  thyroxine.  The  results  of  this 
analysis  are  presented  in  Fig.  17.  Examination  of  these  data 
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Fig.  17.  Thyroxine  (T^)  response  to  600  MHz  radiation, 
deviations  ar ■  not  shown  overlapping.) 
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again  reveals  little  effect  of  low  or  high  intensity  resonant 

radiation  exposure  on  thyroxine  (T^)  levels. 

To  determine  any  long  term  effects  of  short  duration 

irradiation,  the  remaining  nine  animals  were  divided  into  two 

groups.  Five  animals  were  exposed  to  600  MHz  radiation  for 

2 

five  minutes  at  100  mW/cm  and  four  animals  underwent  sham 
irradiation  procedures  as  in  the  above  study.  Average  body 
temperature  rise  during  radiation  and  sham  irradiation  is 
presented  in  Table  10. 


Table  10.  Average  colonic  temperature  rise 

during  radiation  and  sham  irradiation. 
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of  Ihyro ... i.ne.  The  results  of  th  L  ■>  an  ilysis.  Fig.  18,  again 

reveal  little  effect  of  resonant  radiation  exposure  on  T, 

4 

response  in  the  rat  as  compared  to  sham  irradiated  animals. 

In  the  three  studies  described  above,  there  is  a 

gradual  decrease  in  with  time.  This  effect,  however, 

appears  for  both  irradiated  and  sham  irradiated  animals. 

Animals  of  the  first  experiment  subjected  to  repeated  ex- 
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posure  of  low  level  radiation  (23.85  mW/cm  )  did  appear  very 
docile  and  lethargic  upon  handling  compared  to  the  sham  ir¬ 
radiated  controls.  This  may  have  been  due  to  the  daily 
thermal  load  to  which  the  animals  were  subjected.  This  ef¬ 
fect  was  not  observed  for  animals  of  the  second  and  third 
experiments.  Nevertheless,  the  gradual  decrease  in  T^  ap¬ 
pears  to  be  due  not  to  microwave  exposure  but  rather  to  some 
other  factor.  Although  the  evidence  is  by  no  means  conclusive, 
it  seems  apparent,  under  the  experimental  conditions  used, 
that  no  direct  effects  of  resonance  frequency  irradiation  on 
thyroid  function  were  observed  in  the  adult  rat. 

16.  Experiments  have  recently  been  completed  to  determine  the  ef¬ 
fect  of  resonant  and  nonresonant  electromagnetic  power  absorp¬ 
tion  on  rats  performing  a  lever-pressing  task  on  a  variable 
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of  the  animal  parallel  to  the  electric  fields.  All  radiation 

exposures  produced  behavioral  deficits.  The  suppression  of 

behavior  was  greatest  for  exposures  to  600  MHz  radiation.  The 

six  remaining  rats  were  exposed  to  600  MHz  CW  radiation  at 

2 

field  intensities  of  5,  7.5,  10,  and  20  mW/cm  while  perform- 

2 

ing  the  behavioral  task.  The  exposures  at  20  and  10  mW/cm 

produced  suppression  of  behavior,  while  exposures  at  5  and 
2 

7.5  mW/cm  showed  no  measurable  effect  on  performance.  In 
addition  some  rats  were  exposed  to  pulsed  600  mHz  radiation 
(1000  pps  at  3  and  30  microsecond  pulse  durations).  Pulsed 
radiation  conditions  showed  no  measurable  effect  on  rat 
performance.  The  data  show  that  the  resonant  frequency  of 
power  absorption  for  the  adult  Long  Evans  rat  is  600  MHz, 
while  suppression  of  VI-30  second  responding,  due  to  short¬ 
term  radiation  exposure,  requires  a  field  intensity  of 
2 

10  mW/cm  or  greater.  A  copy  of  this  paper  (submitted  to 
Radio  Sri-.n^e  for  editorial  review)  is  attached  as  Appendix  G. 
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BIOLOGICAL  PHANTOM  MATERIALS  FOR  SIMULATING  MAN 
AT  DIFFERENT  FREQUENCIES 

Om  P.  Gandhi  and  Kazem  Sedigh 
Departments  of  Electrical  Engineering  and  Bioengineering 
University  of  Utah,  Salt  Lake  City,  Utah 


Abstract 


The  paper  gives  compositions  of  biological  phantom  materials 
for  simulating  man  over  the  frequency  range  13-230  MHz  (including 
the  important  resonance  regions  of  3  and  65  MHz).  These  materials 
have  been  developed  to  fill  eight  proportionately  scaled  man-shaped 
cavities  of  height  7.6-40.6  cm  for  dosimetric  measurements  at  300, 

400,  600,  915,  and  985  MHz.  Complex  permittivity  (ec)  measurements  at 
the  irradiation  frequencies  were  made  using  a  modified  version  of  the 
coaxial  line  method  used  previously  by  Guy.  Models  reduced  by  B  in 
all  dimensions  allow  simulation  of  man  at  frequencies  given  by 
(experimentation  frequency) /B,  provided  materials  with  ec  at  the 
(lower)  simulated  frequencies  are  used.  Whole-body  "average"  values 
of  ec  at  simulated  frequencies  were  first  calculated  on  the.  basis  of 
65  percent  muscle  and  tissues  of  high  water  content  and  35  percent 
fat,  bone,  and  tissues  of  low  water  content  from  extrapolated  values 
for  these  tissues  from  Johnson  and  Guy's  article.  Several  compositions 
of  salt,  polyethylene  powder,  Superstuff,  and  water  were  first  measured 
for  ec's  at  the  six  experimental  frequencies.  Interpolated  composi¬ 
tions  were  then  measured  and  modified,  if  necessary,  for  desired  per¬ 
mittivities  . 
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BIOLOGICAL  PHANTOM  MATERIALS  FOR  SIMULATING 


MAN  AT  DIFFERENT  FREQUENCIES 

Ora  y.  Gandhi  and  Kazem  Sedigh 
Departments  of  Electrical  Engineering  and  Bioengineering 
University  of  Utah,  Salt  Lake  City,  Utah 


Introduction 

An'  important  aspect  of  electromagnetic  wave  biological-effects 
research  is  the  quantification  of  absorbed  dose  and  its  distribution 
in  biological  systems.  For  humane  the  experimentation  may  be  done  by 
using  either  the  full  scale  or  reduced  scale  models  that  use  materials 
having  the  same  complex  permittivity  ec  =  -  jo/we^  as  that  of  the 

system  that  is  to  be  simulated. 

For  whole  body  exposure  situations,  reduced  scale  models  are 
preferred  for  convenience  of  handling.  From  electromagnetic  field 
theory,  a  body  reduced  by  a  factor  of  3  in  all  dimensions  may  be  used 
to  obtain  absorption  characteristics  of  the  full  body,  provided  an 
irradiation  frequency  scaled  up  by  a  factor  cf  B  if  used.  For  dis¬ 
tribution  of  power  deposition,  high  field  intensities  on  the  order  of 
2 

50-150  mW/cm  are  generally  needed  to  prevent  the  pattern  from  being 
smeared  by  beat  diffusion.  A  limitation  often  encountered  is  that  the 
high  power  sources  needed  for  experimentation  are  not  very  widely 
alterable  in  their  operation  frequencies.  Use  of  carefully  propor¬ 
tioned  different  size  models  at  a  given  experimental  frequency  may 
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therefore  be  used  to  simulate  man  at  frequencies  given  by  (experi¬ 
mental  frequency)/^.  It  is  necessary,  however,  that  for  such  experi¬ 
ments  the  complex  permittivity  for  the  mixtures  used  in  the  various 
models  be  different  and  have  the  values  corresponding  to  the  values 
for  man  at  the  reduced  frequency  that  is  being  simulated.  With  these 
scaling  precautions,  the  distribution  of  power  deposition  in  the 
reduced-scale  model  is  identical  (though  the  magnitudes  are  higher 
by  a  factor  of  f$)  to  that  of  the  full-scale  body. 

To  simulate  lower  frequencies  where  wavelengths  are  consider¬ 
ably  larger  than  the  dimensions  of  individual  parts  of  the  body,  whole 
body  "average'1  values  may  be  used  for  the  e^'s  of  the  homogeneous 
filler  media.  This  is  a  valid  procedure  to  frequencies  on  the  order 
of  200-300  MHz.  At  higher  frequencies  a  more  detailed  representation 
of  the  body  is,  however ,  needed.  Different  compositions  of  biological- 
phantom  materials  should  then  be  used  to  simulate  the  "average"  elec¬ 
trical  properties  of  individual  parts  of  the  body.  It  is  fortuitous 

* 

that  at  such  frequencies  the  depth  of  penetration  in  the  torso  is  on 
the  order  of  a  few  centimeters,*  and  only  the  compositions  of  the  top 
layers  may  therefore  be  considered  to  calculate  the  "average"  e^. 

For  legs,  arms,  head,  etc.,  a  more  careful  determination  of  may  be 
in  order,  based  in  the  first  instance  on  the  distribution  of  power 
deposition  as  though  these  parts  were  isolated  from  the  rest  of  the 

While  detailed  power  distributions  in  the  torso  are  not  known  and 
this  assumption  may  indeed  need  to  be  justified,  the  previous 
experiments  with  spherical  targets*  also  showed  that  the  penetra¬ 
tion  of  waves  at  larger  frequencies  [X  <  diameter  d]  was  super¬ 
ficial  in  nature. 


For  near-field  or  part-body  exposure  situations,  full  scale 
models  (0  -  1)  are  quite  convenient  and  may  therefore,  be  preferred. 
Experiments  of  this  type  are  often  needed  in  arriving  at  the  optimum 
design  for  applicators. 

Materials  and  Methods-  of  Procedure 

As  described  above  it  is  necessary  to  develop  biological-phantom 
mixtures  for  complex  permittivities  at  the  experimental  frequencies 
that  simulate  man  at  lower  frequencies.  For  available  experimental 
frequencies  of  300,  400,  600,  915,  and  985  MHz,  the  e^’s  needed  for 
eight  proportionately  scaled  man-models  have  been  calculated  and  are 
given  in  Table  1  for  an  illustrative  case  of  600  MHz.  Corresponding 
numbers  have  been  calculated  for  the  remaining  four  frequencies  and 
require  e^'s  varying  between  35  and  105  and  a  varying  between  2.3  and 
14  mho/m. 

Using  the  homogeneous  compositions  developed  here,  it  should  be 
possible  to  obtain  dosimetric  measurements  for  man  over  the  frequency 
region  13-230  MHz,  which  includes  the  important  resonance  regions2*3 
of  33  and  65  MHz. 

Nine  mixtures  of  compositions  given  in  Table  2  were  measured  at 
the  five  experimental  frequencies,  using  the  arrangement1*  of  Fig.  1. 

A  slotted  coaxial  line  (General  Radio  type  874  Lb)  was  filled  with 
the  respective  mixtures  and  used  to  measure  the  attenuation  constants 
a  and  propagation  constants  P  through  the  media  with  the  help  of  the 


Calculation  of  the  electrical  propertie 
needed  of  the  biological  phantom  materi 
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Table  2.  The  compositions  of  the  nine  mixtures  used  for 
measurements  of  dielectric  properties. 


Mixture 

Percent  by  Weight 

h2o 

NaCl 

* 

ss 

mk-k 

PEP 

I 

75 

13.0 

12.0 

0 

II 

75 

11.1 

7.9 

6.0 

III 

75 

9.8 

8.2 

7.0 

IV 

75 

9.0 

8.0 

8.0 

V 

75 

5.45 

10.55 

9.0 

VI 

75 

4.45 

10.55 

10.0 

VII 

75 

4.0 

10.5 

10.5 

VIII 

75 

3.9 

10.1 

11.0 

IX 

75 

2.0 

10.0 

13.0 

Superstuff  —  A  gelling  agent,  obtained  from  Wharao  Manufacturing 
Company,  San  Gabriel,  California. 

** 


Polyethylene  pov/der  —  obtained  from  Wedco,  Inc.,  California. 


Fig.  1.  Block  diagram  of  the  setup  for  coaxial  slotted  line  measurement  of  complex  permittivity. 


Hewlett-Packard  network  analyzer  model  8410A.  Dae  to  the  high  loss  of 
the  materials,  measurements  were  taken  at  0.2  cm  intervals  and  l\  to  6 


readings  were  used  for  the  calculations  of  a  and  8.  In  order  to  find 
and  o  of  the  materials,  the  following  equations5  were  used: 


u 


-J-  (G  +  1) 


m  |fr 
“  “  c  1  2 


(G  -  1) 


(1) 

(2) 


g 

where  c  is  the  velocity  of  light  (3  x  10  meter/second)  and 


G  = 


1  +  (a/6 )' 
1  -  (a/8): 


1  -I- 


\  we  e  / 

\  o  r / 


(3) 


Upon  rearranging  the  terras. 


p  _  2(«/(w/c)) 
r  G  -  1 


(4) 


and 


(c2 


a/2 

1)  We  e 


o  r 


„  2,,  ,  ,2  G  -  1 
2a  /(w/c)  u>eo  ' 


(5) 


The  values  of  e  and  o,  for  the  nine  mixtures,  measured  at  the  five 
r 


experimental  frequencies  arc  given  in  Tables  3  and  4.  The  values  of 


For  these  mixtures  with  high  cr/we0 Er,  G  is  very  much  higher  than  unity 
which  results  in  a  and  8  being  comparable  in  magnitude  to  one  an¬ 
other.  Small  errors  in  determination  of  either  a  or  8  results  in  a 
much  larger  error  in  G  and  consequently  Er-  This  method,  though 
reasonably  accurate  for  determination  of  a,  does  not  appear  to  be 
very  good  for  er  for  high  loss  (o/t or  e  >>  1)  materials. 


Table  3.  Values  of  the  dielectric  constant  and  conductivity  of 
the  various  mixtures  measured  at  frequencies  300,  400 
and  600  KHz. 


Mixture 

300  MHz 

400  MHz 

600  MHz 

E1 

a 

£1 

a 

E1 

o 

1 

63.52 

14.09 

71.6 

12.73 

51.80 

15.70 

II 

34.25 

12.05 

48.72 

11.74 

54.65 

12.90 

III 

21.35 

11.73 

67.89 

9.28 

53.50 

12.61 

IV 

46.47 

8.89 

71.52 

8.45 

46.77 

9.98 

V 

39.45 

5.03 

45.36 

5.94 

51.90 

6.26 

VI 

64.37 

5.06 

58.79 

4.20 

55.99 

5.44 

VII 

41.62 

5.1 

33.74 

5.62 

59.91 

4.84 

VIII 

39.09 

4.02 

18.82 

4.68 

55.87 

4.04 

IX 

57.82 

2.59 

62.11 

f 

2.25 

38.64 

3.29 
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33.88  12.78  251.18  44.38  13.62  248.66 


II  52.56  12.98  255.11  57.06  12.16  222.01 


III  49.52  9.74  191.43  51.57  9.95  181.66 


IV  43.12  9.22  181.21  47.25  9.05  165.41 


48.36  6.01  118.12  48.56  6.82  124.52 


VI  51.91  5.27  103.58  50.73  5.32  97.13 


VII  62.4  4.91  96.50  63.53  5.13  93.66 


VIII  52.08  4.28  84.12  48.44  4.13  75.40 


IX  52.89  2.68  52.67  27.05  1.49  27.20 


T 


conductivity  at  different  Rl’  frequencies  are  plotted  in  Fig.  2  for 
different  salt  contents  and  compared  with  the  dc6  and  low  frequency1* 

(f  -  0.1  MHz)  values.  The  conductivity  increases  linearly  with  the 
percentage  of  salt  and  is  relatively  independent  of  the  frequency. 

For  any  required  conductivity,  therefore,  a  percentage  of  salt  needed 
for  the  mixture  can  therefore  be  estimated  from  Fig.  2.  To  date  we 
have  not  developed  any  simple  relationship  for  e^'s  needed  for  the 
simulating  media. 

For  simulation  at  lower  frequencies  (including  the  important 
resonance  regions),  o/(oeo  is  2  to  4  times  e  (see  e.g.  Table  4)  and 
therefore  dominates  the  complex  permittivity  of  the  medium.  Matching 
o/toeo  accurately  with  the  required  values  while  allowing  a  slight  mis¬ 
match  in  therefore  does  not  cause  great  errors  in  dosimetric  measure¬ 
ments  for  simulated  frequencies  less  than  about  200  MHz.  This  point 
has  been  confirmed  by  measuring  distribution  of  absorbed  energy  using 
the  liquid  crystal  fiber  optic  temperature  probe7  for  near-resonant 
condition  (L/A  -  0.417)  for  free-space  irradiation.  The  model  was 

•ft 

filled  with  a  biological  phantom  material  simulant  of  man  at  71.6  MHz. 
The  distribution  pattern  scaled  to  values  for  the  full  scale  body  is 

A 

For  L/A  =  0.417  or  a  5"  model  at  985  MHz  to  simulate  1.75  meter  man 
(at  71.6  MHz),  it  is  necessary  that  the  complex  permittivity  er  - 
ja/u)E0  of  the  filler  medium  be  58.4  -  jl30.6,  computed  on  the  basis 
of  65  percent  muscle,  skin,  and  tissues  with  high  water  content 
and  35  percent  fat,  bone,  and  tissues  with  low  water  content.  The 
composition  used  for  free-space  irradiation  experiments  (5.75  per¬ 
cent  NaCl,  10.25  percent  Superstuff,  9  percent  polyethylene  powder, 
and  75  percent  water)  has  a  measured  permittivity  of  49.5  -  jl31.7. 


Nad  percent 

Fig.  2.  Electrical  conductivity  of  phantom  materials  and  saline  versus  percentage  of  salt 
for  different  frequencies. 


p 


f 


shown  in  Fig.  3.  For  free-space  irradiation  specific  absorption  rates 
(SAR)  considerably  higher  than  the  whole-body  average  are  observed  for 
the  neck,  the  legs,  and  the  elbows,  with  the  lower  torso  receiving 
SAR's  comparable  to  the  average  value  and  the  upper  torso  receiving 
SAR’s  lower  than  the  average  value.  Noting  that  the  values  indicated 
in  Fig-  3  are  relative  to  whole-body  saline-filled  average,  it  is  re¬ 
markable  that  the  weighted  average  of  the  numbers  obtained  with  the 
biological  phantom  mixture  is  not  far  off  from  one,  or  the  value  ob¬ 
tained  v/ith  the  saline-filled  figurine. 

Conclusions 

In  the  foregoing  we  have  given  the  measured  properties  of  some 
biological-phantom  mixtures  for  experimental  frequencies  to  985  MHz. 
From  these  results  it  is  possible  to  obtain  the  salt  content  for  a 
conductivity  required  of  the  biological-phantom  material.  Work  is 
continuing  with  mixtures  of  different  water  and  Superstuff  contents  to 
develop  a  ternary  diagram  that  may  be  of  help  in  obtaining  required 
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ABSTRACT 


Experimental  results  are  described  for  whole-body  absorption  (for  orienta¬ 
tions  £||  L,  £||  L,  and  l5||  L)  of  a  series  of  saline-filled  and  biological- 
phantom  figurines  of  major  lengths  (L)  that  varied  from  0.4  to  1.5  X.  Measure¬ 
ments  are  also  reported  for  several  exposure  conditions  on  the  density  of  energy 
deposition  for  different  regions  of  the  body.  In  the  absence  of  ground  ef¬ 
fects ,  a  whole-body  absorption  cross-section  as  high  as  4.2  times  the  shadow- 
cross  section  has  been  measured  for  E||  L  at  the  resonance  frequency  (in  KHz)  of 
68  x  (1.75/height  1^  of  the  human  in  meters).  Under  these  frequency  and 
polarization  conditions,  a  whole-body-average  energy  deposition  rate  of  (2. IS) 

*  (1 . 75/Ljjj)  watts/kgm  is  projected  for  incident  plane  waves  of  10  mW/cm^.  1  he 
highest  rnte  of  energy  deposition  is  observed  for  the  neck  region  of  the  body 


where  a  value  as  high  as  28.7  times  the  whole  body  average  has  been  measured 
under  conditions  of  F||  f,  resonance.  Maximum  whole-body  absorption  under 
grounded  conditions  is  observed  at  a  frequency  about  one  half  that  of  the  un¬ 
grounded  body.  For  grounded  resonance  condition  an  absorption  cross  section  as 
high  as  14  times  the  shadow  area  has  been  measured.  Here,  while  the  neck  is 
still  one  of  the  hot  spots,  the  zone  of  maximum  power  deposition  is  the  ankle 
region. 


1.  INTRODUCTION 

We  have  previously  [1,  2]  reported  on  the  observation  of  strong  resonance 
of  whole-body  absorption  of  electromagnetic  waves  in  biological  bodies.  In 
the  absence  of  ground  effects,  highest  rate  of  energy  deposition  occurs  for 
fields  polarized  along  the  longest  dimension  of  the  body  for  frequencies  such 
that  the  major  length  is  on  the  order  of  0.4  times  free  space  wavelength  (X)  of 
radiation.  At  resonance  an  effective  absorption  area  of  about  3-4  times  the 
shadow  cross  section  has  been  measured  using  saline-filled  prolate  spheroidal 
"equivalents"  of  man.  In  these  prolate  spheroids,  a  major  to  minor  axes  ratio 
of  6  was  used.  On  the  basis  of  some  preliminary  experiments,  peak  absorption, 
in  the  presence  of  ground  effects,  has  been  observed  for  frequencies  one  half 


Also  consultant  to  the  Division  of  Neuropsychiatry,  Walter  Reed  Army 
Institute  of  Research,  Washington,  D.  C.  20012. 
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as  much  as  for  bodies  isolated  in  free  space. 

In  this  paper  we  report  further  results  on  this  important  phenomenon.  In 
particular,  data  are  presented  on  whole-body  absorption  of  salina-f illed 
figurines  of  major  lengths  L  varying  from  0.3  to  1.5  X.  Measurements  are  also 
reported  on  the  distribution  of  power  deposition  in  biological-phantom  figu¬ 
rines  exposed  to  free-space  irradiation  at  both  985  and  2450  MHz  and  to  fields 
in  the  TOM  parallel  plate  [2]  chamber.  The  results  obtained  under  the  latter 
irradiation  conditions  are  shown  to  be  remarkably  similar  to  those  with  free- 
space  radiation.  The  highest  intensity  of  energy  deposition  in  both  chambers 
is  observed  for  the  neck  region  of  the  body.  The  a  parameter  (defined  as  the 
ratio  of  watts/kgia  of  absorbed  energy  to  raW/cm?  Df  in;  ide.it  field  intensity) 
for  the  neck  region  is  approximately  30  tines  that  observed  fur  the  whole-body 
average.  Initial  results  are  given  on  figurines  with  feet  touching  the  ground 
plane  of  the  parallel-plate  chamber. 


2.  MATERIALS  AND  METHODS 

Cavities  corresponding  to  human  shapes  were  formed  in  styrofoam  by  using 
12.1,  18.4,  22.3,  and  23.5  cm  tall  dolls  as  wolds.  After  styrofoam  had  been 
formed,  the  dolls  were  cut  in  half  around  the  sagittal  plana  for  easy  removal 
of  the  molds.  Three  cavities,  one  resealed  for  filling  with  0.9  percent  saline 
and  two  for  stuffing  with  biological-phantom  materials,  were  formed  for  each  of 
the  four  sizes. 

Different  compositions  of  the  biological-phantom  material  were  used  to 
fill  the  figurines  for  parallel-plate  and  free-space  irradiation  experiments. 
From  our  previous  results  [1],  the  frequency  of  maximum  absorption  for  a  1.75 
meter  tall  human  is  on  the  order  of  63--70  MHz  (L/A  -  0.37-0.4).  In  order  to 
model  human  for  this  frequency  band,  an  "average"  electrical  relative  permit¬ 
tivity  er  was  firs!  calculated  on  the  basis  of  65  percent  muscle,  skin,  and 
tissues  with  high  water  content  and  35  percent  fat,  bone,  and  tissues  with  low 
water  content.  The  value  of  er  needed  for  modeling  for  the  63-70  MHz  band  was 
calculated  [3]  to  be  60.5  -  jl43.  The  use  of  whole-body  "average”  er  is  justi¬ 
fied  because  the  dimensions  cf  individual,  constituents  are  so  much  smaller  than 
the  wavelength  of  radiation  for  this  baud.  For  use  at  500  MHz  in  the  parallel- 
plate  chamber,  the  composition  of  the  biological-phantom  material  was  3.26  per¬ 
cent  Na Cl,  8.74  perc-nt  Superstuff  (obtained  from  Whamo  Manufacturing  Company, 
San  Gabriel,  California),  and  87.0  percent  water.  This  has  a  measured  [4] 
dielectric  constant  of  66  and  a  conductivity  a  of  4,39  mho/meter.  At  the 
experimental  frequency  of  500  KHz,  the  complex  relative  p- "mi ttivity  of  the 
material  is  theref'  tr  =-  66  —  j  o/(Jt0  «  66  -  J157.9.  I  different  composition 

of  8.03  NaCI,  13.21,  rcent  Superstuff,  and  78.74  percent  water  was  used  for 

free-space  irradiation  experiments  at  985  KHz.  For  use  at  this  frequency  er  *=> 
66  -  j!64.3,  once  again  close  to  the  needed  value  of  60.5  -  jl,43. 

The  irradiation  facilities  used  were  the  parallel  plate  TEM  chamber  at 
S  '.I  fv >z  and  free-space  anechoic  chambers  at  985  and  2450  MHz.  For  the  parallel 
I-  if  ch-.a.ber  [2],  the  working  area  consists  of  a  copper  plate  of  width  63.5  cm 
;  >ted  by  a  25.4  cm  clearance  from  the  ground  plane  of  116.8  cm  width.  The 
'  lull  transmission  length  of  198.2  cm  is  occupied  by  two  symmetrical  tapered 
'■  '  of  Cl  cm  axial  length,  the  central  working  area  cf  61  cm  length,  and 

•  (•  c.i  id  sections  connected  to  UC58/U  input  and  output  coaxial  connectors. 

W,  200  to  500  MHz  F.psco  signal  source  was  used  to  feed  pow-r  to  this 
'  i  .  Tl..-  fields  in  the  central  (working)  region  of  the  chc  '•  »r  were  mapped 
*  t  ••  General  Microwave  Corporation  Raham  model  1  (300-18, Ci>0  MHz)  field 
i'  •/  probe.  The  variation  of  the  field  intensity  in  the  working  region  is 
’  in  Figs.  1(a),  (b),  and  (c) .  The  985  and  2450  MHz  free  space  fields 

:  r-ivided  by  Scientific  Atlanta  pyr,  fdal  horns  radiating  into  the 
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rcPi’t-c t iv •'  ancchoic  chambers  at  the  Department  of  Microwave  Research,  Waller 
Kecd  Army  Institute  of  Research ,  Washington,  I).  C. 

The  whole  body  power  absorption  was  obtained  from  the  temperature  rise  of 
saline-filled  figurines.  The  temperature  before  and  after  irradiation  was 
measured  with  a  digital  thermometer  after  shaking  the  styrofoam  cavity  figu¬ 
rines  in  order  to  equalize  the  temperature  over  the  volume  of  the  saline.  To 
determine  the  distribution  of  power  deposition,  biological-phantom  figurines 
were  exposed  to  the  highest  available  radiation  intensities  to  speed  up  the 
heating  process  and  thus  reduce  the  heat  diffusion  to  cooler  areas  of  the  body. 
Temperature  under  irradiation  was  recorded  using  the  liquid  crystal  temperature 
probe  [5]. 

The  incident  field  intensities  used  in  the  experiments  were  considerably 
higher  for  free-space  exposures  on  account  of  a  2000  W  power  source  available 
at  Walter  Reed  Laboratories.  The  field  intensities  used  for  free-space  and 
parallel-plate  irradiations  were  100  and  26.5  mW/cm^,  respectively. 

From  the  measured  increase  in  temperature  AT,  the  density  of  power  absorp¬ 
tion  in  W/kgm  is  calculated  from  the  expression  (41.80  AT  x  specific  heat  of  the 
medium  in  calories/gm/°C) /irradiation  time  in  seconds.  This  divided  by  the 
field  intensity  in  mW/cm^  gives  the.  a  parameter  either  for  the  whole  body  (for 
saline-filled  figurines)  or  for  different  parts  of  the  body  (for  biological- 
phantom  figurines).  The  values  of  the  a.  parameter  so  calculated  are  then  re¬ 
duced  by  the  model  scaling  factor  (height  of  the  human/length  of  the  figurine) 
to  obtain  the  corresponding  values  for  humans. 

3.  RESULTS 

Free-Space  Irradiation  Experiments 

The  whole  body  average  value  of  a  defined  in  W/kgm  of  absorbed  power 
density  divided  by  incident  field  intensity  in  mW/cm^  is  calculated  (by  scaling) 
from  the  measured  results^at  98^  and  at  2450  MHz  and  is  plotted  in  Fig.  2.  For 
each  of  the  orientations  £||  L,  k||  L,  ^nd  nj]  ,  two  distinct  exposure  configura¬ 
tions  are  possible.  (The  vectors  E,  H,  and  k  are  along  the  electric  and  mag¬ 
netic  fields  and  along  the  direction  of  propagation,  respectively;  L  is  along 
the  major  length  L  of  the  body.)  For  electric  field  along  the  major  length, 
for  example,  these  are: 

1.  Power  propagating  from  front  to  back. 

2.  Power  propagating  from  arm  to  arm. 

A  5  to  10  percent  larger  vhol e-body  ab  rption  is  found  for  case  2  for 
«lectric  polarization..  The  results  for  thin  most  absorbing  configuration  are 
tlso  plotted  in  Fig.  2.  Table  1  summarizes  the  projected  a  for  the  six  distinct 
orientations  under  respective  conditions  of  maximum ^absorpt ion.  The  most  dif¬ 
ference  in  absorption  is  found  for  k-polnrization  (k||  1.)  where  a  50  percent 
Increase  in  overall  absorption  is  measured  for  electric  field  from  arm  to  arm 
®s  compared  to  the  case  where  the  electric  field  is  from  front  to  back  of  the 
b"dy.  The  results  observed  with  saline-filled  figurines  are  qualitatively 
* t-'- .iar  to  those  observed  with  rectangular  parallelepipeds  [6]  pi  biological- 
pnsntou  material.  At  (63)  x  (1.75/l.m)  MHz  which  is  the  frequency  of  maximum 
i'di  body  power  deposition  (see  Fig..  2),  the  a  x  (I,  /l. 75)  values  for  k|(  I.  and 
I-  orientations  are  0.017  and  0.011,  respectively, W«md  both  of  these  values 
*rc  core  than  an  order  of  magnitude  smaller  than  the  value  for  electric  polar- 
nation  of  0.216  (W/kgm)/ (r.iW/cn^) .  From  a  parameter  the  relative  absorption 
c°dficient  S  (defined  as  the  electromagnetic  absorption  cross  section  divided 


Projected  rate  of  energy  deposition  in  a  human  of  height  meters 
(scaled  froa  measurements  performed  on  saline-filled  doll9  under  free- 
space  irradiation  conditions  at  985  and  2400  MHz).  The  standard 
deviation  of  three  measurements  for  each  of  the  data  points  is  less 
than  3  percent. 


electromagnetic  absorption  in  a  human  of  height  L  meters 


is  agrees  with  our  previously  reported  condition  of  kL  ■  2t rL/X  *  a/b.  For  man  a/b  »  6-6.5,  giving  therefore  an 


r 


by  the  shadow  area  of  the  body)  is  calculated  from  the  expression: 

_  ^  ^ _ wcip.ht  of  the  body 

shadow  area  of  the  body 

The  value  of  the  S  coefficient  thus  calculated  for  e||  f.  and  "h||  I.  orientations 
is  also  shown  in  Fig.  2.  At  resonance  an  effective  absorption  area  of  4.2 
times  the  shadow  cross  section  is  measured.  For  the  post-resonance  region  to 
frequencies  on  the  order  of  3  to  4  times  the  resonance  frequency  fr,  p  whole- 
body  absorption  reducing  as  (f/fr)-^  from  the  resonance  value  may  be  seen  from 
Fig.  2. 


The  measured  rates  of  temperature  increase  are  used  to  calculate  a's  for 
various  parts  of  the  body  in  a  human  for  different  exposure  parameters.  The 
values  of  a  multiplied  by  the  field  intensity  in  mW/ctu^  give  the  rates  of 
deposition  in  W/kgm.  Figures  3-5  give  the  rates  of  energy  deposition  for 
10  mW/crn^  incident  fields  for  prescribed  exposure  conditions.  While  Figs.  3-5 
give  an  idea  of  the  rates  of  energy  deposition,  the  temperature  rise  for  the 
various  parts  in  an  actual  human  would,  of  course,  depend  upon  several  complex 
factors  such  as  blood  circulation,  convective  cooling,  etc. 


Experiments  in  the  Parallel -PI ate  Radiation  Chamber 


The  TEM  radiation  chamber  was  used  to  measure  the  whole  body  a  parameter 
of  the  four  sizes  of  the  saline-filled  figurines  at  the  highest  power  source 
frequency  of  500  MHz.  The  a  parameters  so  obtained  were  scaled  to  the^values 
for  humans  and  these  are  plotted  in  Fig.  6  as  a  function  of  h/X.  For  F.||  I, 
orientation  the  results  are  noteworthy  in  two  respects: 

1.  Maximum  power  absorption  is  observed  for  L/X  -  0.37.  This  is  In 
agreement  with  the  condition  for  maximum  absorption  observed  with 
free-space  experiments  (Fig.  2).  This  result  is  surprising  in  that 
the  parallel  plate  chamber  seeuis  to  simulate  the  free-space  plane 
waves  rather  than  the  ground  plane  effects.  As  will  be  pointed  out 
later,  physical  contact  with  the  ground  plate  does  indeed  allow  this 
chamber  to  simulate  ground  effects. 


2.  The  value  of  a  at  resonance  observed  with  the  parallel  plate  chamber 
is  about  80  percent  of  the  value  obtained  from  free-space  irradiation 
experiments.  In  view  of  the  fact  that  different  frequencies,  figurine 
sizes,  and  calibrating  field  probes  were  used  for  the  two  radiation 
conditions,  the  two  results  should  be  considered  in  remarkably  good 
correlation. 


While  we  have  always  been  convinced  of  the  validity  of  using  a  parallel 
plate  transmission  line  to  provide  a  medium  of  plane  waves,  and  the  above  re¬ 
sults  are  a  further  confirmation  of  our  feelings,  doubts  have,  however,  been 
expressed  by  some  researchers  in  the  field  regarding  this  transmission  line. 
Also  a  handicap  of  the  parallel  plate  line  is  the  rather  limited  space  avail¬ 
able  between  the  two  plates,  and  this  hinders  its  application  to  frequencies 
where  the  target  size  may  be  comparable  to  or  larger  than  one  half  the  free- 
space  wavelength.  On  account  of  these  reasons,  we  have  designed  [2]  and  in¬ 
stalled  a  monopole-above-ground  radiation  chamber.  The  ancchoic  chamber  it- 
formed  by  lining  a  10’  x  10'  x  7.5'  Faraday  screen  room  with  sheets  of  micro- 
wave  absorbing  material  of  pyramidal  (Eccosorb  VHP-18)  and  planar  (F.ccor.orb 
AN-79)  varieties.  The  radiator  consists  of  a  quarter-wave  mor.opole  above 
ground  mounted  in  a  45°  corner  reflector.  Due  to  the  high  gain  of  the  antenna, 
field  intensities  on  the  order  of  13  ir.W/cm^  have  been  measured  at  2.4  meters 
from  the  radiator  for  100  watts  input  to  the  antenna.  An  MCL  model  11145 
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Pav  -  2.15  •  (1.75/Lro) 
(saline-filled) 
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FREE-SPACE  IRRADIATION 

FiC.  3.  Power  deposition  in  W/kgm  for  1.73  n  human  under  F.  jj  f.  resonance 

condition  (I./X  “  0.4,  67.9  •  (1.75/1^)  MHz).  Incident  fields:  10  mW/ 

cm^. 


For  a  human  of  height  L  meters,  the  numbers  marked  alongside  the 
figure  should  be  multiplied  by  (1.75/L^). 
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Pav  .  0.67  •  (1.75/L„> 
(sal ine-filled) 


FREE-SPACE  IRRADIATION 


Power  deposition  in  W/kgra  for  1.75  n  human  (LA  -  0.985,  169  •  (1.75/ 
Ljj)  MHz).  Incident  fields:  10  mW/cra^. 

For  a  human  of  height  i^,  meters,  the  numbers  marked  alongside  the 
figure  should  be  multiplied  by  (1.75/L-,). 


(saline-filled) 


Power  deposition  in  V/kgm  for  1.75  m  man  for  £|]  1.  resonance  condition 
(L/X  «  1,  169  •  (1.75/1^)  MHz).  Incident  fields:  10  mW/cm^. 

For  a  human  of  height  1^,  meters,  the  numbers  marked  alongside  the 
figure  should  be  multiplied  by  (1.75/L.,). 
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1000  W  output  40C-809  MHz  RF  amplifier  ir.  presently  on  order,  and  this  would  be 
used  in  the  continuing  dosimetry  experiments.  While  detailed  data  from  this 
chamber  are  not  available  at  this  time,  preliminary  experiments  have  given  u 
values  fairly  comparable  to  those  of  Fig.  6. 

Figures  7-11  give  the  rates  of  energy  deposition  in  various  parts  of  the 
body  determined  from  the  parallel  plate  chamber  for  prescribed  exposure  condi¬ 
tions.  At  least  three  measurements  w.ue  taken  to  obtain  the  values  indicated 
on  these  diagrams.  A  standard  deviation  of  less  than  5  percent  for  hotter 
regions  of  the  body  and  a  value  of  less  than  10  percent  for  the  cooler  regions 
were  calculated  from  the  experimental  data.  From  Figs.  7-11  the  important 
features  to  note  are: 

1.  Maximum  rate  of  energy  deposition  once  again  is  observed  for  the  neck 
for  K ||  L  resonance  and  near-resonance  conditions.  For  this  region  an 
energy  deposition  of  nearly  10  times  the  whole  body  average  is  ob¬ 
served.  This  ratio  is  considerably  lower  than  28.7  observed  with 
free-space  irradiation  (Fig.  3)  on  account  of  '  w  field  intensities 
used  to  date  with  parallel  plate  chamber  experiments.  This  would 
cause  a  lower  power  deposition  rate  resulting  in  heat  being  carried 
away  from  the  "hot"  neck  region  to  the  cooler  torso  and  head  regions 
of  the  body.  With  the  future  availability  of  higher  field  intensities 
(~150  mW/cti2),  it  is  felt  that  energy  deposition  rates  comparable  to 
those  observed  under  free-space  irradiation  conditions  may  indeed  be 
obtained  with  parallel  plate  TEM  chambers. 

2.  For  E|]  L  orientation,  a  more  detailed  examination  of  the  heating  of 
leg  shows  a  higher  rate  of  energy  deposition  in  the  knee  and  ankle 
regions  (Figs.  8  and  9).  These  also  are  the  zones  of  smaller  cross- 
sectional  areas  than  the  rest  of  the  leg.  The  plot  of  a  as  a  function 
of  the  cross-sectional  area  over  the  entire  length  of  the  leg  is  shown 
in  Fig.  12.  A  rate  of  energy  deposition  varying  inversely  as  the 
cross-sectional  area  is  observed.  If  a  constant  current  were  flowing 
through  the  entire  leg,  a  (cross-sectional  area)-2  type  dependence  of 
a  would  be  expected.  However,  even  from  antenna  theory,  a  diminishing 
value  of  current  is  anticipated  as  one  approaches  the  ankle  region. 
This  may  well  be  the  reason  for  the  experimentally  observed  (cross- 
sectional  area)-*-  dependence  of  a. 

3.  A  detailed  examination  of  the  heating  pattern  of  the  anu  (Figs.  8  an  1 
9)  shows  the  maximum  rate  of  deposition  in  the  elbow  region.  It 
should  be  mentioned  that  the  arms  are  in  the  plane  of  the  body  in  our 
figurines.  In  Figs.  3,  4,  and  9,  the  arms  are  shown  protruding  only 
for  clarity  of  marking  the  measured  values. 

4.  For  k ||  L  orientation  (Figs.  4  and  11),  maximum  rates  of  energy  deposi¬ 
tion  are  observed  for  the  elbow  followed  by  the  neck  region. 


Power  Deposition  in  the  Presence  of  Ground  Effects 

Experiments  were  performed  with  the  18.4  cm  saline-filled  figurine  ground¬ 
ed  to  the  bottom  plate  of  the  parallel  plate  chamber  by  means  of  a  thin  wire 
probe  forming  the  electrical  connection  through  a  rubber  cork  at  the  base  of 
the  figurine.  On  account  of  an  increased  rate  of  energy  deposition  under  these 
conditions,  only  two  minutes  of  irradiation  time  (as  against  five  minutes  in 
experiments  with  ungrounded  figurines)  was  found  to  be  adequate  for  these 
exposures.  The  values  of  a  calculated  from  the  measurements  over  the  frequency 
band  220-500  MHz  are  plotted  in  Fig.  13.  Also  shown  in  the  same  figure  are  the 
values  projected  for  humans  of  different  lirights  and  the  relative  absorption 
coefficient  S.  Peak  power  absorption  is  observed  for  a  frequency  slightly 
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PARALLEL-PLATE  RADIATION  CHAMBER 


Fig.  7.  Power  deposition  in  W/kgra  for  1.75  ni  human  (L/X  ~  0.37,  63.5  •  (1.75/ 
L^)  MHz).  Incident  fields:  10  roW/cm^. 

For  a  human  of  height  Ljj,  meters,  the  numbers  marked  alongside  the 
figure  should  be  multiplied  by  (1.75/L™). 
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(saline-filled) 


PARALLEL-PLATE  RADIATION  CHAMBER 


Pig.  8.  Power  deposition  in  W/kgra  for  1.75  meter  human  (L/A  "  0.31,  52.6  • 
(1.75/Lpj)  KHz).  Incident  fields:  10  mV/cn^. 

For  a  human  of  height  meters,  the  numbers  marked  alongside  the 
figure  should  be  multiplied  by  (1.75/I.m). 
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PARALLEL-PLATE  RADIATION  CHAMBER 


Power  deposition  In  W/kgra  for  1.75  m  human  (L/X  =  0.31,  52.6  •  (1.75/ 
Lp,)  MHz).  incident  fields:  10  mW/cm^. 

For  a  human  of  height  meters,  the  numbers  marked  alongside  the 
figure  should  bn  multiplied  by  (1.75/b-). 
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Pay  *  0.38  •  (1 . 75/Lni) 
(saline-filled) 


PARALLEL-PLATE  RADIATION  CHAMBER 


Fig.  10.  Power  deposition  in  W/kgra  for  1.75  m  Human  (I./ A  **  0.2,  34,6  •  (1.75/ 
1^)  MHz).  Incident  fields:  10  roW/cra^. 

For  a  human  of  height  1^,  meters,  the  numbers  marked  alongside  the 
figure  should  be  multiplied  by  (1.73/1^). 
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Fig.  11.  Power  deposition  in  W/kgm  for  1.75  m  man  (1,/A  -  0.39,  67.1  •  (1.75/ 
Ljj)  MHz).  Incident  fields:  10  mW/cra  . 

For  a  human  of  height  Ljj,  meters,  the  numbers  marked  alongside  the 
figure  should  be  multiplied  by  (1.75/1.-,). 
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*  Value  for  various 
regions  of  the  leg 

o  Value  for  various 
regions  of  the  arm 


Cross-sectional  area  A  of  the  figurine  in  cm 


Fig-  12.  Plot  of  a  as  a  function  of  the  cross-sectional  area  for  arms  and  legs 
(18.4  cm  figurine;  500  MHz,  L/X  *=■  0.31). 
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the  data  roinca  is  shown. 


lover  than  one  half  the  value  for  ungrounded  conditions.  At  resonance  an  ef¬ 
fective  absorption  area  of  nearly  14  times  the  shadow  cross  section  is  measured 
which  is  about  A  times  the  value  measured  at  resonance  under  ungrounded  condi¬ 
tions  in  the  parallel  plate  chamber.  Both  these  results  are  predicted  approxi¬ 
mately  by  antenna  theory.  The  derivations  giving  an  S  of  A. 2  for  free-space 
irradiation  and  a  four-times  enhancement  in  the  absorption  cross  section  of  an 
isolated  versus  a  grounded  body  are  illustrated  in  Figs.  1A  and  15. 

Judging  from  the  success  of  the  antenna  theory  in  explaining  most  of  the 
whole  body  absorption  observations  under  grounded  and  ungrounded  conditions, 
highly  enhanced  values  of  S  may  result  for  bodies  in  the  presence  of  180°,  90°, 
and  other  corner  reflectors.  Preliminary  experiments  in  the  monopole-above¬ 
ground  chamber  have  shown  this  to  be  true  indeed.  S  parameters  on  the  order  of 
20  end  higher  have-  been  observed  to  date.  Further  experiments  are  in  progress 
to  completely  understand  these  effects. 

The  distribution  of  power  deposition  under  grounded  resonance  conditions 
has  been  determined  using  the  biological-phantom  figurine  with  feet  touching 
the  ground  plate  of  the  parallel  plate  chamber.  The  rates  of  energy  deposition 
projected  for  a  human  exposed  to  10  mW/cm^  fields  are  6hown  in  Fig.  16.  The 
points  to  note  are: 

1.  The  whole  body  average  value  is  fairly  close  to  that  obtained  with 
saline-filled  figurine  grounded  by  means  of  a  wire  probe,  lending 
some  credence  to  this  method  of  grounding  for  saline-filled  figurines. 

2.  While  the  neck  is  still  one  of  the  hot  spots,  the  region  of  maximum 
power  deposition  shifts  to  the  lower  half  of  the  leg  under  grounded 
conditions. 

3.  Consistent  with  the  previous  observations  (Figs.  8  and  9)  on  the 
detailed  distribution  in  the  leg  region,  the  zones  of  maximum  power 
deposition  are  the  ankle  and  the  knee,  which  also  correspond  to  the 
areas  of  minimum  cross  section  in  the  leg. 

4.  CONCLUSIONS 

In  the  abeenc.e  of  ground  effects ,  a  whole  body  absorption  cross  section  as 
high  as  A. 2  times  the  shadow  cross  section  has  been  measured  for  electric 
polarization  at  the  resonance  frequency  (in  MHz)  of  (68)  x  (1.75/belght  L,,,  of 
the  human  in  meters).  Under  these  frequency  and  polarization  conditions,  a 
whole-body-average  energy  deposition  rate  of  (2.15)  x  (1,75/1^)  watts/kgra  is 
projected  for  incident  plane  wave  fields  of  10  rnW/cm^.  At  higher  frequencies 
(to  f  =  3-4  fr),  the  whole  body  absorption  reduces  as  (f/fr)-I  from  the 
resonance  value.  The  distribution  of  power  deposition  is  given  for  the  various 
regions  of  the  body  for  several  exposure  conditions.  For  E||  f.  resonance  condi¬ 
tions,  the  highest  rate  of  energy  deposition  is  observed  for  the  neck  region  of 
the  body,  and  a  value  as  high  as  28.7  times  the  whole  body  average  is  measured. 

Maximum  whole  body  absorption  under  grounded  conditions  is  observed  at  a 
frequency  about  one  half  that  of  the  ungrounded  body.  For  grounded  resonance 
condition  an  absorption  cross  section  as  high  as  1A  times  the  shadow  area  has 
been  measured.  Antenna  theory  has  been  used  to  explain  this  highly  enhanced 
absorption  cross  section.  The  distribution  of  power  deposition  is  determined 
for  grounded  resonance  condition.  While  the  neck  is  still  one  of  the  hot  spots, 
the  zone  of  maximum  power  deposition  is  the -ankle  region.  Judging  from  the 
success  of  antenna  theory  in  explaining  most  of  the  whole  body  absorption  ob¬ 
servations  under  grounded  and  ungrounded  condition,  highly  enhanced  deposition 
rates  are  anticipated  for  bodies  in  the  presence  of  180°,  90°  ,  and  other 
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A  biological  body  isolated  in  free  space 
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For  humans  a/b  =  6.34  which  gives  an  S  of  4.21  in  close  agreement  with 
experimentally  observed  values  for  human-shaped  figurines  under  free-sp;.ce 
irradiation. 


Fig.  14.  Derivation  of  the  S  parameter  of  a  biological  bcdy  Isolated  in  free 
space. 
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Body  of  shadow  area  A  in 
free  space. 

Resonance  condition  L  =  0.5  A 

EH  absorption  cross  section 
"  4.2  A 


Two  bodies,  each  of  shadow 
area  A,  placed  end  to  end. 

Resonance  condition  L  =  0.25  A 

EM  absorption  cross  section 
-  4.2  x  2A 


One  half  body  placed  above  ground. 

Resonance  condition  L  -  0.5  A 

EM  absorption  cross  section  =  8.4  A 

(From  antenna  theory,  the.  gain  of  a 
tnonopole  above  ground  is  twice  as 
much  as  that  of  a  twice-length  dipole.) 


A  body  of  shadow  area  A  above  ground. 

Resonance  condition  L  =  0.25  A 

EM  absorption  cross  section  = 

8.4  x  2A 


A  four-times  enhancement  at  resonance  for  conditions  I  and  IV  may 
be  seen  while  the  resonance  frequency  is  reduced  to  one  half  the  value. 

Fig.  15.  Comparison  of  the  S  parameter  for  a  body  in  free  space  with  the  vale 
for  the  same  body  in  the  presence  of  a  ground  plane. 
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Fig.  16. 


Pav  =  0.7  (1.75/Lm) 
(saline-filled) 


Energy  deposition  in  mW/gm  for  1.75  m  tall  human  under  resonance 
condition  (T./X  “  0.15,  25.2  MHz)  for  feet  touching  ground;  incident 
power:  10  mW/cra^. 

For  a  human  of  height  1^  meters,  the  numbers  marked  alongside  the 
figure  should  be  multiplied  by  (1.75/Lg,)* 
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corner  reflectors. 


5.  REFERENCES 

1.  Gandhi,  0.  P.  Polarisation  and  frequency  effects  on  whole  animal  absorp¬ 
tion  of  RF  energy.  Proc.  IEEE*  62,  1171-1175,  1974. 

2.  Gandhi,  0.  P.  Conditions  of  strongest  electromagnetic  power  deposition  in 
man  and  animals.  IEEE  Trans.  KTT,  KTT-23,  1021-1029,  1975. 

3.  Johnson,  C.  C.,  and  A.  W.  Guy.  Nonionizing  electromagnetic  wave  effects  in 
biological  materials  and  systems.  Proc.  IEEE,  60,  246-253,  1972. 

4.  Guy,  A.  U. ,  C.  C.  Johnson,  J.  C.  Lin,  A.  F.  Emery,  and  K.  K.  Kranirg. 
Electromagnetic  power  deposition  in  man  exposed  to  KF  fields  and  the  associated 
thermal  and  physiologic  consequences.  USAF  School  of  Aerospace  Medicine, 

Brooks  Air  Force  Base,  Texas,  Report  SAM-TR-73-13,  December  1973. 

5.  Rozzell,  T.  C.,  C.  C.  Johnson,  C.  H.  Durney,  J.  L.  Lords,  and  R.  C.  Olsen. 

A  nonperturbing  temperature  sensor  for  measurements  in  electromagnetic  fields. 
J.  Microwave  Power,  9,  241-249,  1974. 

6.  Gandhi,  0.  P.  Orientational  and  frequency  effects  on  whole  animal  absorp¬ 
tion  of  RF  energy.  J.  Indian  Inst,  of  Sci.,  56,  315-364,  1974. 


ACKNOWLEDGMENT 

The  technical  support  of  J.  H.  Jacobi  and  P.  B.  Brown,  Department  of 
Microwave  Research,  Walter  Reed  Arr;>  Institute  of  Research,  is  gratefully 
acknowledged.  This  work  was  supported  by  the  U.  S.  Army  Medical  Research  and 
Development  Command  under  Contract  DAMD  17-7A-C--4092. 


r 


A. 


Radio  Soienae ,  Vol,  12,  No.  6  (S) 
November-December  1977. 
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Deposition  of  electromagnetic  energy  in  animals  and  in  models  of  man  with  and  without  grounding  and 

reflector  effects 
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Generalized  curves  are  given  tor  rales  of  whole-body  absorption  of  electromagnetic  energy  by 
models  of  human  beings  as  a  function  of  frequency.  Exposures  were  made  in  free  space,  or  when  the 
model  made  electrical  contact  with  the  ground.  The  effects  of  the  presence  of  reflecting  surfaces  also 
we  e  analyzed.  Peaks  of  absorption  with  and  without  a  ground,  are  projected,  respectively,  to  be 
(31  to  34  MHz)  and  (62  to  68  MHz)  X  (1.75/height  of  model  in  meters).  Rates  of  energy  deposition 
are  given  for  models  of  man  and  for  animals  subjected  to  radiation  at  a  power  density  of  10  mW/cm* 
for  the  various  conditions  of  exposure.  At  resonance,  values  of  whole-hodv  absorption  as  high  as 
4,077  to  8,154  watts  for  the  adult  human  being  are  predicted.  The  times-to-convulsion  of  ~  100-g 
rats  at  power  densities  of  3  to  20  mW/ctn2  confirmed  predictions  of  extremely  high  rates  of  absorp¬ 
tion  in  the  presence  of  reflecting  surfaces. 


I.  INTRODUCTION 

We  have  previously  reported  on  the  observation  of  a 
strong  resonance  in  the  whole-body  absorption  of  elec¬ 
tromagnetic  waves  by  biological  bodies  [ Gandhi ,  1974, 
1975|.  For  plane-wave  irradiation  in  free  space,  the 
highest  rate  of  energy  deposition  occurs  in  fields  that  arc 
polarized  along  the  longest  dimension  of  the  body 
(E||L)  and  for  frequencies  such  that  the  major  length 
(/.)  is  approximately  0..36  to  0.40  times  the  frec-spacc 
wavelength  (X)  of  radiation.  Peaks  of  whole-body  ab¬ 
sorption  for  two  other  configurations  (major  length 
oriented  along  the  direction  of  wave  propagation  or 
along  the  vector  of  the  magnetic  field)  were  also  re¬ 
ported  [ Gandhi .  1974,  1 9 7 5 1  for  l./\  on  the  order  of 
l./4tth  where  2tth  is  the  weighted  averaged  circumference 
of  the  animals. 

A  handicap  of  past  measurements  has  been  that  com¬ 
mercially  available  dolls  that  have  been  used  as  molds 
for  models  are  not  properly  proportioned.  Using  care¬ 
fully  sculptured,  accurately  scaled  figurines  [Drcvfuss, 

1 967 1 ,  measurements  have  now  been  extended  to  cover 
the  region  0.2  <  /./X  <  .3.3.  Specific  absorption  rates 
(SARs)  are  given  in  this  paper  under  thice  configura¬ 
tions  in  free  space  for  the  frequency  range  0.5  to  8.7 
times  the  E||L  resonance  frequency  Jr.  which,  for 
man.  is  65  X  ( 1 ,75/lcngtli  m  meters)  MHz.  From  these 
results,  antenna  theory  is  used  to  develop  an  empirical 
equation  for  the  SAR.  foi  the  E||L  orientation  of 
models,  at  resonant  and  suprarcsonant  frequencies.  A 
comparison  with  the  lethality  data  of  Schrot  and 


Hawkins  [  1975)  confirms  the  validity  of  the  equation 
for  rats  and  mice  of  differing  sizes. 

Results  are  also  presented  for  a  simulated  human  being 
with  its  feet  in  electrical  contact  with  the  ground.  In 
such  an  exposure,  a  resonant  frequency  one-half  as  high 
with  a  peak  absorption  twice  that  of  free-space  irradia¬ 
tion  is  observed. 

Energy  deposition  rates  are  given  for  models  of  man 
and  for  animals  that  were  exposed  to  radiation  at  power 
densities  of  10  mW/cm2  under  various  conditions  of 
exposure. 


2.  ELECTROMAGNETIC  ENERGY  ABSORPTION  IN 
MAN  IN  FREE  SPACE 

The  procedure  for  obtaining  the  mass-normalized  rate 
of  electromagnetic  energy  deposition  (SAR)  through  the 
use  ol  reduced-scale  models  has  been  detailed  in  earlier 
publications  [Gandhi.  1975;  Gandhi  ft  ai.  1975].  The 
SAR  in  W/kg  is  calculated  from  the  temperature  rise, 
T.  of  7.6.  10.2.  12.7.  15.2.  20.3,  25.4.  3.3.0,  and  40.6 
cm  high,  saline-filled  (0.9 %  NaCI)  figurines  from  the 
expression  (4.180  T  X  specific  heat  of  the  medium 
in  calorics/g/r'C (/irradiation  time  in  seconds.  The  SARs 
so  calculated  are  divided  by  the  model's  scaling  factor 
(height  of  the  human/length  of  the  figurine)  to  obtain 
the  value  for  a  human  being.  At  least  three  measure¬ 
ments  were  taken  to  calculate  the  averaged  SARs  that 
are  given  in  this  paper.  Unless  otherwise  stated,  the 
standard  deviation  of  the  measurements  is  less  than 
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±  3%.  The  values  obtained  witli  saline-filled  figurines 
were  found  to  agree  with  the  average  of  the  values  for 
individual  parts  of  the  body  for  tissue-simulating 
figurines,  wlich  lends  some  credence  to  this  method  of 
arriving  at  the  whole-body  absorption  curves  for  man. 
A  more  accurate  procedure  that  is  yet  to  be  carried  out 
is  calorimetric  measurement  of  whole-body  energy  depo¬ 
sition  in  figures  that  are  filled  with  materials  that  dupli¬ 
cate  the  differing  conductive  and  dielectric  properties  of 
various  tissues. 

Curves  for  whole-body  absorption  by  models  of  man 
exposed  to  radiation  in  free  space  are  given  in  Figures 
1  and  2.  For  each  of  the  indicated  polarizations,  the 
orientation  resulting  in  a  maximum  rate  of  energy 
deposition  was  used  [Gandhi  et  ah,  1975],  These 
orientations  correspond  to  propagation  normal  to  the 
sagittal  plane  ~  from  arm  to  arm  (rather  than  from  front 
to  back)  -  for  E  1 1 L  orientation;  E  normal  to  the  sagittal 
plane  (rather  than  from  front  to  back)  for  k  1 1  C  orienta¬ 
tion;  and  E  normal  to  the  sagittal  plane  for  H  |  (  C  orien¬ 
tation. 

To  date,  the  experimental  data  have  been  obtained  to 
8.7  times  the  E||C  resonant  frequency  which  lies 
between  62  to  68  MHz.  The  quantities  plotted  in 
Figures  1  and  2  are: 

1.  The  relative  absorption  coefficient  S,  which  is  a 
measure  of  the  efficiency  of  energy  absorption: 

Microwave  absorption 
,  _  cross  section 

Physical  cross 
section  of  the  body 

(Whole-body  rate  of  energy  /  (Field  intensity 
absorption  in  watts)  /  in  watts/cm2) 

Physical  cross  section  of  the  body  in  cm2 

2.  The  specific  absorption  rate  (SAR)  in  W/kg. 

3.  The  SAR  as  compared  to  basal  metabolic  rate 
BMR). 

In  calculating  the  SAR-to-BMR  ratios  for  Figures  1 
and  2,  a  basal  metabolic  rate  [ Guyton ,  1968]  of  100 
kilocalories/hr  is  assumed  which  corresponds  to  1 .66 
W/kg  for  a  70-kg  man.  A  higher  metabolic  rate  (MR)  of 
1 .66  to  1 .96  W/kg  is  assumed  to  occur  during  normal 
activity  such  as  standing  while  relaxed,  dressing  and 
undressing,  etc.;  use  of  any  number  in  this  range  would 
not  appreciably  alter  the  SAR/MR  ratios  shown  in 
Figures  1  and  2. 

For  the  most  highly  absorbing  E||L  orientation,  the 
whole  body  absorption  curve  may  be  discussed  in  terms 
of  five  regions  of  frequency: 

Region  1  --  Frequencies  well  below  resonance  (7,/X  < 
0.1  -0.2).  An  /2-typc  dependence  derived  theoretically 
and  checked  experimentally  by  Durney  and  coworkers 
[Massoudi  et  at.,  1975] . 

Region  2  --  Subresonant  regio  i  (0.2  <  L/\  <  0.36). 
An  /2'7S  to  /3  dependence  of  rate  of  energy  deposition 
has  been  experimentally  observed  for  this  region. 

Region  3  -  Resonant  region  (/./X  —  0.36-0.4).  A 
relative  absorption  cross  section  Srcs  on  the  order  of 


0.665  L/2b  (derivable  also  from  antenna  theory  [Gandhi 
et  al.,  1975] )  has  been  measured  for  this  region,  where 
L  is  the  major  length  of  the  body  and  2nb  is  its  weighted 
averaged  circumference.  For  an  adult  human  being, 
LI2b  —  6.34,  Srcs  —  4.21.  This  corresponds  to  SAR* 
(/-,n/1.75)  —  2.15  for  an  incident  power  density  of  10 
mW/cm2 .  For  a  man  of  height  1 .75  meters,  the  resonant 
frequency  is  on  the  order  of  62  to  68  MHz  and  at 
resonance,  the  SAR  is  1 .3  times  the  BMR. 

Region  4  -  Supraresonant  region,  i.e.,  to  frequencies 
on  the  order  of  1 .8  Sres  times  the  resonant  frequency  ft 
(for  human  beings  this  covers  the  region  fr<f<  8  /r). 
A  whole-body  absorption  rate  reducing  as  (flfr)~l  from 
the  resonant  value  has  been  observed. 

Region  5  -  /  »  /r  region.  The  S  parameter  should 
asymptotically  approach  the  “optical”  value,  which 
is  (1  -  power  reflection  coefficient)  or  about  0.5. 

In  comparing  the  graph  of  Figure  1  with  those  of 
Figure  2  (drawn  to  the  same  scale),  the  following  points 
are  noteworthy: 

1.  For  frequencies  (4-5)  /r  <  f  <  (8-9)  /r,  there  is 
little  distinction  between  averaged  rates  of  energy 
absorption  for  the  various  polarizations.  For  humans, 
this  corresponds  to  (frequency  in  MHz)  X  (height  in 
meters/1.75),  to  a  range  of  frequencies  between  250 
and  550  MHz. 

2.  The  resonances  for  k||C  and  H||C  orienta¬ 
tions  are  not  very  sharp.  Indeed,  for  the  H  1 1  C  orienta¬ 
tion,  the  SAR  gradually  reaches  a  peak  value  and  stays  at 
that  value  at  higher  frequencies. 


3.  ENERGY  DEPOSITION  IN  RATS  AND  MICE  FOR 
FREE-SPACE  IRRADIATION 


Aiuenna  theory  and  the  observed  frequency  depend¬ 
ence  that  is  shown  in  Figure  1  have  been  used  to  develop 
the  following  empirical  equations  for  rats  and  mice  and 
for  other  animals  that  may  be  represented  by  an  equiva¬ 
lent  prolate  spheroid  of  dimensions  I.  and  2b  along  the 
major  and  minor  axes,  respectively. 


Resonant  frequency  f(  =  1 1 A/L 
Sres  =  0.665  Lf2b 

where 

L  3 

L/2b  =  0.724  - - - 

l  mass  of  the  animal  in  g  j 


(1) 

(2) 

(3) 


For  the  supraresonant  frequency  region  ft  <  f  < 
(1.8  S[es)/r,  and  at  an  incident  power  density  of  10 

niW/cm! , 


SAR  in  W/kg  =, 


59.54 


GHz  Mass  in  8 


(4) 


The  equations  have  been  used  to  generate  the  general¬ 
ized  curves  for  rats  and  mice  (Figure  3)^under  free-space 
irradiation  for  the  polarization  (  E  1 1 L  )  that  produces 
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Eig.  1.  Whole-body  averaged  SARs  and  the  relative  absorption  coefficient,  S,  of  saline-filled  models 
of  human  beings  during  microwave  irradiation  in  the  free  field.  The  E'-field  vector  is  parallel  with  the 
long  axis  of  the  model.  Direction  of  propagation  is  normal  to  sagittal  plane.  Legend:  O  =  prolate 
spheroids  of  (L/2b)  =  6,  at  710  MHz;x  =  human-shaped  dolls  at  500  MHz  (parallel-plate  transmission 
line);  •  =  human-shaped  dolls  at  987  and  2450  MHz;  a  =  accurately  proportioned  figurines,  7  to  41 
cm  in  length,  at  987  MHz;  A  =  the  proportioned  figurines  at  2450  MHz. 


the  highest  SAR.  A  cross-sectiona!  enhancement  factor 
Sres  of  2.16  (given  from  equation  (2)  [also,  Gandhi  et 
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fig.  2.  Whole- body  averaged  SARs  of  saline-filled  models  of 
human  beings  during  microwave  irradiation  in  the  free  field. 
Dashed  line:  best  fit  with  data-points  where  propagation  vector, 
k.  is  parallel  with  long  axis  of  the  model  (from  head  to  toe)  and 
/  -field  vector  is  normal  to  sagittal  plane.  Continuous  line:  best 
fit  with  data  points  where  //-field  vector  is  parallel  with  long 
axis  of  model  and  /  -field  vector  is  normal  to  sagittal  plane. 
Legend  •  ~  human-shaped  dolls  at  987  and  2450  MHz; 
v  —  proportioned  figurines  at  2450  Mil/;  f»  =  human-shaped 
dolls  at  987  and  2450  MHz;  =  proportioned  figurines  at  2450 
MHz. 


al.,  1975]  for  an  aspect  ratio  L/2b  =  3.25  equivalent 
prolate  spheroid)  is  assumed  in  the  graphics  of  Figure  3. 
Also  shown  in  the  figure  are  the  values  calculated  from 
the  lethality  data  of  Schrol  and  Hawkins  [1975].  In 
calculating  the  SAR  from  their  data,  it  is  assumed  that: 
(a)  an  elevation  of  temperature  of  7.4  °C  results  in  a 
convulsion,  and  (b)  The  &T  occurs  in  an  adiabatic  man¬ 
ner;  i.e  ,  the  rate  of  energy  deposition  is  high  enough  so 
that  there  is  no  olher  mechanism  for  gain  or  loss  of 
heat  from  the  animal’s  body. 

Because  of  the  high  power  density  (150  mW/cm2)  used 
in  the  experiments  by  Schrot  and  Hawkins,  the  second 
assumption  may  be  quite  justified.  From  inspection  of 
Figure  3,  it  can  be  seen  that  there  is  a  good  correlation 
between  empirical  and  experimental  values. 


4.  WHOLE-BODY  ABSORPTION  BY  MAN  IN 
F.LFCTRICAL  CONTACT  WITH  THE  GROUND 

Whole-body  absorption  by  a  model  of  a  human  being 
with  feet  in  electrical  contact  with  the  ground  is  shown 
in  Figure  4.  The  results  are  projected  from  measurements 
with  10. 2,  12.7.  15.2.  and  20.3  cm^height)  saline-filled 
figurines  tha’  were  exposed  to  E  1 1  L  radiation  at  differ¬ 
ent  frequencies  in  the  monopole-above-ground  radiation 
chamber  \Gandhi,  1975].  This  chamber,  which  uses  a 
radiator  that  is  a  quarter-wave  above  ground,  provides  in 
conjunction  with  figurines  a  proper  simulation  of 
grounding  effects  on  energy  absorption.  Comparing  the 
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I  ig.  3.  Theoretical  projections  (continuous  lines)  of 'whole-body 
SARs  of  mice  and  rats  of  varying  mass  and  axial  length  for 
microwave  irradiation  in  the  free  field.  Data  points  are  based  on 
experiments  in  the  free  field,  or  on  irradiation  in  the  monopole- 
a bo ve- ground  chamber. 


results  to  those  of  free-space  radiation,  peak  absorption 
occurs  at  a  frequency  about  one-half  the  value  for  the 
ungrounded  condition.  The  new  resonant  frequency  is 
projected  to  be  (31-34)  X  1 .75//. nl  MHz.  At  resonance, 
the  SAR  is  about  twice  thar  of  irradiation  in  free  space. 


5.  ENERGY  DEPOSITION  IN  THE  PRESENCE  OF 
FLAT  REFLECTORS 

Judging  from  the  success  of  antenna  theory  in  predicting 
resonant  frequencies  and  absorption  cross  sections 
of  different  bodies  under  grounded  and  ungrounded 
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Fig.  4.  Projected  SARs  for  a  human  being  of  height  L  (in 
meters)  as  a  function  of  /./X.  Feet  arc  in  conductive  contact 
with  ground.  Sample  size  (AO  is  shown  in  parentheses  in  body 
of  figure. 


conditions,  highly  enhanced  SARs  are  expected  foi 
resonant  bodies  near  Hat  (180°)  reflectors  and  near  OCF 
and  other  corner  reflectors  [Jasik,  1961).  Gains  or  en¬ 
hancements  of  absorption  that  are  associated  with  half¬ 
wavelength  dipoles  in  the  presence  of  a  typically  en¬ 
countered  corner  reflector  are  given  in  Figure  5.  To 
test  the  validity  of  antenna  theory  in  predicting  enhance¬ 
ments  of  rates  of  energy  absorption,  experiments  were 
performed  with  a  near-resonant  (/./X  =  0.417)  figurine 
placed  at  different  distances  in  front  of  a  fiat  reflector 
of  dimensions  2.5  X  X  2.5  X.  The  enhancement  factors 
relative  to  the  free-space  value  and  corresponding 
SARs  are  plotted  in  Figure  6.  Also  plotted  in  the  figure 
are  the  values  of  l:1  that  would  be  expected  for  a 
standing  wave  if  the  target  were  absent.  By  comparing 
Figures  5  and  6.  one  may  note  that  the  degree  of  en¬ 
hancement  of  the  SAR  is  given  more  accurately  by 
antenna  theory  than  by  standing-wave  theory.  As  pre¬ 
dicted  by  antenna  theory  (Figure  5),  a  higher  SAR  is 
observed  for  d  =  0.125  X  Ilian  for  0.25  X  because  of 
coupling  of  the  target  to  its  image.  Also,  a  maximum 
gain  of  4.7  was  observed  experimentally.  This  value  is  in 
excellent  agreement  with  antenna  theory  and  is  1 77; 
higher  than  the  value  of  4.0  predicted  on  the  basis  of 
standing-wave  theory. 

Additional  points  to  note  in  Figure  6  are:  (a)  Non-zero 
energy  deposition  at  d  =  0.5  X.  which  may  be  ascribed 
to  absorption  due  to  magnetic  fields  at  this  setting, 
and  (b)  An  enhancement  factor  of  less  than  4  at  d  = 
0.75  X.  which  may  be  due  to  insufficiency  of  the  reflec¬ 
tor's  dimensions.  From  antenna  theory,  a  reflector  with 
a  width  larger  than  4  d  (in  the  //  plane)  is  recommended 
to  simulate  an  infinite  plane. 

The  important  features  of  the  above  results  have  been 
confirmed  by  experiments  that  were  performed  with 
100-g  Wistar  rats;  the  rats  were  placed  at  different  loca¬ 
tions  in  front  of  the  flat  reflector  and  were  subjected  to 
incident  fields  at  their  resonance  frequency  of  987 
MHz.  Tiic  details  and  the  results  of  these  experiments 
will  be  given  elsewhere  [Gandhi  and  Hunt,  1977). 

The  enhancement  factor  as  a  function  of  L/X  was 
measured  (Figure  7)  with  different  sizes  of  figurines  at 
987  and  2450  MHz.  Because  of  the  fairly  broad  lateral 
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Fig.  5.  Gain  in  rate  of  energy  received  by  a  half* wavelength 
antenna  located  at  various  distances  from  90  -comer  and  180- 
flat  reflectors. 
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Fig.  6.  The  data  points  are  measured  values  of  whole-body 
averaged  SARs  of  a  model  of  a  human  being  of  height  L  (in 
meters)  that  is  exposed  near  a  flat  reflec.or.  The  values  of  H2  in 
the  absence  of  the  target  are  shown  by  the  dashed  line. 


dimensions  of  the  larger  targets,  a  complete  set  of  data 
was  possible  only  for  d  =  3  X/4.  The  points  to  note  in 
Figure  7  are:  (a)  The  enhancement  in  SAR  at  frequen¬ 
cies  above  resonance  diminishes  and  asymptotically 
approaches  a  value  of  about  2,(b)  At  frequencies  below 
resonance,  a  higher  gain  in  SARs  from  the  corresponding 
free  space  values  (Figure  1)  is  observed,  and  (c)  For 
d  =  0.125  X,  a  reduced  enhancement  factor  is  observed 
for  smaller  targets,  which  is  reminiscent  of  the  standing- 
wave  point  of  view. 


6.  ENHANCED  ABSORPTION  CAUSED  BY  90° 
CORNERS  OF  REFLECTING  SURFACES 

Using  the  near-resonant  (L/X  =  0.417)  saline-filled 
figurine,  energy  deposition  rates  were  measured  at 


different  locations  along  the  diagonal  of  the  90°  corner 
of  several  reflectors,  each  of  dimensions  3.5  X  X  2  X(in 
the  E  direction).  The  measured  enhancement  factors 
relative  to  the  free-space  value  and  the  corresponding 
SARs  for  man  are  plotted  in  Figure  8.  Also  shown  for 
comparison  in  the  same  figure  are  the  values  of  E 2  as 
calculated  from  standing-wave  theory  (the  kind  of 
values  that  will  be  registered  by  a  so-called  power- 
density  meter  that  is  calibrated  for  free-space  measure¬ 
ments).  As  with  the  flat  reflector,  the  results  here  also 
conform  to  predictions  based  on  antenna  theory  (Figure 
5).  Enhancement  factors  as  large  as  27  are  observed  for 
d/\  =  1.5.  Another  point  of  agreement  with  antenna 
theory  is  a  plateau  in  the  enhancement  factor  at  a  value 
of  approximately  10.5  for  d/\  <  0.45.  Reduced  gains  of 
hoi  spots  that  are  closer  to  the  comer  compared  to 
d  =  1 .5  X  are  predicted  from  antenna  theory  [Jasik, 
1961 )  because  of  mutual-impedance  effects. 

A  measured  enhancement  factor  of  27  (Figure  8)  in 
the  SAR  is  considerably  larger  than  the  value  of  16  that 
is  predicted  by  standing-wave  theory  and  is  even  higher 
than  the  value  of  20  that  is  given  by  antenna  theory  for 
a  90  corner  of  a  semi-infinite  reflector.  Gains  higher 
than  those  for  semi-infinite  reflectors  have  previously 
been  observed  [Cottony  and  Wilson,  1958]  for  comer- 
reflector  antennas  of  appropriate  combinations  of 
width  and  length.  It  turns  out  that  our  reflector’s  dimen¬ 
sions  of  3.5  X  X  2  X  are,  fortuitously,  the  result  of  one 
of  the  super-gain  combinations.  There  are  many  other 
combinations  of  width  and  length  [Cottony  and  Wilson, 
1958]  that  may  also  result  in  high  SARs. 

To  confirm  predictions  of  highly  enhanced  SARs, 
experiments  were  performed  with  male  Wistar  rats 
(104  ±  8.3  g)  that  were  placed  0.44  X  and  1.5  X  from 
the  corner  of  the  reflectors.  The  d  =  0.44  X  location  was 
the  minimum  distance  to  the  corner  possible  with  the 
relatively  microwave-transparent  rat  holder  (Figure  9) 
that  we  used  in  our  experiments.  Also  for  the  d  =  0.44  X 
experiments,  90°  corner  reflectors,  each  of  dimensions 
2  X  X  2  X  (rather  than  3.5  X  X  2  X),  were  used.  From 
antenna  theory,  reflectors  of  dimensions  2  d  (=  0.88  X) 
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I  ig.  7.  I  nhunccd  SARs  of  a  body  near  a  flat  reflector  as  a  func¬ 
tion  of  the  body's  length  in  relation  to  wavelength  of  incident 
radiation. 


Fig.  8.  Whole-body  averaged  SARs  of  a  model  of  a  human  being 
near  a  90°-corner  reflector.  The  dashed  line  gives  values  of 
h.  in  the  absence  of  the  model. 
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Fig.  9.  Semischematic  drawing  of  polycarbonate  chamber  that 
was  used  in  irradiation  of  young  rats. 


X  0.6  X  would  be  quite  adequate  to  simulate  the  semi¬ 
infinite  sheets  that  would  ideally  be  needed  for  90° 
corners. 

The  mean  times-to-convulsion  together  with  the 
standard  deviations  for  four  animals  for  power  densities 
of  3  to  20  mW/cm2  are  shown  in  Figure  10.  For 
d  =  0.44  X,  the  mean  time-to-convulsion  of  225  seconds 
(for  ~  100-g  animals)  during  irradiation  at  15  mW/cm1 
corresponds  to  the  average  of  260  seconds  that  was 
observed  by  Schrot  and  Hawkins  [  1 975 1  for  free-space 
irradiation  at  150  mW/cm2.  It  is  apparent,  therefore, 
that  the  presence  of  a  90°-  corner  reflector  causes  an 
enhancement  in  the  SAR  by  a  factor  of  almost  1 1  as 
predicted  from  antenna  theory  [Jasik,  1961].  The  rela¬ 


tive  enhancements  at  the  two  locations  may  be  com¬ 
pared  by  noting  that,  at  power  densities  of  4.5  and  6 
mW/cm  for  d  =  1 .5  X,  times-to-convulsion  correspond, 
respectively,  to  9.85  and  13.02  mW/cm2  for  d  =  0.44  X. 
For  a  target-to-corner  separation  distance,  d,  of  1.5  X, 
a  deposition  rate  larger  by  a  factor  of  2.18  as  compared 
with  d  =  0.44  X  is  consequently  obtained  from  lethality 
experiments.  This  result  is  amazingly  close  to  the  rela¬ 
tive  enhancement  of  2.15  observed  between  the  two 
locations  with  man-shaped  dolls  (Figure  8). 

The  enhancement  factors  measured  for  different 
values  of  L/X  for  a  90°-  corner  reflector  are  shown  for 
d/\  =  0.44,  1.5,  and  2.5  in  Figure  11.  As  with  a  flat 
reflector,  the  enhancement  in  SAR  at  supraresonant 
frequencies  diminishes  and  asymptotically  approaches  a 
value,  in  this  case,  of  about  8.  Also,  for  closer  spacing 
to  the  corner,  a  reduction  in  the  enhancement  factor 
occurs  for  smaller  values  of  L/\.  An  asymptotic  value  of 
9  would  be  projected  from  the  simple-minded  argument 
that  for  targets  large  compared  to  a  wavelength,  a  three- 
sided  rather  than  the  four-sided  exposure  is  more 
realistic.  The  wave  that,  for  smaller  targets,  would 
bounce  off  the  corner  is  largely  masked  for  large  values 
of  /,/X. 

From  the  point  of  view  of  standing-wave  theory,  the 
fields  at  different  locations  in  the  vicinity  of  a  90° - 
corner  reflector  can  be  written  (sec  Figure  12)  as: 

E  =  E,  +  E2  +  E3  +  E4  =  E-  E 

’  (5) 

+e-vkd  _  2  e-ikd  cosky] 

where  k  =  27r/X.  Peaks  of  the  E  field  are  consequently 
anticipated  not  only  for  axial  locations  but  also  for 
several  spots  elsewhere  (some  of  which  are  shown  in 
Figure  1 2)  in  the  region  in  front  of  the  reflecting  corner. 

In  order  to  verify  the  existence  of  hot  spots  elsewhere 
in  the  region  of  a  90°-comer  reflector,  whole-body 


Fig.  10.  Convulsive  latencies  (±  SD)  of  rats  (averaged  mass  near 
100  g)  as  a  function  of  measured  power  density  of  incident  987- 
MHz  radiation  for  two  distances  of  the  animal  from  the  corner 
of  the  reflector.  (The  dashed  line  and  the  data  for  d  =  3  X/ 2 
were  obtained  with  a  corner  reflector  of  dimensions:  3.5  X  X 
2.0  X). 
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tig.  II.  Enhancement  of  SARs  in  a  body  near  a  90°-corncr 
reflector. 
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I  ip.  12.  A'- field  maxima  near  a  9(1  -corner  reflector, 
circles  denote  areas  where  a  biological  body  will  absorb  maximal 
quantities  of  energy. 


SARs  were  measured  (see  Figure  Id)  at  c//A  =  1.0  with 
gpd  without  the  26.5”  tilt  of  the  re  Hectors  around  the 
E  axis.  The  “cold"  region  without  the  tilt  becomes  a 
fairly  "hot”  region  upon  tilting.  Another  implication  of 
this  result  is  that  regions  associated  with  high  rates  of 
energy  deposition  arc  also  strongly  dependent  on  direc¬ 
tionality  of  the  incident  field. 

In  view  of  these  observations  and  also  since  the  hot 
spots  may  shift  rather  readily  upon  placement  of  other 
targets  in  close  proximity,  the  entire  region  may  be  con¬ 
sidered  as  one  that  is  potentially  capable  of  creating 
large  field  enhancements. 


( r  |  Oil*  I) 

\  ,j  d  /X  -- 1  0 


Both  for  fla^  and  for  90°Acorner  reflectors,  targets 
placed  in  k  1 1  L  and  in  H  1 1  L  orientations  gave  SARs 
that  were  considerably  enhanced  as  compared  with  their 
respective  free-space  values.  In  general,  SAR  enhance¬ 
ments  commensurate  with  target-averaged  E  values 
were  measured  for  L/ A  on  the  order  of  0.5  or  less. 
For  larger  /./A,  enhancements  lower  than  those  predicted 
from  a  simple-minded  /-.'-averaged  basis  were  measured 
and  these  are  ascribed  to  significant  perturbation  of  the 
Helds  by  the  targets. 


7.  CONDITIONS  OF  HIGHEST  RATES  OF  ENERGY 
DEPOSITION  IN  MAN  AND  ANIMALS 

As  discussed  above,  some  of  the  highest  SARs  in 
models  of  man  and  animals  are  observed  for  frequencies 
close  to  the  region  of  resonance.  Fairly  strong  enhance¬ 
ments  are  possible  at  higher  frequencies  in  the  presence 
of  reflecting  surfaces,  but  the  enhancement  of  SARs 
here  is  considerably  less  than  that  which  occurs  at 
resonant  frequencies.  Enhancement  factors  just  as  large 
or  somewhat  larger  arc  observed  for  frequencies  below 
resonance,  but  the  net  SARs  here,  too,  are  not  as  large 
because  of  the  rapid  drop-off  (see  Figure  l)in  the  SAR 
on  the  low  frequency  side  of  the  resonant  curve. 

In  order  to  appreciate  fully  the  reflector-caused  en¬ 
hancements  of  SARs,  it  should  be  mentioned  that 
reflectors  with  dimensions  2  d  X  1 .2  L  (in  the  E  direc¬ 
tion)  arc  considered  adequate  to  simulate  infinite  planes. 
Furthermore,  the  reflecting  surfaces  need  neither  be 
good  conductors  nor  solid  in  construction  to  cause 
enhancements.  Indeed,  surfaces  of  insulating  material 
with  conducting  rods  (oriented  along  E  fields)  that  are 
spaced  0.1  A  act  quite  effectively  as  solid  conducting 
surfaces  [Jasik,  1961). 

Whole-body  rates  of  energy  absorption  for  models  of 
man  and  for  animals  subjected  to  10  mW/cm2  incident 
fields  have  been  calculated  for  various  exposure  condi¬ 
tions  and  are  tabulated  in  Table  1.  The  rates  of  energy 
deposition  at  resonant  frequencies  in  the  presence  of 
a  ground  plane  or  a  reflector  markedly  increase  and, 
with  grounding  and  reflection  in  combination  they 
are  truly  staggering. 

8.  CONCLUSIONS 
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I  iL’.  13.  Whole-body  averueed  SARs  ol  models  of  human  beinys 
lor  an  E||L  orientation  in  a  90°-eorner  reflector.  Numerical 
values  in  parentheses  are  enhancement  factors  relative  to  the 
tree  field  SAR. 


Generalized  curves  have  been  given  for  whole-body 
absorption  of  electromagnetic  energy  by  models  of 
human  beings  as  a  function  of  frequency  for  free  space, 
for  electrical  contact  with  ground,  and  for  the  pre¬ 
sence  of  reflecting  surfaces.  Based  on  these  curves, 
generalized  graphs  for  specific  absorption  rates  for 
rats  and  mice  are  drawn  and  the  results  are  shown  to 
be  in  good  agreement  with  experimental  data.  Tremen¬ 
dous  enhancements  in  rates  of  energy  deposition  result 
for  targets  that  are  placed  in  proximity  to  reflecting 
surfaces.  Energy  deposition  rates  are  given  for  models  of 
man  and  for  animals  subjected  to  fields  at  a  power 
density  of  10  mW/cm2  under  various  conditions  of 
exposure;  SAR  values  as  high  as  35  to  70  times  the  BMR 
are  predicted  for  adult  human  beings  at  resonant  fre- 
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TABLL  1.  Rate  of  whole-body  energy  absorption  in  watts  by  man  and  animals  ji  10  mVi/cnr . 

70  kg  Man  100  g  Rat  400  g  Rat  I  70  kg  Man  100  g  Rat  400  g  Rat 

(1.75  m  Height)  (1 1.5  cm  Length)  (18.9  cm  Length)  (175mHeight)  (1 1  5  cm  Length)  (18  9  cm  Length) 


38  Watts 
(no  reflectance) 
0.54  W/kg 

19  Watts 

(50%  reflectance) 
0.27  W/kg 


I.  At  resonance  for  free  space. 


151  Watts 
2.16  W/kg 
f  -  62-68  MHz 


0.8  Watts 
8  W/kg 
f  -  987  MHz 


2  Watts 
5  W/kg 
f  -  600  MHz 


II.  At  resonance  for  conditions  of  electrical 
contact  with  the  ground  plane 


2  *  151  =  302  Watts  _ _ 

4.31  W/kg 

f  =  31-34  MHz  _  _ 

III.  At  resonance  for  placement  in  front 
of  a  flat  reflector. 


4.7  »  151  =  710  Watts 
10.14  W/kg 
f  -  62-68  MHz 


d  -0.125  X 

i  3.8  Watts 
38  W/kg 
f  -  987  MHz 


quencies.  The  times-to-convulsion  of  ~  100-g  rats  for 
incident  waves  of  power  densities  of  3  to  20  mW/cm2 
confirm  some  of  the  predictions  of  enhanced  absorption 
in  the  presence  of  reflecting  surfaces. 
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IV  At  resonance  for  placement  in  a  90° 
corner  reflector 


d  -  1.5  X 

27  x  151  =  4077  Watts  21  6  Watts 
58  24  W/kg  216  W/kg 

f  -  62-68  MHz  *  ~  987  MHz 


54  Watts 
135  W/kg 
f  =-  600  MHz 


V.  At  resonance  in  electrical  contact  with  ground 
plane,  in  front  of  a  flat  reflector 


d  =  0 125  X 

2  x  710  =  1420  Watts _  _ 

20.28  W/kg 

f  ~  31-34  MHz  _ _  _ 

VI.  At  resonance  in  electrical  contact  with  ground 
plane,  in  a  90°  corner  reflector. 


— ^|dh — 


9.4  Watts 
23.5  W/kg 
f  -  600  MHz 


d  =  1.5  X 


2x  4077  =  8154  Watts 
116.48  W/kg 
f  -  31-34  MHz 
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Radiators  for  microwave  biological  effects  research  -  waveguide  slot  array  with  constant 

radiation  intensity 
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A  longiluHinal-stiunl  slot  attay  with  reflectors  Ini'  been  designed,  built,  and  tested.  The  antenna 
is  easily  matched  near  the  design  frequency.  The  radiation  pattern  has  an  intensity  within  ±  0.5  dB 
over  a  total  angle  of  15.6°  in  the  11  plane  and  15.5°  in  the  E  plane.  Gain  is  19.7  dB  which  is 
3.9  dB  above  that  of  a  horn  with  a  similar  work  area.  By  using  energy-absorbing  material,  it  is  possible 
to  produce  a  pattern  in  which  the  fields  drop  suddenly  from  the  nearly  constant  values  to  nearly 
zero. 
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1 .  INTRODUCTION 

Microwave  biological  studies  with  large  targets  or 
multiple  targets  require  a  large  working  area  that  pro¬ 
vides  near-uniform  exposures  to  energy.  At  high  micro¬ 
wave  frequencies  where  high-power  sources  are  not  read¬ 
ily  available,  it  is  desirable  to  use  a  radiator  that  provides 
a  constant  intensity  beam  so  the  working  area  may 
receive  almost  all  of  the  radiated  energy. 

A  longitudinal-shunt  slot  array  with  reflectors  has  been 
designed,  built,  and  tested,  and  has  been  found  to  have 
characteristics  that  arc  suitable  for  use  in  microwave 
biological  studies.  Design  theory  and  characteristics  are 
described  in  the  following  sections. 


2.  DESIGN  THEORY 

If  finile  width,  thickness,  and  conductivity  of  the  wall 
are  neglected,  then  the  radiation  of  a  longitudinal 
resonant  slot  is  identical  to  that  of  a  half-wave  dipole 
of  same  orientation  if  E  and  H  fields  arc  interchanged 
\Compton  el  at.,  1969] .  A  collinear  array  of  half-wave 
dipoles  may  therefore  be  designed  to  obtain  a  desired 
pattern  in  the  plane  of  the  array,  and  the  design  param¬ 
eters  may  then  be  used  for  a  longitudinal-shunt  slot 
array  with  a  similar  pattern  of  radiation  intensity. 

A  half-wave  dipole  antenna  that  is  located  at  the 
origin  and  oriented  on  the  z  axis  will  produce  radiation 
with  electric  field  intensity 

E=  {jnfm  cos[(rr/2)cos0  |  /(2rrrosin0) } 

(I) 

•[exp/(wf-*ro)]0 

where  r?  is  the  intrinsic  impedance  of  free  space,  lm  is 
the  peak  current,  w  is  the  angular  frequency,  k  is  the 
magnitude  of  the  frcc-spacc  propagation  vector,  r0 
is  the  distance  from  the  point  of  origin  to  the  point 
of  observation,  and  0  is  the  angle  front  the  positive  z 
axis  to  the  point  of  observation. 


If  a  collinear  array  is  formed  with  an  odd  number  N 
of  symmetrically  fed  half-wave  dipoles  that  have  center- 
to-center  spacing  d  and  peak  current  in  the  nth  dipole 
/„,  then  superposition  using  equation  (1)  gives  the 
total  E  field 


/i?cos  [(tr/2)cos0] 
rrrosin0 


[exp/(cjf  —  tr0)J 


(2) 


t/V— 1  )/2 

•  [  £  /  cos(/ilr/cos0)]0 

n=0  " 

Equation  (2)  was  used  to  design  a  collinear  array  of  half- 
wave  dipoles  having  a  radiation  intensity  that  is  nearly 
constant  within  a  specified  beam  width  and  is  nearly 
zero  elsewhere.  Weighted  least  squares  were  used  to 
reduce  the  overshoot  and  ripple  found  in  a  Fourier 
synthesis  of  a  discontinuous  pattern. 

Let  the  desired  pattern  be  given  by 

hA  ^  ff  -  {  rr  +  6 

- [exp/(u>/  —  *ro)]0  for - <0<  - 

{r0  2  2 

(3) 

0  elsewhere 


Define  the  weighting  factor  as 

1 

f 


f0,  IzJ  <  o  <  ZJtl 


(4) 


A  elsewhere 

For  weighted  least  squares,  we  wish  to  minimize  the 
quantity: 
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f  cos[(n/2)cos0]  (A^ — 1)/2 

5=  /  { -  2  /  cos(«A:rfcos0)}‘!  M0 

'  sin  5  «=0  " 


(ir+»/2 

r  COs[(7t/2)cOS0}  (/V  — 1)/2 

+  /  {/l  —  — T— -  2  /  cos  (nkd  cos  6)  }2  dd 

J  sin0  n=o  " 

(r-W 


v 


(5) 


3 


3 


) 


n 

f  cos[(rr/2)cos0]  (N—\  )/2 

+  /  { - . — — —  —  2  I  cosfnkd  cos  0)}~  &d 8 

J  sm0  n=  0  " 

{-n+m 


Set  dS/d/m  =  0  and  use  symmetry  of  the  integrands  in 
equation  (5)  about  0  =  rr/2  to  obtain 


(ir-0/2 


/ 


cos[(tr/2)cos0] 


{ - : 


sin  8 


{N—l)/2 

Y  cos  (mkd  cos  0)  2^  ln  cos  (nkd  cos  0)&d6 


jr/2 

-/ 


A  cos[(tt/2)cos8] 
sin  8 


cos  (mkd  cos  0)  c/0 


(rr-t)/2 


cos[(7r/2)cos]  (At— 1)/2 

{ - - - J'2  cos  (mkd  cos  0)  2  /  cos  (n/crf  cos  6)dd  =  0 

sin  0  n=0  " 


Using  all  possible  values  of  the  index  m  in  equation  (6), 
there  are  (N  +  l)/2  linear  equations  involving  the 
(N  +  l)/2  unknown  currents.  The  system  of  equations 
may  be  written  using  matrix  notation 


Mmn I  = B  A 

mn  it  fn 


(7) 


The  matrix  elements  in  equation  (7)  are  given  by 
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(tr— f)/2 

-  / 
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{■ 


cos[(ir/2)cos0] 
sin  0 


}2  cos  ( mkd  cos  0)  cos  (nkd  cos  0)&dd 


n/2 

r  cosf(»r/2)cos0J  , 

+  J  { — -r^ - }  cos  (mkd  cos  0)  cos  (nkd  cos  8)d$ 

(sr-€)/2 
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if/  2 

/-  cos[(rr/2)cos0] 

B  =  /  - cos  (mfaf  cos  0)  (9) 

m  J  sin  9 

(ff-*)/2 

A  FORTRAN  program  was  written  in  which  quadra¬ 
tures  in  equations  (8)  and  (9)  were  performed  using 
Simpson’s  rule  with  due  respect  for  the  indeterminant 
form  of  the  integrand  in  equation  (8)  at  9  =  0.  The 
matrix  equation  (7)  is  then  solved  with  L  -  U  decompo¬ 
sition. 

The  beam  width  £  was  chosen  as  12.5°  to  be  small 
enough  for  the  radiation  to  approximate  a  plane  wave, 
yet  give  enough  divergence  for  convenient  chamber  di¬ 
mensions.  For  a  resonant  array  the  center-to-center 
spacing  of  slots  should  be  an  integral  multiple  of  X^/2, 
so  set  d  =  \gl 2,  or  kd  =  4.049  for  a  design  frequency  of 
10.4  GHz  with  an  RG-52/U  waveguide. 

The  ANPLOT  [Gandhi,  1974]  program  was  used  to 
find  the  radiation  pattern  for  a  series  of  collinear  dipole 
arrays  with  currents  satisfying  the  matrix  equation  (7) 
with  £  =  12.5°,  kd  =  4.049,  and  various  values  of 
N  and  A.  A  value  of  A  =  0.002  was  found  to  provide  a 
suitable  compromise  in  reducing  overshoot  and  ripple 
with  fast  rolloff  outside  the  chosen  beam  width.  The 
value  N  =  11  was  found  to  be  the  minimum  number 
of  dipoles  required  to  give  a  resonable  approximation  of 
the  desired  pattern.  Table  1  gives  the  normalized  cur¬ 
rents  found  in  the  solution  for  1 1  dipoles  with  A  = 
0.002. 

Next,  a  longitudinal-shunt  slot  array  was  designed  that 
produces  a  pattern  of  radiation  intensity  similar  to  that 
of  the  optimized  dipole  array. 

The  normalized  conductance  of  a  longitudinal  slot  at  a 
distance  2f  from  the  axis  of  the  broad  face  of  a  wave¬ 
guide  is  given  by  [Ehrlich,  1961  ] : 

g  =  2.09(Af/A)(a/b)cos2  (7rX/2X^)sin2  (nX/a)  ( 1 0) 

The  currents  in  the  fourth  pair  of  dipoles  are  inconven¬ 
iently  small,  so  slots  were  not  used  in  those  positions. 
Currents  in  the  remaining  slots  have  the  same  ratio  as 
those  in  the  dipole  array  if 

go  gi  -gi  Si-gi  =  fo7  'lt7 :122 :l37  :/s2  (11) 

For  matching  at  the  design  frequency, 

go+2g,  +2g2+2g3+2gs  =  l  (12) 


The  conductances  gn  in  Table  1  were  found  by  solving 
the  simultaneous  equations  given  by  equations  (1 1)  and 
(12).  Knowing  the  gn,  the  locations  X„  are  found  from 
equation  (10)  and  are  given  in  Table  1 .  The  slots  are  not 
collinear,  but  patterns  found  with  ANPLOT  suggest  that 
radiation  in  the  H  plane  of  the  slot  array  is  not  signifi¬ 
cantly  altered  from  that  for  a  collinear  array.  The  slot 
lengths  in  Table  1  arc  for  resonant  slots  with  width  of 
0.0625  inches  with  the  values  of  X„  [Ehrlich,  1961]. 

A  sketch  of  the  designed  longitudinal-shunt  slot  array 
is  given  in  Figure  I .  A  sliding  short  is  used  to  allow  tun- 


TABLE  ).  Current,  conductance,  location,  and  length  of  slots. 


n 

I 

n 

*„ 

X 

D 

Inches 

1 

n 

Inches 

0 

1.000000 

0.244381 

0.1781 

0.5678 

1 

0.908742 

0.201978 

0. 1598 

0.5656 

2 

0.679020 

0.112768 

0.1163 

0.5576 

3 

0.366270 

0.032811 

0.06139 

0.5508 

4 

-0.0163023 

-- 

- 

- 

5 

-0.331113 

0.030132 

0.05900 

0.3497 

ing  instead  of  using  a  fixed  sho.;  at  the  theoretical 
location  of  3X^/4  from  the  center  of  the  last  slot. 

The  radiation  pattern  from  the  slot  array  just  described 
is  acceptable  in  the  H  plane  but  the  pattern  is  slowly 
varying  in  the  E  plane.  Reflectors  are  therefore  used  to 
improve  the  pattern  in  the  E  plane. 

Reflectors  may  be  attached  to  the  sides  of  the  wave¬ 
guide  as  shown  in  Figure  2.  Let  a'  be  the  distance  be¬ 
tween  the  reflectors  in  the  plane  of  the  aperture  and  <p 
be  the  angle  from  a  normal  to  the  plane  of  the  array  to 
the  observation  point  in  the  E  plane.  If  we  consider 
only  diffraction  in  the  E  plane,  then  from  [Sherman, 
1970]  we  can  define 

U  =  (jr/a'A)sin0  (13) 

and  the  radiation  intensity  will  follow 

P  ~  (sinu/u)  2  (14) 

If  we  require  the  radiation  intensity  at  <t>  =  ±  6°  to  be 
down  only  ten  percent  from  the  peak  value,  then 
a  /X  =  1.703.  The  reflector  design  is  shown  in  Figure  2. 
A  photograph  of  the  slot  array  with  reflectors  is  given 
in  Figure  3. 

3.  ANTENNA  CHARACTERISTICS 

Figure  4  gives  the  VSWR  characteristics  of  the  slot 
array  without  reflectors.  The  sliding  short  was  adjusted 
for  minimum  reflection  at  each  frequency  of  measure¬ 
ment.  At  the  design  frequency  of  10.4  GHz,  the  mini¬ 
mum  VSWR  was  found  to  be  1.01  with  a  short  posi¬ 
tion  0.34  cm  further  from  the  center  of  the  last  slot 
than  the  theoretical  value  of  3X^/4. 

Figure  5  gives  the  VSWR  characteristics  of  the  slot 
array  with  reflectors.  The  bottom  curve  gives  the  VSWR 
with  the  short  adjusted  for  minimum  reflection  at  each 
frequency  of  measurement.  At  the  design  frequency 
of  10.4  GHz,  the  minimum  VSWR  was  1 .16  with  a  short 
position  0.50  cm  further  from  the  center  of  the  last 
slot  than  the  thoretical  value  of  3X^/4.  Tuning  with 
the  sliding  short  is  relatively  sharp  so  that  a  VSWR  of 
2.36  was  found  using  the  theoretically  shorted  position. 
Figure  5  also  gives  the  VSWR  characteristics  with  the 
short  fixed  at  positions  for  minimum  reflection  at  10.0, 
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Fig.  1.  Design  of  the  longitudinal-shunt  slot  array. 


10.2,  10.4,  and  10.S  GHz.  A  sliding  short  with  a  set 
screw  is  recommended  for  use  with  the  antenna.  All 
of  the  pattern  measurements  that  are  described  in  the 
following  were  measured  with  the  sliding  short  adjusted 
for  minimum  reflection. 

All  pattern  measurements  were  made  in  a  plane  at  a 
distance  of  2  meters  from  the  plane  of  the  array  in  an 
anechoic  chamber.  The  reflectors  were  used  for  all 
patterns. 

Figure  6  shows  both  the  theoretical  and  measured 
patterns  in  the  plane  of  the  array  or  in  the  //  plane.  The 
radiation  intensity  varied  by  less  than  ±  0.5  dB  over  an 
angle  of  15.6°.  The  theoretical  pattern  was  found  by 
neglecting  the  reflectors  and  approximating  the  slots 
with  half-wave  dipoles  as  described  in  Section  2. 

Figure  7  shows  both  the  theoretical  and  measured 
patterns  in  the  E  plane.  The  radiation  intensity  varied 
by  less  than  ±  0.5  dB  over  an  angle  of  15.5°.  The  theor¬ 
etical  pattern  was  found  by  assuming  diffraction  from  an 
aperture  having  uniform  irradiation  as  described  in 
Section  2. 

The  gain  was  determined  by  a  substitution  method  in 
which  the  same  crystal  detector  was  used  to  measure 
both  transmitted  and  received  values  of  power.  An 
open-ended  RG-52/U  waveguide  was  used  as  the  receiv¬ 
ing  antenna  at  a  distance  of  2  meters  from  the  plane 
of  the  array.  Attenuation  was  introduced  so  that  the 
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Fig.  2.  Design  of  reflectors  used  with  the  slot  array. 


two  measurements  of  power  were  made  within  1  dB. 
Gain  of  the  array  wilh  reflectors  was  found  to  be  19.7 
dB.  The  array  with  reflectors  has  radiation  intensity 
down  less  than  1  dB  from  peak  over  a  total  angle  of 
15.6°  in  the  //  plane  and  15.5°  in  the  E  plane.  The  two 
angles  define  the  work  area  if  we  will  tolerate  only  ±  0.5 
dB  variation  in  radiation  intensity.  An  ideal  radiator 
that  produces  a  constant  radiation  intensity  over  a  rec¬ 
tangle  limited  by  the  two  angles  and  zero  elsewhere 
would  have  a  gain  of  22.3  dB  so  the  gain  of  the  array 
with  reflectors  is  2.6  dB  below  ideal. 

An  optimal  horn  is  defined  as  one  where  the  dimensions 
of  the  aperture  are  adjusted  to  give  maximum  gain  for 
given  slant  lengths  of  flare  in  the  //and  E planes  [Jakes, 
1961],  Using  design  curves  given  by  Jakes,  an  optimal 
horn  will  produce  a  pattern  with  the  same  angles  for 
±  0.5  dB  variation  with  b/\  =  2.08,  aj\  =  2.87,  and  a 
gain  of  15.8  dB.  Gain  of  the  array  with  reflectors  is 
consequently  3.9  dB  above  the  horn. 


4.  INCREASING  CONTRAST  WITH  ENERGY  - 
ABSORBING  MATERIAL 

There  are  biological  applications  in  which  it  is  desired 
to  have  a  region  of  constant  irradiation  and  a  region  of 
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Fig.  3.  Slot  array  with  rcflecton. 
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near  zero  irradiation  with  a  sharp  common  boundary. 
Tests  have  been  made  in  which  a  4-inch  thick  sheet  of 
energy-absorbing  material  is  used  to  interrupt  part  of  the 
beam  from  the  slot  array. 

Figure  8  gives  the  //-plane  pattern  with  and  without 
interruption  by  a  sheet  of  energy-absorbing  material. 
Figure  9  gives  the  E-  plane  pattern  found  with  and  with¬ 
out  interruption  by  a  sheet  of  energy-absorbing  material. 
In  both  figures  diffraction  effects  broaden  the  boundary 
between  the  irradiated  and  nonirradiated  regions  to  an 
inch  or  two.  All  measurements  were  made  immediately 
behind  the  plane  of  the  dark  side  of  the  sheet  of  energy- 
absorbing  material. 


5.  CONCLUSIONS 

A  longitudinal-shunt  slot  array  with  reflectors  has 
been  designed,  built,  and  tested.  The  antenna  is  easily 
matched  with  a  sliding  short  at  the  design  frequency  and 
does  not  require  an  additional  tuner.  The  antenna 
produces  a  pattern  with  a  radiation  intensity  that  is 
constant  within  ±  0.5  dB  over  a  total  angle  of  15.6° 
in  the  H  plane  and  15.5°  in  the  E  plane.  Gain  is  3.9  dB 
above  that  of  a  typical  horn  with  a  similar  work  area. 

All  design  parameters  are  given  for  the  10.4-GHz 
antenna  but  the  design  theory  is  readily  applicable  to 
other  microwave  frequencies. 


Fig.  7.  /■/-plane  pattern  of  slot  array  with  reflectors  at  10.4 
Fiji.  5.  Reflection  characteristics  of  slot  array  with  reflectors.  CHz  at  a  distance  of  2  meters. 
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Fig.  8.  //-plane  pattern  of  slot  array  with  reflectors  at  10.4  Kig.  9.  /'-plane  pattern  of  the  slot  array  with  reflectors  at 
GHz  at  a  distance  of  2  meters  with  and  without  Eccosorb.  10.4  GHz  at  a  distance  of  2  meters  with  and  without  Eccosorb. 
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Abstract — When  pulse  functions  are  used  in  moment-method  solu¬ 
tions,  failure  to  allow  for  variation  of  the  field  within  each  cell  limits 
the  maximum  usable  electrical  size  of  the  cells.  Appreciable  error  is 
expected  for  |  fc  |  /  >  2  in  one  or  two  dimensions  and  |fc|/  >  v/6inthe 
three-dimensional  problem  where  /  is  the  side  of  a  cell  and  k  is  the 
propagation  constant  in  the  material. 

I.  Introduction 

N  ELECTROMAGNETICS,  discretization  for  transfor¬ 
mation  of  an  integral  equation  to  a  matrix  equation  is 
often  accomplished  using  pulse  functions  as  a  basis  [1]  [4]. 
When  using  pulse  functions,  the  scatterer  is  partitioned  into 
a  number  of  cells  N,  where  N  is  large  enough  that  complex 
permittivity  and  the  complex  time-independent  electric  field 
may  be  assumed  constant  within  each  subvolume. 

In  a  given  problem,  it  is  common  to  estimate  the  maxi¬ 
mum  usable  cell  size  and  try  several  values  of  N  to  test  for 
apparent  convergence  of  the  resulting  solutions  [5].  An 
approximate  upper  bound  on  cell  size  may  be  found  from 
observation  of  the  oscillatory  nature  of  the  kernel  in  the 
integral  equation,  which  suggests  that  the  cell  size  not  exceed 
A0  /5,  where  20  is  the  free  space  wavelength  [6],  It  is  the  object 
of  this  paper  to  establish  a  significantly  tighter  upper  bound 
on  cell  size  in  moment-method  solutions  with  scatterers 
when  pulse  functions  are  used  as  the  basis. 

II.  Evaluation  of  the  Bound  on  Cell  Size 

Consider  a  source-free  region  of  space  in  which  the 
dielectric  properties  are  homogeneous,  linear,  and  isotropic. 
All  fields  are  assumed  to  have  exp  (jtot)  time  variation.  We 
may  set  up  a  local  Cartesian  coordinate  system  at  any  point 
in  the  region  and  require  that  the  homogeneous  wave 
equation  be  satisfied: 

V2£  +  k2E  =  0.  (1) 

Six  points  are  chosen  a  distance  S  from  the  origin  on  halves 
of  each  of  the  three  Jocal  axes./j  will  represent  the  value  of 
one  component  of  £,  say  Ex ,  at  the  ith  point.  Let 

6f  =fi  -L  (2) 

where /„  is  the  value  of  Ex  at  the  origin. 

If  S  is  small  enough  so  that  there  is  little  variation  in  Ex . 
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the  difference  equation  approximation  of  (1)  may  be  used 
with  the  result 

i  <5 V,  +  k2S2f0  =  0.  (3) 

/=  1 

The  greatest  (in  absolute  value)  of  the  5/j  must  satisfy 

M>i^!  f4) 

\f\  “  6  •  (  ) 

When  (4)  predicts  large  variations,  we  cannot  expect  the 
prediction  to  be  quantitative,  but  we  may  safely  infer  that 
substantial  variation  of  the  fields  will  always  occur  within 
the  volume  of  a  cell  containing  the  six  points  if  the  predicted 
fractional  variation  has  a  value  of,  say,  one-half,  for  which 

|*|S  =  aA  (5) 

Restricting  our  attention  to  cubical  cells,  we  note  that  the 
smallest  cube  containing  the  six  points  has  side 

l  =  s/2S  (6) 

so  large  variations  in  the  fields  are  expected  for 

1*1/*  a/6.  U) 

The  corresponding  result  for  square  cells  with  side  /  in  the 
two-dimensional  problem  or  linear  cells  of  length  /  in  the 
one-dimensional  problem  is 

|k|l>2.  (8) 

In  an  efficient  moment-method  solution  such  as  in  [1]— [4], 
variation  of  the  Green’s  function  within  each  cell  is  closely 
approximated  so  that  the  primary  source  of  error  is  imper¬ 
fect  representation  of  the  fields  by  the  basis.  If  pulse  func¬ 
tions  are  used,  it  is  assumed  that  the  fields  are  constant 
within  each  cell.  Then  appreciable  error  is  expected  if  (8)  is 
satisfied  in  a  one-  or  two-dimensional  problem  or  if  (7)  is 
satisfied  in  a  three-dimensional  problem.  The  smoothing 
property  of  the  integral  operator  causes  the  error  in  the 
solution  to  be  somewhat  less  than  may  be  anticipated  fora 
simple  basis,  but  we  may  still  expect  that  (7)  and  (8)  give  a 
reasonable  upper  bound  for  cell  size. 

III.  Examples 

The  bound  of  (7)  is  tighter  than  that  found  from  consider¬ 
ation  of  the  oscillatory  nature  of  the  kernel  for  scatterers 
having  a  relative  permittivity  cr  such  that  |er|  >  3.8.  In 
order  to  illustrate  this,  a  couple  of  biological  applications 
using  dielectric  properties  of  muscle  [7]  will  be  considered 
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TABLE  1 

Values  of  the  Indu  ed  Electric  Field  at  the 
Center  oe  Muscle  Cubes  at  2450  MHz 
as  Given  in  [4]:  £,/e0  =  47.0,  a  =  2.21  mho/m 


Nuabec  ! 
of  Cells 

M* 

l,Xo 

|E|  Center 
(V/o) 

Cube  one 

1 

1 

6.37 

0.144 

0.0739 

a  m  X 

27 

1/3 

2.12 

0.0481 

0.0922 

Cube  two 

1 

1/4 

1.59 

0.0361 

0.0592 

a  -  X/4 

27 

1/12 

0.531 

0.0L20 

0.0556 

TABLE  II 

Errors  in  the  Average  Specific  Absorbed  Power 
Density  for  an  Infinite  Cylinder  of  Muscle 
with  a  20-  by  20-cm  Cross  Section  at  100  MHz: 
£|  /to  =  68.0,  a  =  0.890  mhom 


Number 
of  Cells 

i 

»/xt 

i 

1/Xo 

Error 

4 

0.367 

2.76 

0.0334 

-34 . 90? 

9 

0.24  5 

1.84 

0.0222 

-16.402 

16 

0.183 

1.38 

0.0167 

-8.642 

25 

0.147 

1.11 

0.0133 

-4.962 

since  the  complex  permittivity  has  a  relatively  large 
magnitude. 

Table  I  uses  the  results  given  in  [4]  for  two  muscle  cubes 
having  sides  of  a  =  7,eand  7,t/4,  respectively,  at  2450  MHz, 
where  ).c  is  the  wavelength  in  the  material.  No  analytical 
solution  is  available  for  comparison,  but  variation  of  the 
calculated  electric  field  with  cell  size  suggests  that  there  is 
significant  error  in  the  single-cell  solution  for  the  larger 
muscle  cube.  The  small  values  of  //A0  suggest  that  the 
oscillatory  nature  of  the  kernel  should  contribute  little  error 
in  the  four  calculations.  Values  of  \k\l  and  (7)  suggest  that 
significant  error  should  be  found  using  one  cell  with  the 
larger  cube,  as  is  observed. 

For  the  second  example,  Richmond’s  method  [1]  was  used 


to  calculate  the  average  specific  absorbed  power  density  in 
an  infinite  cylinder  of  muscle  with  a  20-  x  20-cm  square 
cross  section.  A  plane  wave  at  100  MHz  with  a  power  density 
of  1  mW/cm2  incident  normal  to  one  of  the  four  congruent 
flat  surfaces  was  used  for  TM  excitation.  The  average 
specific  absorbed  power  density  found  using  100  cells  is 
0.05534  mW/cm  \  which  was  used  as  a  standard,  differing  by 
0.39  percent  from  the  value  found  using  81  cells.  Table  II 
gives  the  results  found  using  fewer  cells.  Note  that  significant 
error  is  encountered  when  (8)  is  satisfied  even  though  lj/.0  is 
so  small  that  the  oscillatory  nature  of  the  kernel  contributes 
negligible  error. 

IV.  Conclusions 

When  pulse  functions  are  used  in  moment-method  solu¬ 
tions,  failure  to  allow  for  variation  of  the  field  within  each 
cell  limits  the  maximum  usable  electrical  size  of  the  cells. 
Appreciable  error  is  expected  for  |  A:  j  /  >  2  in  one  or  two 
dimensions,  and  |  k  \  I  >  v/6  in  three-dimensional  problems, 
where  /  is  the  side  of  a  cell.  The  new  upper  bounds  for  cell  size 
are  significantly  tighter  than  those  found  from  the  oscilla¬ 
tory  nature  of  the  kernel  if  the  scatterer  has  a  large  relative 
permittivity.  The  new  upper  bounds  have  been 
demonstrated  with  two-  and  three-dimensional  solutions. 
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ABSTRACT 


Two  new  methods,  termed  "Plane  Wave  Correction"  and  "Cylindri¬ 
cal-Cell  Correction",  are  presented  for  improving  convergence  of 
moment-method  solutions  in  electromagnetics.  Convergence  is  improved 
because  the.  calculations  include  approximation  for  the  variation  of 
the  field  within  each  cell. 
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I.  INTRODUCTION 

In  electromagnetics,  discretization  for  transforming  an  integral 
equatio.  to  a  matrix  equation  is  commonly  accomplished  using  pulse 
functions  as  a  basis. 1-1>  More  elaborate  buses  such  as  modal  fields 
may  be  justified  in  a  particular  problem  but  are  avoided  in  the  matrix 
formulation  for  the  general  scattering  body. 

Discretization  with  pulse  functions  requires  that  the  scatter¬ 
ing  body  be  partitioned  into  a  number  N  of  cells,  where  N  is  large 
enough  that  complex  permittivity  and  the  complex  time-independent 

electric  field  may  be  assumed  constant  within  each  subvolume.  In  a 

2 

numerical  solution  the  required  storage  is  proportional  to  N  .  With 

Gaussian  elimination,  solution  of  the  matrix  equation  requires  a 

3 

number  of  operations  that  are  proportional  to  N  .  If  N  is  large  so 
that  matrix  storage  and  solution  of  the  matrix  equation  are  the 
dominant  factors,  the  cost  of  computation  is  proportional  to  the  product 
of  storage  and  time,  or  N"*.  Procedures  for  improving  convergence  with 
the  general  scattering  body  are  of  interest  since  even  a  13  percent 
reduction  in  the  required  number  of  cells  would  result  in  a  cost  sav¬ 
ings  of  50  percent  when  large  K  is  required. 

This  paper  describes  two  new  procedures  which  have  been  found 
to  improve  convergence  in  solution  of  the  two-dimensional  problem  of 
TM  excitation  of  an  infinite  cylinder  of  arbitrary  cross  section  shape. 
Convergence  is  improved  because  the  calculations  include  part  of  the 
variation  of  the  field  within  each  cell. 
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II.  REVIEW  OF  EARLIER  SOLUTIONS  USING  PULSE  FUNCTIONS 


For  TM  excitation  of  an  infinite  cylinder,  both  incident  and 

—V 

scattered  fields  have  E  directed  parallel  to  the  infinite  dimension 
(z  axis).  The  electric  field  is  described  by  a  scalar  integral  equa¬ 
tion: 


Ez(x,y)  =  Ez(x,y)  - 


jk 

T 


~  ||  ^r(>:’,y')  -  ljE^Cx'.y'jK^^k^^dx'dy' 


(1) 


where  (x,y)  and  (x’,y')  are  coordinates  of  the  observation  point  and 

source  point,  respectively;  p  is  the  distance  between  the  two  points; 

E*  and  E  are  incident  and  total  electric  field  intensity,  respectively, 
z  z 

with  e^<0v_  time  dependence;  e^x'.y')  is  complex  permittivity  at  the 

. 

source  point;  k  =  «oJp  e  ,  and  H'  fk  p)  is  the  Hankel  function  of  zero 
o  ;  o  o  o  \  c  ) 

order. 

The  discrete  analog  of  Eq.  1  consists  of  the  N-by-N  system  of 
linear  equations: 


N 

l 

n=l 


A  E  **  E 
ran  zn  zra 


m  =  1,  2,  — ,  N 


(2) 


Richmond  performed  the  discretization  using  pulse  functions.1 
Each  square  cell  was  approximated  by  a  circle  of  radius  a  having  an 
equal  area  to  simplify  integration  of  Bessel  functions.  His  expres¬ 
sions  for  the  matrix  elements  follow: 


A 

mm 


1  + 


jtrk  a 

1  +  - _ °— 

2 


(3) 
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(4) 


In  this  method,  the  electric  field  intensity  is  considered  to 
be  constant  over  the  area  of  a  cell,  an  approximation  that  results  in 
significant  error  if  the  cells  are  not  very  small.  Larger  cells  can 
be  used  if  the  variation  of  the  electric  field  is  included  in  the 
calculations.  A  method  of  approximating  the  variation  of  the  electric 
field  called  ’’plane  wave  correction"  is  described  in  the  next  section. 
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III.  PLANE  NAVE  CORRECTION 


Variation  of  is  initially  unknown,  but  we  may  make  the 
approximation  of  representing  the  fields  within  each  cell  by  a  super¬ 
position  of  plane  waves.  All  member  plane  waves  are  required  to  have 
propagation  vectors  of  magnitude,  k  corresponding  to  the  frequency  of 
the  incident  wave  and  local  complex  permittivity  of  the  dielectric  body. 
If  cu  and  represent  orientation  of  propagation  and  amplitude  of  the 
ith  component,  respectively,  the  superposition  of  plane  waves  may  be 
written  as: 


-jkr’  cos  ( 6 '  -  a.] 

Ez(r’,G’)  =  l  Bi  e  V  (5) 

i 

For  small  cell  size  we  may  expand  the  exponential  in  Eq.  5  to  obtain 
the  approximation 

E2(r,,0')  =  l  \l  -  jkr’  cos  (8*  -  a±)  -  k2r’2  cos2  (G'  -  ajj 

i  L 

(6) 

If  complex  permittivity  is  constant  within  a  cell  and  the  approxima¬ 
tion  of  circular  cells  is  used  in  Eq.  1,  then  the  diagonal  matrix 
elements  are: 

jk2  fe  1)  3  271 

■ 1  +  n  W  1  r’dr’  J  E2(r''e,)  “o2’  (v) d9’  <?> 

z  0  0 

From  Eqs.  6  and  7,  after  performing  the  0*  integrations. 
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I 
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I 


A  =  1 
mm 


jk2  (e  -  l\  ar 

*-sfft-Ll  \‘A» 

4  3  i  X 


*  ,2  ,2 

—  k  r 


Hq2)  (kor>)  r’dr' 


(8) 


The  summation  in  the  integrand  in  Eq.  8  has  a  factor  that  is 
independent  of  the  summation  index  and  may  ?)e  removed  from  the  summa¬ 
tion.  Consequently,  the  diagonal  matrix  elements  have  no  dependence 
on  the  unknown  coefficients  in  the  expansion,  a  fortunate  occurrence 
indeed.  Using  a  dummy  variable  in  Eq.  8, 


k  a 
o 


k  a 
o 


A  =1  +  A?  ( £  -  1\  [  H^2^(x)xdx  -  Ap  (c  -  l'W  [  H^2^(x)x^dx 

mm  2V.  rm  )  J  o  7  8  ^  n  /nnjo 


(9) 


Performing  the  integrations  over  the  Hankel  functions  gives 
for  the  diagonal  elements: 


mm 


=  1  + 


frm-1)  (erm  +  i) 


jrk  a 
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H 


rm 
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(ko* 


A 

/  i 


(koa) 
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il 

4 


(10) 


In  evaluating  the  off-diagonal  matrix  elements,  suppose  that 
the  source  point  P’  is  in  the  nth  cell  and  the  observation  point  P  is 
at  the  center  of  the  mth  cell,  as  shown  in  Fig.  1.  Let  s  and  r*  be 
distances  from  the  center  of  the  nth  cell  to  the  observation  point  and 
source  point,  respectively.  If  the  angle  between  s  and  r’  is  8’, 
then  the  distance  between  the  source  point  and  observation  point  is 
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Fig.  1.  Notation  used  in  evaluation  of  off-diagonal  matrix  elements 

=  Jr'2  +  s2  -  2r’s  cos  8’ 

1  mn  mn 


mn 


(ID 


If  is  constant  within  the  nth  cell  and  the  approximation  of  circular 


cells  is  used  in  Eq.  1,  then  the  off-diagonal  matrix  elements  are: 
2  .  a  2ir 


jk 


mn  4 


-eW)  j  ’t'dr'  \  EZ(r’>6,)  “o2>  (koP»n)  <’9'  (12) 


From  Graf's  addition  theorem,5 


H(2)  fk  p  >  “  l  H^;  fk  s  )  J.  fk  r’V  cos  *9'  (13) 

o  V  o  mnj  k  0  ran)  v  °  / 

From  Eqs.  11,  12,  and  13, 
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mn 
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Jko(ern  “  X)  H£’  (koSnm)  r  , 
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A  F.z  (0,0’) 
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6')  cos  HO  d0‘ 


(14) 
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From  Eqs.  6  and  14, 


jk2  (c  -  i\  J"  H0(2)  (ks  \ 
J  o  1,  rn  )  ot-co  ^  '  °  mn' 


Inn 


x  a*-—** 
4 1 s, 


a 

{ (V)  r'dt' 


/-  ii 

|  l  Bi  fl  -  jkr 1  cos  ^6 1  -  a.)  -  |  k2r’^  cos2(0’  -  o±)j 


cos  W’d0' 


Performing  the  O'  integrations  in  Eq.  15, 
jk2  (c  -  r\  fk  s  ) 

A  =  -°-A  rn  _JLjL  -i-o-nffll  V  B.  2r  J  (k  r')  r'dr’ 
mn  4  V  B.  j  J  i  o  \  o  ) 

J  1  0  L  J 

jk2  (e  -  l')  H(2)  (ks  j 

_  .1™ _ 1-2 - LP-JSUl  h  i  k2r’2|  J  ( k  r’)  r'dr’ 

4  J  B.  J  ^  i  2  J  o  V.  o  / 

4*  1  /\ 


2k2  (e  -  ll  Hp)  (l>  )r  , 

+  ■  °  L»U  Jjc-i-SJgal  h  B4.kr-  cos  «d  ^  (ke)  r'dr' 

V  1  0 

1  (26) 


The  summations  in  all  three  integrands  of  Eq.  16  contain  factors 
independent  of  the  summation  index  that  may  he  removed  from  the  summa¬ 
tions.  Also  changing  to  a  dummy  variable. 
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Define 


y  B.  coo  a. 
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(18) 


Performing  the  integrations  in  Eq.  17  and  using  Eq.  18, 
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(19) 


Since  only  the  last  member  in  Eq.  19  is  dependent  on  the  unknown 
coefficients  in  the  plane  wave  expansion,  we  wish  to  delete  the  term 
and  use  the  first  three  members  for  the  off-diagonal  matrix  elements: 
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The  use  of  Eq.  20  is  justified  only  if  the  last  member  of  Eq.  19  is 

small.  Fractional  error  in  A  due  to  neglecting  the  last  member  of 

mn  ** 


Eq.  .19  is: 

# 

*  Hl(2) 

(k  s  ]  k  a  J  (k  a) 

V  o  ran/  o  2  \  c  J 

mn  H<2V  s  ) 

o  Vo  mny 

[(*-  re 

J  (k  a)  +  k  a  Jofk 

o  /  J  oV  o  J  2  o  2V  o  J 

(21) 


Since  the  magnitude  of  the  ratio  of  the  Hankel  functions  is 
greatest  for  small  argument,  we  set  s^  =  4  a  for  adjacent  cells. 
Expanding  the  Bessel  functions  for  small  arguments  gives: 


|  v(c — |  k  a 

|v  I  <  -1--- — ?. - 

Bn  '^"OvrivO 


(22) 


where  T  =  e^1  »  1.78107,  with  y  =  Euler's  constant. 

In  the  plane  wave  correction  method,  Eqs.  10  and  20  are  used 
for  the  matrix  elements  and  the  solution  corresponds  to  the  value  of 
E^  at  each  cell  center.  Numerical  considerations  suggest  that  deletion 
of  the  last  member  in  Eq.  19  is  the  prir  ipal  source  of  error,  so  that 
it  is  possible  to  use  Eq.  22  with  the  condition  number  of  the  matrix 
to  give  an  error  bound  for  the  solution.  Since  no  error  bound  is 
available  in  the  method  using  pulse  functions  that  is  described  in 
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Section  II,  tlie  plane  wave  correction  has  the  advantage  of  an  expres¬ 
sion  for  the  error  bound,  in  addition  to  improving  convergence.  A 
second  method  of  approximating  the  variation  of  the  electric  field 
called  "cyJ.indrical-cell  correction"  is  described  in  the  next  section. 


IV.  CYLINDRICAL-CELL  CORRECTION 


In  the  methods  described  in  Sections  II  and  III,  the  scattering 
body  is  divided  into  cells  that  are  approximated  by  circular  cylinders 
in  calculation  of  the  matrix  elements.  The  method  of  cylindrical-cell 
correction,  which  is  presented  for  the  first  time  in  this  paper, 
emphasizes  the  properties  of  such  a  model  of  the  original  body. 

Let  the  cross  section  of  the  scatterer  be  divided  into  square 
cells,  each  having  side  s.  For  the  model,  replace  each  cell  with  a 
circular  cylinder  having  radius  a,  and  assume  properties  of  free  space 
between  the  cylinders.  For  the  chosen  geometry6  the.  effective  relative 
complex  permittivity  Eref ^  of  the  model  and  the  complex  permittivity 
of  each  cylinder  are  related  by 

ereff  =  1  +  “V  V£r  "  X)  (23) 

s  x 

If  we  require  that  the  effective  relative  complex  permittivity  of  the 
model  equal  that  of  the  original  scatterer,  then 


If  the  electrical  size  of  the  cylinders  is  sufficiently  small, 
the  inhomogeneities  of  the  model  will  have  no  significant  effect,  so 
that  scattering  from  the  model  satisfying  Eq.  24  will  duplicate  that 
of  the  physical  scatterer.  The  model  is  useful  since  we  may  approxi¬ 
mate  the  variation  of  the  electric  field  within  each  cylinder.  A 


I 


circular  cylinder  with  TM  excitation  wil I  have  internal  fields  given 
by 
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00 
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(kr)  cos  ^nO  +  c^j 


(25) 


where  the  b  and  c  are  determined  by  the  excitation.  If  the  cylinder 
n  n 

is  sufficiently  small,  the  zero  order  term  will  dominate,  so  we  may 
use  the  approximation 


E  =  b  J  (kr) 
z  o  o 


(26) 


From  Eqs.  7  and  26  the  diagonal  matrix  elements  are  given  by: 
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where  k'  n  if?"  k  . 

\  r  o 

Performing  the  9'  integration  and  using  a  relationship  from 
Ref.  7  for  the  r'  integration,  we  obtain  the  following  expression  used 
to  calculate  the  diagonal  matrix  elements: 


A  =  fk’a  J'  (k'a)  H^2)(k  a)  -  k  a  J  (k’a)  lh(2)  (k  a)  (28) 
mm  2  [  1  o  \  o  /  oo  1V.°/J 


From  Eqs.  12  and  25  the  off-diagonal  matrix  elements  are  given  by: 


Jk,  (e;„  - i) 


rar. 


2it 

r'dr’  J  (k’r’)  H(2)  (k  P  )  d0’ 
J  Jo  o  o  mn,' 


(29) 
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Using  Graf's  addition  theorem5  and  notation  of  Fig.  1: 


A  =  "7 
ran  4 


z  a  z?f 

CO  •  • 

-(c1  ~  1)  ]  H„(2)fk  s  \  J  fk  r*)  J  (k’r')r'dr’  f  cos  5,0  ’  d0 

V  rn  /  M  o  rnj  J  U  o  J  o  J 


Performing  the  0'  integration,  which  collapses  the  summation,  we  get 


cl 

A  =  J-  k2  fe'  -  l)  H(2)  fk  s  \  J  (k  *  r ' )  J  fk  r'N  r’dr' 
tin  2  o  V  rn  Jo  \,  o  mj  j  o  o  J 


Using  a  relationship  from  Ref.  6  for  the  r’  integration,  we  obtain  the 
following  expression  used  to  calculate  the  off-diagonal  matrix  elements 

A  =  K(2)  fk  s  \fk’a  J  (k  *  a)  J  fk  a)  -  k  a  J  (k'a)  J.fk  a)}  (32) 

mn  2  o  o  mn/L  1  o  /  o  o  lv.  o  /J 


m  ^  n 


In  the  cylindrical-cell  correction  method,  Eos.  28  and  32  are 

used  for  the  matrix  elements  and  the  solution  corresponds  to  the  value 

of  E  at  each  cell  center.  A  solution  is  defined  for  any  ratio  of  a/s, 
z 

but  only  two  cases  have  been  studied: 


S 

Tangent  cylindrical  cells:  a  =■  tj-  -  1  +  —  -  1^ 


Overlapping  cylindrical  cells:  a  =  ~4r  c'  «*  c 

nr  r  r 


For  very  large  or  very  small  values  of  a/s,  the  structure  of  the  model 
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differs  significantly  from  that  of  the  scattered  so  that  the  solution, 
which  actually  corresponds  to  scattering  from  the  composite  of  circular 
cylinders,  will  differ  from  that  for  the  desired  scatterer.  Advantages 
of  the  method  will  be  seen  in  the  example  in  Section  VI. 
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V.  CALCULATION  OF  POWER  ABSORPTION 


With  both  plane  wave  correction  and  cylindrical-cell  correction, 

the  solution  corresponds  to  the  value  of  Ez  at  each  cell  center.  The 

plane  wave  correction  leaves  the  variation  of  E  within  each  cell 
r  z 

undefined,  but  with  cylindrical-cell  correction  Eq.  26  describes  the 

approximate  variation  of  E^  within  each  circular  cylinder  of  the  model. 

If  the  value  of  power  absorption  is  desired,  then  1/2  cE^ E*  may 

be  computed  within  each  cell  using  the  values  of  Ez  found  using  pulse 

functions  or  plane  wave  corrections.  In  this  section  the  approximate 

variation  of  E  within  each  cell  is  used  for  an  improved  expression 
z 

for  power  absorption  if  cylindrical-cell  correction  is  used. 

Total  power  absorption  in  one  circular  cylinder  is  given  by: 


P  = 


a  2n 

H 

0  0 


1  a’E  E*  rdrdO 

2  z  z 


(33) 


From  Eq-  24  the  conductivity  of  the  cylinder  is  related  to  that  of  the 
scatterer  by 


a* 


2 

_s _ 

2 

Tra 


a 


(34) 


Substituting  Eqs.  26  and  34  in  Eq.  33  gives: 


Using  a  relationship  from  Ref.  6  for  the  r'  integration: 


P  = 


s2ob  b*  [k'a  J  (k'*a)  J.(k'a)  -  k'--a  J  (k'a)  J. (k'*a)l 
oof  O  .1  O  J.  J 

- - - y - r - ~i -  (36) 

a2  [(k'T  -  (k* *)  J 


Then  the  average  absorbed  power  density  within  the  cell  is  given  by: 


ob  b*  [k'a  J  (k’*a)  J,(k'a)  -  k'*a  J  (k'a)  J,(k'*a)l 
_  O  o  l _ o _ 1 _ _ _ O  _ ]_ _ j 

[(k'a)2  -  (k'*a)2 


(37) 


Simplifying  Eq.  37  we  obtain  the  result 


P ,  =  ^bb* 
d  2  o  o 


Ira  |  k'a 

_ i _ 


J*(k'a)  J^k'a)] 


Re(k'a)  Im(k'a) 


(38) 


where  b^  is  the  calculated  value  cf  found  using  cylindrical-cell 
correction. 
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VI.  NUMERICAL  EXAMPLE :  THE  INFINITE  CYLINDER  MODEL  OF  HAN 

The  infinite  cylinder  has  been  suggested  for  use  as  a  model  for 
body  extremities  or  the  cheat  cavity  in  evaluation  of  biological 
hazards  from  electromagnetic  radiation. R  Approximation  of  an  extremity 
or  the  v.'hole  body  by  an  infinite  cylinder  may  be  made  using  the  ratio 
of  volume  to  length  for  the  cross-sectional  area. 

The  model  considered  in  this  example  is  an  infinitely  long 
cylinder  with  square  cross  section  of  20  by  20  cm,  having  TM  excita¬ 
tion  by  a  plane  wave  incident  perpendicular  to  one  of  the  four 
congruent  faces.  The  body  is  homogeneous  with  values  of  permittivity 
and  conductivity  that  are  typical  for  muscle,  skin,  and  other  tissues 
having  high  water  content.9 

The  methods  described  in  this  paper  have  been  used  to  calculate 
average  specific  absorbed  power  for  the  model.  Suitability  of  the 
model  can  only  be  justified  at  high  frequencies  where  end  effectu  may 
be  neglected,  but  calculations  have  been  made  over  a  wide  range  of 
frequencies  to  allow  comparison  of  convergence  with  the  different 
methods. 

Figure  2  illustrates  the  frequency  dependence  of  average 
specific  absorption  rate  (SAR)  for  the  model.  The  va.lues  obtained 
using  the  accepted  method  of  pulse  functions  with  300  cells  are  the 
standard  for  comparison.  Values  obtained  using  pulse  function?  with 
81  cells  were  found  to  differ  from  those  with  100  cells  by  less  than 
1  percent  for  frequencies  up  to  200  MHz.  Values  found  using  9  cells 
with  pulse  functions,  plane  wave  correction,  and  cylindrical-cell 
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S  100  cells  pulse  functions 

©  9  cells  pulse  functions 

~  9  cells  plane  wave  correction 

i 

+  9  cells  tangent  cylindrical  cells 


X  9  cells  overlapping  cylindrical  cells 


FREQUENCY  (m. 


Variation  of  average  specific  absorption  rate  with  fre 
quency  for  the  infinite  cylinder  model  of  man. 


* _ 1  __ 


correction  with  tangent  and  overlapping  cells  are  also  given  in  the 
figure.  For  9  cells  all  four  calculations  are  in  good  agreement  with 
the  standard  at  low  frequencies,  where  the  cells  have  small  electrical 
size.  Calculations  with  the  new  methods  appear  to  have  comparable 
accuracy  and  appear  to  offer  a  significant  improvement  in  convergence 
when  compared  to  pulse  functions. 

From  Eq.  22,  |fi  }  <  3.3,  7.0,  or  18.1  percent  with  9  cells  at 
50,  100,  or  200  MHz,  respectively.  Then  significant  error  in  the 
calculations  with  plane  wave  corrections  is  not  predicted  until  the 
frequency  is  about  200  MHz ,  at  which  the  curve  deviates  from  the 
standard  in  Fig.  2.  It  appears  that  Eq.  22  is  usable  for  determining 
maximum  usable  cell  size  for  calculations  at  a  given  frequency. 
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VII.  CONCLUSIONS 


w 

t 


Two  methods  have  been  found  which  improve  convergence  in  solu¬ 
tion  of  the  two-dimensional  problem  of  TM  excitation  of  an  infinite 
cylinder  of  arbitrary  cross  section  shape.  The  extension  to  the 
three-dimensional  problem  has  not  yet  been  completed. 

Convergence  is  improved  because  the  calculations  include  ap¬ 
proximation  for  the  variation  of  the  field  within  each  cell.  The  new 
methods  give  expressions  for  the  matrix  elements  that  are  different 
from  those  found  using  pulse  functions,  but  values  of  the  matrix 
elements  are  in  good  agreement  for  electrically  small  cells. 

The  matrix  elements  found  with  the  new  methods  do  rot  require 
a  substantial  increase  in  computational  effort.  For  a  homogeneous 
scatterer,  if  equal  size  cells  are  used,  the  new  methods  require  one 
time  computation  of  several  additional  Bessel  functions,  but  time  for 
calculation  of  the  matrix  elements  ir  dominated  by  recalculation  of 
the  zero-order  Hankel  function  which  is  nee  1  with  or  without  the 
corrections  made  in  the  new  methods.  If  the  scatterer  is  not 
homogeneous  or  if  different  cell  sizes  are  used,  then  time  for  calcula¬ 
tion  of  each  matrix  element  is  approximately  doubled  in  the  new  methods. 

Time  spent  in  calculation  of  the  matrix  elements  is  proportional  to 
2 

N  ,  whereas  time  spent  in  solving  the  matrix  equation  is  proportional 

3 

to  N  ,  so  for  large  numbers  of  cells  there  is  still  no  significant 
increase  in  computation  time  with  the  new  methods. 
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APPENDIX  G 


Behavioral  and  thermal  effects  of  microwave  radiation  at  resonant  and  nonresonant  wavelengths 
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Behavioral  and  thermal  effects  of  radiating  an  animal  with  differing  wavelengths  of  microwave 
energy  at  the  same  power  density  were  investigated  in  the  first  of  two  studies.  Five  Long-Evans  rats 
were  trained  to  perform  a  lever-pressing  task  and  were  rewarded  with  food  on  a  variable  interval 
schedule  of  reinforcement.  Rats  were  individually  exposed  in  random  order  to  400-,  500-,  600-, 
and  700-MHz  CW  radiation  at  a  power  density  of  20  mW/cm1 2  with  the  long  axis  of  the  rat’s  body 
parallel  to  the  vector  of  the  electric  field.  Radiation  at  all  wavelengths  produced  rises  of  body  tem¬ 
perature  and  stoppage  of  lever  pressing.  The  averaged  rise  in  body  temperature  was  greatest  and 
work  stoppage  was  most  rapid  during  exposures  at  600  MHz.  In  the  second  study,  six  rats  were 
exposed  in  random  order  to  600-MHz  CW  radiation  at  power  densities  of  S,  7.5,  10,  and  20  mW/cm2 
while  performing  the  same  behavioral  task.  Exposures  at  10  and  20  mW/cm2  resulted  in  work 
stoppage,  while  exposures  at  5  and  7.5  mW/cm2  did  not.  In  addition,  three  of  the  rats  were  sub¬ 
sequently  exposed  while  responding  to  600-MHz  pulsed  radiation  (1000  pps,  3-  or  30-/rs  pulse  dura¬ 
tions  at  a  peak  power  density  of  170  mW/cm2  (averaged  0.51  and  5.1  mWcm2).  No  work  stoppage 
occurred  to  pulsed  radiation.  Taken  in  sum,  the  data  show  that  the  mature  Long-F.vans  rat  is  resonant 
at  a  frequency  near  600  MHz  while  work  stoppage  during  short-term  exposures  to  600-MHz  radiation 
occurs  at  a  power  density  between  7.5  and  10  mW/cm2. 


1.  INTRODUCTION 

Biological  effects  of  microwave  radiation  have  gen¬ 
erated  considerable  interest  in  recent  years.  Much  of  this 
interest  has  focused  on  variables  that  control  the  amount 
of  energy  absorbed  by  the  living  animal.  Recent  investi- 
tations  have  shown  that  such  variables  as  the  animal’s 
mass  and  geometry  and  the  wavelength  of  radiation  are 
potent  determiners  in  production  of  differential  absorp¬ 
tion  of  microwave  energy  [Justesen  and  King,  1970; 
Gandhi,  1974;  Githens  etal.,  \915,Schrot  and  Hawkins, 
1976;  D’Andrea  et  al.,  1976] .  The  peak  of  a  differential 
curve  of  microwave  absorption  has  been  termed  reso¬ 
nance.  A  clear  demonstration  of  resonance  would  re¬ 
quire  nonmonotonic  rates  of  absorption  across  a  range 
of  frequencies,  the  maximum  occurring,  by  definition, 
at  the  resonant  frequency.  Recently,  D’Andrea  et  al. 
[1976]  demonstrated  an  increasingly  larger  rise  of 
body  temperature  and  an  increasing  tendency  for 
behavioral  disruption  to  occur  in  animals  as  frequency  of 
radiation  was  increased  from  200  to  500  MHz.  However, 
500  MHz  is  below  the  resonant  frequency  of  the  rat. 
The  present  study  is  part  of  a  program  to  delineate  a 
predicted  disruption  of  behavior  at  a  frequency  near 
600  MHz,  which  approximates  the  resonant  wavelength 
of  the  long  axis  of  the  mature  rat. 

2.  METHOD  AND  MATERIALS 

2.1.  Rationale. 

Eleven  rats  were  trained  to  perform  a  lever-pressing 
task.  The  task  was  performed  in  a  radiation  chamber  for 


a  reward  -  a  Noyes  food  pellet  -  on  a  variable  interval 
(30-second  average)  schedule  of  reinforcement.  Over 
many  sessions  one  group  of  five  rats  was  exposed  in 
random  order  to  400-,  500-,  600-,  and  700-MHz  CW 
radiation  at  a  power  density  of  20  mW/cm2 .  In  addition, 
the  remaining  rats  were  exposed  in  random  order  to 
600-MHz  CW  radiation  at  20,  10,  7.5,  and  5  mW/cm’ 
to  determine  latencies  of  work  stoppage  Subsequently, 
three  animals  of  this  group  were  exposed  to  pulsed 
600-MHz  radiation  at  3-  and  30-microsecond  pulse 
durations  and  170-mW/cm2  peak  power  density. 

2.2.  Subjects. 

Eleven  male  Long-Evans  rats  were  obtained  from 
Simonsen  of  California  and,  at  commencement  of  study, 
had  body  masses  ranging  from  420  to  450  grams.  All  of 
the  rats  had  been  partially  deprived  of  food  until  their 
averaged  body  mass  stabilizated  at  85  percent  of  that 
before  deprivation.  All  of  the  rats  were  maintained  at 
this  level  of  Noyes  45-mg  pellets  and  carefully  metered 
supplemental  feeding  one  hour  after  each  session  with 
Simonsen  4.5  rat  chow.  Each  animal  was  housed  separ¬ 
ately  in  a  35  by  30  by  16  cm  cage.  Each  rat  was  given 
free  access  to  water  in  the  home  cage  and  was  kept  on 
a  1 2/ 12-hour  light/dark  cycle  (on,  0700  hours;  off, 
1900  hours).  Ambient  temperature  and  relative  humid¬ 
ity  in  the  colony  were  maintained  at  23  ±  1°C  and  10 
to  30%  respectively.  The  averaged  length  from  snout  to 
base  of  tail  of  the  animals  was  19  cm. 

2.3.  Apparatus. 

The  microwave  radiation  chamber  (Figure  1)  consisted 
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Fig.  1.  Monopole-above-ground  radiation  chamber. 


of  a  Filshield  grounded  Faraday  cage  (3  by  3  by  3 
meters).  The  interior  of  the  cage  was  partially  lined 
with  Eccosorb  AN-79  material.  One  wall  and  one  half 
of  the  floor  and  ceiling  of  the  chamber  opposite  a  radiat¬ 
ing  antenna  were  covered  with  Eccosorb  VHP- 18  pyra¬ 
midal  material.  The  remaining  wall  was  covered  by  a 
2-mm  thick  copper  sheet  to  form  a  ground  plane. 
A  single  brass  monpole  antenna  (1/4  X  length,  1.2-cm 
diameter)  was  inserted  through  the  copper  sheet.  A  cop¬ 
per  sheet  (3  m  by  0.5  m  by  2  mm  thick)  was  formed 
into  a  45°  reflector  and  was  mounted  on  the  copper  wall 
33  cm  behind  the  antenna.  An  animal-response  cage 
(Figure  2)  that  was  constructed  of  Plexiglas  was  located 
143  cm  from  the  antenna.  The  cage  was  held  on  the 
vertical  wall  by  a  Styrofoam  support  that  was  con¬ 
structed  of  2.5-cm  thick  Styrofoam  sheet. 

The  chamber  was  equipped  with  a  50-watt  spotlight 
and  a  sensitive  RCA  television  camera  that  was  coupled 
to  an  external  video  monitor.  A  ventilation  fan  was 
mounted  on  top  of  the  chamber  and  circulated  ait  from 
the  room  through  the  chamber  at  a  rate  of  approxi¬ 
mately  eight  cubic  meters  per  minute. 

Since  standard  metallic  caging  and  behavioral  operands 
cannot  be  used  in  a  microwave  field,  an  alternate  system 
of  monitoring  and  reinforcing  lever-pressing  behavior 
within  the  radiation  chamber  was  devised.  A  plastic 
response  cage  (Figure  2)  was  constructed  using  two, 
24  by  16  cm,  2-cm  thick  Plexiglas  panels.  The  panels 
were  joined  at  top  and  bottom  by  32  rods  (5-mm 
dia  by  13.3-cm  long)  of  Plexiglas  equally  spaced  to 
form  a  grid  floor  and  top.  Two  pieces  of  Plexiglas  sheet 
(12.7  by  15.2  cm)  were  used,  one  as  an  end  panel  and 
one  as  a  sliding  guillotine  door  for  the  response  cage. 
A  plastic  lever  was  mounted  on  a  Plexiglas  rod  (5  mm  by 
13  cm)  at  one  end  of  the  cage.  One  end  of  the  rod  was 


cemented  to  the  cage,  while  the  other  end  was  free  to 
move.  Depression  of  the  lever  was  detected  remotely  by 
its  interrruption  of  a  beam  of  light  that  normally  passed 
via  Dupont  Crofon  fiber  optic  material  to  and  from 
the  plastic  response  cage.  A  sensitive  photocell  relay  in 
conjunction  with  timing  circuitry,  digital  printing 
counters,  and  cumulative-response  recorders  served  to 
detect,  to  program  and  to  record  events  of  the  experi¬ 
ments.  A  BRS/LVE  pellet  feeder  was  mounted  outside 
the  chamber.  A  Noyes  food  pellet  (45  mg)  was  delivered, 
when  programmed,  to  the  response  cage  by  a  tube  of 
flexible  Tygon  plastic  (1.27  cm  dia  by  1.5  m  length). 

The  microwave  source  used  was  an  Epsco  tunable 
signal  generator  (200  to  750  MHz)  with  a  minimum 
power  output  of  100  watts.  For  the  exposures  to  pulsed 
radiation,  an  MCL-1114  one-kilowatt  linear  amplifier 
(400  to  800  MHz)  was  used  in  conjunction  with  the 
Epsco  source.  A  Hewlett-Packard  No.  214A  pulse 
generator  was  used  to  drive  the  signal  source  at  1000 
pulses  per  second.  Radiation  frequency  was  monitored 
by  a  DANA  7580  frequency  counter.  A  Sierra  Philco 
168B  power  meter  was  used  to  monitor  the  forward 
and  reverse  powers  at  the  monopole  antenna.  A 
Weinschell  double-stub  tuner  was  used  to  match  tire 
signal  source  to  an  antenna  (1/4  X  for  each  frequency 
used)  and  to  reduce  reflected  power  to  a  minimum.  The 
power  density  of  the  microwave  field  within  the  plastic 
response  cage  in  the  radiation  chamber  at  the  position 
normally  occupied  by  the  rat  was  measured  in  mW/cm5 
by  a  General  Microwave  Corporation  Raham  Model  1 
(300  to  18,000  MHz)  probe.  The  probe  was  inserted 
into  the  radiation  chamber  from  the  top  with  the 
length  of  the  probe  perpendicular  to  the  £'-field  vector 
and  was  remotely  monitored  outside  of  the  radiation 
chamber  before  and  after  each  exposure  of  a  rat. 

2.4.  Procedure. 

All  rats  were  trained  to  press  the  lever  for  food  pellets 
on  a  variable  interval  schedule  of  reinforcement.  The 
clock  time  used  for  the  schedule  was  three  seconds,  with 
a  probability  of  0.10  of  delivery  of  a  food  pellet.  With 
this  arrangement  a  food  pellet  could  be  delivered  at  in- 
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tervals  of  three  seconds  or  longer,  the  averaged  time  of 
delivery  being  30  seconds.  Each  rat  was  scheduled  daily 
for  a  one-hour  session  of  lever  pressing,  until  a  stable 
rate  of  responding  was  achieved.  Lever-pressing  behavior 
was  considered  stable  when  the  total  number  of  lever 
presses  during  a  session  differed  by  less  than  15%  from 
the  total  of  a  previous  session.  On  the  average,  approxi¬ 
mately  14  sessions  were  required  to  achieve  stable 
responding. 

Once  stable  responding  was  observed  the  animals  were 
randomly  divided  into  two  groups.  Each  of  the  animals 
of  the  first  group  ( n  =  5)  was  individually  exposed  in 
random  order  to  400-,  500-,  600-,  and  700-MHz  CVV 
radiation  at  a  power  density  of  20  mW/cm2  while 
responding  in  the  radiation  chamber.  Each  of  the  ani¬ 
mals  of  the  second  group  (n  =  6),  also  while  responding 
in  the  radiation  chamber,  was  individually  exposed  in 
random  order  to  CW  radiation  at  600  MHz  at  20,  10, 
7.5,  and  5  mW/cm2 .  In  addition,  three  of  the  animals  in 
the  second  group  were  subsequently  exposed  to  600- 
MHz  pulsed  radiation  (1000  pulses  per  second  at  3- 
and  at  30-ps  pulse  durations;  170  mW/cm2  peak  power 
density,  0.51  and  5.10  mW/cm2  averaged.  The  long  axis 
of  the  rat’s  body  in  each  exposure  was  parallel  to  the 
vector  of  the  electric  field. 

Four  to  six  retraining  sessions  were  conducted  for  each 
rat  after  each  exposure  to  microwaves  to  maintain 
stable  rates  of  responding  and  to  check  for  possible 
carry-over  of  behavioral  effects  from  sessions  of  prior 
exposure.  Exposures  to  microwaves  began  after  the  fifth 
minute  of  each  session  and  were  terminated  at  the  end 
of  the  first  one-minute  period  during  which  responding 
fell  below  one-third  of  the  rat’s  normal  rate  of  respond¬ 
ing,  as  determined  by  data  from  the  previous  daily  ses¬ 
sion.  If  a  rat  responded  consistently  throughout  the 
session,  as  defined,  irradiation  was  continuous  for  55 
minutes.  The  primary  measure  was  time  to  work  stop¬ 
page  after  onset  of  irradiation  until  the  60th  second  of 
the  minute  during  which  responding  fell  below  the 
33%  criterion.  Each  rat  was  removed  fiom  the  response 
cage  and  radiation  chamber  at  stoppage  of  work  or  at 
session’s  end,  whichever  occurred  first. 

Measurements  of  colonic  temperature  were  made  with 
a  BAT-8  digital  electronic  thermometer,  first  at  the 
beginning  of  a  radiation  session  and  again  at  the  time  of 
work  stoppage  (or  at  sessions’s  end)  to  determine  the 
amount  of  microwave-induced  heating.  Measurement 
of  colonic  temperature  was  also  made  before  and  after 
several  training  and  retraining  sessions  on  a  random 
basis.  Ambient  temperature  of  the  radiation  chamber 
was  controlled  by  room  temperature  (21°  to  22  °C) 
as  monitored  on  the  BAT-8  thermometer.  Relative 
humidity  of  the  radiation  chamber  and  surrounding 
room  was  monitored  by  a  Mason  hygrometer  during 
radiation  sessions  (27%  mean,  2.7%  SE).  During  all 
training  and  radiation  sessions,  each  animal  was  moni¬ 
tored  on  the  closed-circuit  television. 

3.  RESULTS 

Frequency  of  microwave  radiation  proved  to  be  an 
important  variable  in  heating  of  the  rat  and  consequent 


stoppage  of  work.  For  the  rats  exposed  to  400-,  500-, 
600-,  and  700-MHz  CW  radiation  at  20  mW/cm2,  the 
most  rapid  stoppage  of  work  (Figure  3)  and  the  highest 
elevations  of  body  temperature  (Figure  4)  occurred  at 
600  MHz.  Averages  of  both  the  rise  in  temperature  and 
the  time  to  work  stoppage  reflect  a  nonmonotonic  rela¬ 
tion  to  frequency  of  radiation.  In  other  words,  the  rate 
of  absorption  of  energy  increased  as  did  the  likelihood 
of  behavioral  disruption  as  frequency  increased  to  600 
MHz  and  then  declined  at  the  highest  frequency  of  700 
MHz. 

Irrespective  of  frequency,  animals  exposed  to  radiation 
at  20  mW/cm2  exhibited  evidence  of  vasodilation,  be¬ 
came  immobilized,  and  were  prone  upon  removal 
from  the  radiation  chamber.  The  animals  were  clear¬ 
ly  stressed  by  heat.  However,  as  reported  earlier 
[Justesen  and  King,  1970;  D’Andrea  et  al.,  1976],  the 
animals  quickly  regained  an  upright  posture  and  mobil¬ 
ity  within  minutes  after  removal  from  the  radiation 
chamber.  At  the  time  of  work  stoppage,  nearly  all 
animals  engaged  in  licking  behavior,  possibly  to  induce 
evaporative  cooling,  which  continued  after  the  animal 
was  removed  from  the  radiation  chamber  and  was 
placed  in  its  home  cage. 

A  repeated-measures  analysis  of  variance  [Kirk,  1968] 
was  used  to  evaluate  the  effects  of  the  treatments.  The 
analysis  revealed  a  significant  difference  among  the 
radiation  frequencies  with  respect  to  time  to  work 
stoppage  (F  —  16.41,  df  =  3/12, /*<  0.01).  A  Newman- 
Keuls  test  on  the  ordered  means  indicated  a  significant 
difference  in  times  to  work  stoppage  between  600  MHz 
and  500  MHz  (P  <  0.05),  between  600  MHz  and  700 
MHz  (P  <  0.05),  and  between  400  MHz  and  all  other 
frequencies  (Ps  <  0.01).  The  difference  between 
exposures  at  500  and  700  MHz  was  not  significant 
(P>0.05). 

As  an  indicant  of  energy  dosing,  the  rise  of  body  tem¬ 
perature  of  the  exposed  animal  was  divided  by  the 
duration  of  exposure  in  minutes.  Moderate  rises  of 
body  temperatures  were  also  observed  during  training 
sessions.  These  data  are  presented  in  Figure  4  with  the 
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Fig.  3.  Mean  latencies  (±  SE)  of  time  to  work  stoppage  as  a 
function  of  radiation  frequency  at  20  mW/cm*  power  density. 
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Fig.  4.  Rise  of  colonic  temperature  (±  SE)  per  unit  of  time  as 
a  function  of  frequency  of  radiation  at  20-mW/cm2  power 
density.  Shaded  area  represents  averaged  increment  of  body 
temperature  (±  SE)  of  rats  that  responded  during  training 
sessions  in  the  absence  of  radiation. 


shaded  area  representing  averaged  increases  of  body 
temperature  during  training.  A  repeated-measures 
analysis  of  variance  revealed  significant  differences 
across  radiation  frequencies  with  respect  to  heating 
(F  =  6.06,  df  —  3/12,  P  <  0.01).  A  Newman-Keuls  test 
of  ordered  means  revealed  significant  differences  of 
temperature  increments  for  all  possible  comparisons 
(ft  <0.01). 

The  measures  of  time  to  work  stoppage  and  body 
heating  followed  each  other  very  closely  and,  of  course, 
inversely.  As  body  temperature  increased,  time  to 
work  stoppage  became  shorter.  The  coefficient  of  the 
product-moment  correlation  between  time  to  work 
stoppage  and  rise  of  body  temperature  for  all  radia¬ 
tion  exposures  is  —  0 .78  (at  df  =  1 8 ,  P  <  0 .00 1 ). 

A  comparison  was  made  between  the  total  of  responses 
during  a  session  before  and  that  after  exposure  of 
the  rats  to  radiation.  This  was  done  to  determine  if 
a  session  of  radiation  influenced  the  animals’  responding 
during  the  subsequent  session  of  retraining,  as  compared 
to  responding  during  the  training  session  prior  to  expo¬ 
sure.  A  repeated-measures  factorial  analysis  of  variance 
[Kirk,  1968)  revealed  no  differences  between  the 
average  of  total  responses  for  the  session  before  and  that 
after  the  exposure  to  radiation  (F  =  0.01,  df  —  1/28, 
P  >  0.05). 

A  comparison  was  also  made  between  the  time  to  work 
stoppage  for  each  rat  as  a  function  of  the  order  in  which 
treatments  were  given.  This  was  done  to  determine  if 
the  order  (which  was  putatively  random)  of  the  radiation 
sessions  may  have  introduced  an  avoidance  contingency, 
i.e.,  may  have  resulted  in  a  rat’s  learning  to  stop  working 
in  order  to  avoid  radiation.  An  analysis  of  variance 
[Kirk,  1968]  revealed  but  chance  differences  between 
the  variable  of  ordering  and  that  of  work-stoppage 
latency  (F=  0.30,  df=  3/16,  P>  0.05). 

For  the  rats  exposed  to  600-MHz  CW  radiation  at 


power  densities  of  20  and  10  mW/cml,  stoppage  of 
work  invariably  occurred.  As  shown  in  Figure  5,  the 
averaged  time  to  work  stoppage  and  variability  (i.e., 
magnitude  of  standard  errors)  both  tended  to  increase 
as  power  density  decreased.  The  earliest  work  stoppage 
was  observed  in  animals  exposed  to  radiation  at  20 
mW/cm2.  Exposures  at  10  mW/cm2  required  nearly 
twice  the  averaged  duration  to  reduce  the  animals’ 
response  rate  to  the  work-stoppage  criterion.  Exposures 
to  5  and  7.5  mW/cm2  produced  little  evidence  of  work 
stoppage.  A  repeated-measures  analysis  of  variance 
[Kirk,  1968]  was  used  to  evaluate  the  work-stoppage 
variable  as  a  function  of  power  density.  The  analysis 
revealed  a  highly  significant  difference  (F=  35.29, df = 
3/15,  P  <  0.01).  A  Newman-Kculs  test  on  the  ordered 
means  of  time  to  work  stoppage  revealed  a  significant 
difference  in  time  to  work  stoppage  between  20 
mW/cm2  and  all  other  exposures  (ft  <  0.01),  between 
10  mW/cm2  and  7.5  mW/cm2  (P<  0.05),  and  between 
10  mW/cm2  and  5  m\V/cm2  (P<  0.05).  There  was  little 
difference  between  exposures  at  7 5  mW/cm2  and 
5  mW/cm2  (F>0.05). 

The  body  heating  induced  by  the  exposures  to  600- 
MHz  radiation  at  different  power  densities  (Figure  6) 
also  followed  a  nearly  monotonic  function.  The  temper¬ 
ature  differential  was  greatest  at  20  mW/cm2  and 
decreased  by  nearly  one  half  at  10  mW/cm2.  Body 
temperatures  at  5  and  7.5  mW/cm2  were  only  slightly 
higher  than  those  of  baseline  measures  (shaded  area, 
Figure  6).  A  repeated-measures  analysis  of  variance 
revealed  a  significant  difference  in  body  temperature  as  a 
function  of  power  density  (F  =  5.89,  df  =  3/15, 
P  <  0.01).  A  Newman-Keuls  test  on  the  ordered  means 
of  body  temperature  revealed  a  significant  difference 
in  increase  of  body  temperature  as  a  function  of  power 
densities  between  20  mW/cm2  and  lower  power  densities 
of  {P  <  0.05).  All  other  comparisons  were  not  significant . 
The  coefficient  of  the  product-moment  correlation 
between  times  to  work  stoppage  and  temperatures  across 
power  densities  is  -  0.83  (at  df  =  22,  F<  0-01). 

A  comparison  was  again  made  between  the  total 
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Fig.  5.  Mean  latency  (±  SF.)  of  time  to  work  stoppage  as  a 
function  of  power  density  at  600  mHz. 
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Fig.  6.  Rise  of  cotonic  temperature  (±  SE)  per  unit  of  time 
as  a  function  of  power  density.  Frequency  of  radiation  =  600 
MHz.  Shaded  area  represents  average  and  increment  of  body 
temperatures,  (±  SE)  of  rats  that  responded  during  training 
sessions  in  the  absence  of  radiation. 


number  of  responses  that  were  recorded  during  the  ses¬ 
sion  before  and  that  after  exposures  to  radiation  to 
determine  if  a  session  of  radiation  influenced  the  ani¬ 
mals’  responding  during  subsequent  retraining.  A 
repeated-measures  factorial  analysis  of  variance  [Kirk, 
1968]  revealed  no  differences  between  the  average  of 
total  responses  for  a  session  before  and  that  after  the 
radiation  sessions  (F  =  0.07,  df  —  1/35,  P  >  0.05/. 

There  was  no  reliable  effect  of  pulsed  radiation  (1000 
pps,  170  mW/cm2  peak  power)  on  suppression  of  lever¬ 
pressing  behavior  at  either  3-  or  30- microsecond  pulse 
durations.  All  three  of  the  animals  exposed  under  these 
conditions  responded  at  near  normal  rates  throughout 
the  radiation  session.  One  animal  stopped  responding 
near  the  end  of  a  session  of  pulsed  radiation  at  3  micro¬ 
seconds.  The  data  on  body  heating  for  these  exposures, 
however,  indicate  but  little  elevation  of  body  tempera¬ 
tures. 


4.  DISCUSSION 

Our  study  was  conducted,  first,  to  investigate  the  con¬ 
tribution  of  radiation  frequency  on  the  rate  of  heating 
of  the  living,  active  rat  and,  second,  to  investigate  the 
effects  that  such  radiation  may  have  in  altering  a  well- 
trained  behavior.  Considering  the  four  frequencies  of 
radiation  that  were  studied,  our  data  indicate  that  600 
MHz  approximates  the  resonant  frequency  of  the  adult 
rate  with  a  body  mass  of  380-400  grams  and  averaged 
length  (less  the  tail)  of  19  to  20  cm  for  the  long  axis  -- 
and  with  this  long  axis  oriented  parallel  to  the  vector  of 
the  electric  field.  The  proximity  to  resonance  is  reflected 
both  by  the  time  to  work  stoppage  and  by  measures 
of  body  heating,  both  of  which  yielded  nonmonotonic 
functions.  Our  findings  of  resonance  near  600  MHz  are 


in  good  agreement  with  previous  work  on  the  calculated 
resonant  frequency  of  this  animal  [Gandhi,  1975; 
Johnson  et  al,  1976].  It  was  of  some  interest  to  us 
that  the  tail  of  the  rat  did  not  seem  to  contribute  to  the 
factor  of  resonance. 

Our  use  of  the  rise  of  colonic  temperature  (AT)  in  the 
behaving  animal  was  only  as  an  indicant  of  energy  dos¬ 
ing.  Not  only  will  the  intact  mammal  dissipate  thermal 
energy,  emotional  reactions  to  the  environment  can  pro¬ 
duce  psychogenic  “fever”  [cf .  Justesen  et  al,  1974  with 
Lu  et  al.,  1977].  We  have  previously  observed  (unpub¬ 
lished)  that  the  active  alert  animal  as  compared  with 
the  anesthetized  animal  has  a  much  higher  body  temper¬ 
ature.  The  indexing  of  absorbed  microwave  energy  by 
body  temperature,  as  described  in  this  report,  should 
only  be  construed,  therefore,  as  a  relative  measure. 

Radiation  at  a  frequency  of  600  MHz,  given  a  55- 
minute  maximum  period  of  radiation,  produced  disrup¬ 
tion  of  lever- pressing  behavior.  As  measured  by  the 
work-stoppage  criterion,  disruption  occurred  at  power 
densities  of  20  and  10  raW/cm' .  Disruption  of  behavior 
was  not  evident  at  lower  power  densities  at  7.5  and 
5  mW/cm2 .  However,  our  data  do  not  infirm  the  possi¬ 
bility  that  biological  effects  can  occur  at  power  densi¬ 
ties  below  1 0  mW/cm2  . 

The  use  of  pulsed  energy,  albeit  at  an  averaged  power 
density  of  5  mW/cm2 ,  showed  no  reliable  effect  on  time 
to  work  stoppage.  This  is  not  surprising  since  CW  radia¬ 
tion  at  the  same  power  density  did  not  disrupt  behavior. 
The  “package”  (dose)  of  energy  per  radiated  pulse  was 
about  half  that  found  by  Guy  et  al  [1975]  to  be 
necessary  for  threshold  stimulation  of  an  acoustic 
response.  Given  full  resonance,  which  may  not  have 
been  observed  by  us,  and  given  the  use  of  higher  peak 
values  of  power  density  than  used  by  us,  it  is  possible 
that  behavioral  disruption  will  occur  to  pulsed  radiation 
at  lower  averaged  power  densities.  Too,  the  greater  sen¬ 
sitivity  of  other  behavioral  tests  -  such  as  conditional 
suppression  -  is  much  more  likely  to  reveal  effects  of 
radiation  at  low  power  densities  [King  et  al.,  1971], 

As  noted  earlier  [D’Andrea  et  al.,  1976],  no  reliable 
carry-over  of  behavioral  effects  was  observed  from  one 
radiation  session  to  the  next.  Visual  observation  of  the 
animals  during  responding  in  subsequent  retraining 
sessions  did  not  indicate  any  attempts  to  escape  from 
the  response  cage,  as  one  might  see  with  rats  upon  their 
return  to  an  operant  chamber  where  electric  foot  shock 
has  previously  been  administered.  The  lack  of  carry¬ 
over  suggests  that  the  most  reliable  effect  of  the  radia¬ 
tion  was  sufficient  heating  of  the  animal  to  suppress 
behavior  but  the  radiation  was  not  intense  enough  or 
of  sufficient  duration  to  produce  noticeable  physical 
damage. 
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